Chapter 1—Introduction 1-1 


SOLID STATE TESLA COIL 
by 


Dr. Gary L. Johnson 
Canon City, Colorado 


Some years ago I developed an interest in Tesla coils. I was teaching a senior elective 
course at Kansas State University where we talked about power MOSFETs and topics related 
to high voltages and currents. I decided to use a Tesla coil as a class project. We would talk 
about design aspects, then design, build, and test a coil. The best description of the results of 
that plan was fiasco, or maybe disaster. We had some sparks, but none where they belonged. 
That was one of the most humiliating experiences of my career. 


I learned several things from that experience. One is that the Tesla coil is more complex 
than I had thought. Another was that there seemed to be a mismatch between theory and ex- 
periment. At that time, at least, people would go through pages of high powered mathematics 
and quit without giving an example of how to use all the formulas. Experimentalists would 
sometimes make fun of the theorists, and give rules-of-thumb on how to make long sparks. 
It was like I was hearing a debate on whether the best cooks use recipes or not. My mother 
never used a recipe and I always enjoyed her cooking. However, my own talents are such that 
if Iam to cook anything fit to eat, I need a recipe. 


This book is written for people like me, challenged when faced with doing something 
without a recipe or complete set of instructions. I will throw in things learned from other 
Tesla coil builders, but will quickly admit that when it comes to making long sparks, there 
are many who are far better than I. 


I started asking questions about Tesla coils that any electrical engineer would ask. These 
include: 


1. What is the input impedance? 


2. What are the fractions of input power that are dissipated in the spark itself, in elec- 
tromagnetic radiation, the coil wire, the coil form, the toroid, the spark gap, and other 
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circuit components? 


3. Are there circuit models that allow these questions to be answered on the computer 
before building and testing devices in open air? 


4. What are the differences between Tesla coils driven by or through spark gaps, vacuum 
tubes, or solid state devices? 


5. What are the important factors in producing long sparks (energy per bang, power input 
at spark inception, rate of change of power, the coil, the toroid, etc.)? 


One would expect the answers to these questions to come from a mix of theory and 
experiment. One would develop a theory or model and then go to the laboratory to measure 
parameters and check performance. The theory would then be adjusted to reflect experimental 
observations. 


We now review a little Tesla coil history and look at the ‘simplest’ model, the lumped 
circuit element model. 


1.1 History 


Nikola Tesla (1856 - 1943) was one of the most important inventors in human history. He 
had 112 U.S. patents and a similar number of patents outside the United States, including 30 
in Germany, 14 in Australia, 13 in France, and 11 in Italy. He held patents in 23 countries, 
including Cuba, India, Japan, Mexico, Rhodesia, and Transvaal. He invented the induction 
motor and our present system of three-phase power in 1888 [20]. He invented the Tesla coil, 
a resonant air-core transformer, in 1891. Then in 1893, he invented a system of wireless 
transmission of intelligence. Although Marconi is commonly credited with the invention of 
radio, the U.S. Supreme Court decided in 1943 that the Tesla Oscillator patented in 1900 
had priority over Marconi’s patent which had been issued in 1904 [15]. Therefore Tesla 
did the fundamental work in both power and communications, the major areas of electrical 
engineering. These inventions have truly changed the course of human history. 


After Tesla had invented three-phase power systems and wireless radio, he turned his 
attention to further development of the Tesla coil. He built a large laboratory in Colorado 
Springs in 1899 for this purpose. The Tesla secondary was about 51 feet in diameter. It was 
in a wooden building in which no ferrous metals were used in construction [15]. There was 
a massive 80-foot wooden tower, topped by a 200-foot mast on which perched a large copper 
ball which he used as a transmitting antenna. The coil worked well. There are claims of 
bolts of artificial lightning over a hundred feet long, although Richard Hull asserts that from 
Tesla’s notes, he never claimed a distance greater than 43 feet. From photographic evidence, 
the maximum may have been closer to 22 feet [12]. 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 1—Introduction 1-3 


Tesla then abandoned the Colorado Springs Laboratory early in 1900, having learned 
what he needed from that facility, and also having become somewhat unpopular as a result 
of frequently knocking the local sub-station off line. 


Since that time, it appears that no one has built a Tesla coil of both the size and perfor- 
mance of the Colorado Springs coil. Apparently the only coil of that size was built by Robert 
Golka at Wendover Air Force Base in Utah [8] and later moved to a facility near Leadville, 
Colorado [9, 19]. The original purpose of this coil was to produce artificial lightning for 
testing the effects of lightning striking aircraft in flight. Golka determined that the average 
voltage produced in Utah was about 10 MV, with the highest voltage observed being 25 MV. 
Operation was spectacular, even if not quite at the level of the Colorado Springs coil. 


When Golka’s coil was moved to Leadville, however, it performed very poorly. Golka and 
his associates were basically unable to properly tune the coil. There has been considerable 
speculation over the reasons for the difference in performance, but one problem seems to 
be that we did not have adequate theoretical models for the design and operation of Tesla 
coils. What appeared to be minor differences in location and construction caused a major 
decrease in performance. The number of variables was simply too large to allow for a purely 
experimental optimization of performance before the coil was dismantled and moved early in 
1990. 


Some work on theoretical models has been performed by high energy physicists [6, 5, 
1, 17, 18]. They are interested in high voltage capacitor discharges for research in plasma 
physics and in the production of pulsed particle or radiation beams. The most common way 
of producing such high voltage discharges is the Marx circuit, in which capacitors are charged 
in parallel to a lower voltage and then discharged in series through a number of airgaps. The 
Marx circuit requires the capacitor bank to be divided into sub-banks well-insulated from each 
other and from ground. A Tesla coil offers an alternative method of charging the high voltage 
capacitors. Discharges are reported in the range of 100 kA at 1 MV, with one report of 2.5 
MV [5]. These models are all lumped parameter models. 


There are a number of experimenters who build Tesla coils as a hobby. The Tesla Coil 
Builders Association has several hundred members and a quarterly newsletter published by 
Harry Goldman [7]. Harry has announced plans to stop publishing the newsletter at the end of 
2001. The Tesla Coil Builders of Richmond has been a very active local group [11], although 
their leader Richard Hull has recently become interested in other activities. A number of 
manuals are available on how to build coils [16, 4, 5]. The one by Lee [16] is especially 
well illustrated with pictures of capacitors and other components that might be needed for a 
moderate sized Tesla coil. There is an Internet listserv (www.pupman.com) that has about 
700 subscribers, which has been very helpful to me. 


The International Tesla Society was formed in Colorado Springs, Colorado, around 1984. 
They held a large conference every other year, covering the full range of concepts related to 
Nikola Tesla, including the Tesla coil. Bill Wysock would often attend and bring his large 
coils for display. One year someone set up a smallish coil in the hotel ballroom. It did a great 
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job of overriding the elevator controls and the fire alarm system. The next year we were at a 
different hotel (for reasons unknown to me) and the Tesla coils were out at the far side of the 
parking lot. The conferences were a lot of fun. I miss them since the organization folded in 
the late 1990s. 


The brothers James and Kenneth Corum have done considerable work on distributed 
models of Tesla coils in the past few years [2, 3]. They argue that lumped parameter models 
are not adequate for all situations. Sometimes a distributed circuit analysis must be made. In 
this case, the Tesla coil secondary and another component called the extra coil are considered 
as sections of transmission lines. This explains some of the effects in an elegant manner. They 
have written a sophisticated computer program, TCTUTOR, to analyze Tesla coils. They 
have also performed considerable historical research into Tesla’s notes made at his facility in 
Colorado Springs [21]. 


The Tesla coil community is divided over the issue of lumped versus distributed models. 
A majority favors the lumped model approach. Some are outspoken in their belief that 
distributed models are useless at best and just plain wrong on important issues. I confess to 
being somewhere in the middle on this controversy. James Corum and I both have our Ph.D.s 
in electromagnetic theory, so I can mostly understand what he says, and I therefore have a 
natural orientation to the distributed approach. In my eyes, I am like a Baptist pastor of a 
50 person congregation and James is like Billy Graham. That is, I hold him in awe. I have 
heard the Corums speak several times, and have gotten caught up in their knowledge and 
excitement. 


On the other hand, I cannot honestly say that TCTUTOR has been helpful to me in 
building and understanding Tesla coils. I can see significant problems with distributed models, 
which will be discussed later. And James, like many bright people, has a tendency to talk 
down to us slow ones. This puts some people off, of course. 


In this book we will look at both lumped and distributed models. We will point out 
difficulties with both. We will look at some data, and ask which approach does best in 
describing reality. 


1.2 Classical Tesla Coil 


A classical Tesla coil contains two stages of voltage increase. The first is a conventional iron 
core transformer that steps up the available line voltage to a voltage in the range of 12 to 50 
kV, 60 Hz. The second is a resonant air core transformer (the Tesla coil itself) which steps up 
the voltage to the range of 200 kV to 1 MV. The high voltage output is at a frequency much 
higher than 60 Hz, perhaps 500 kHz for the small units and 80 kHz (or less) for the very large 
units. 


The lumped circuit model for the classical Tesla coil is shown in Fig. 1.1. The primary 
capacitor C is a low loss ac capacitor, rated at perhaps 20 kV, and often made from mica or 
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polyethylene. The primary coil D1 is usually made of 4 to 15 turns for the small coils and 1 to 
5 turns for the large coils. The secondary coil [2 consists of perhaps 50 to 400 turns for the 
large coils and as many as 400 to 1000 turns for the small coils. The secondary capacitance 
C2 is not a discrete commercial capacitor but rather is the distributed capacitance between 
the windings of Lz and the voltage grading structure at the top of the coil (a toroid or sphere) 
and ground. This capacitance changes with the volume charge density around the secondary, 
increasing somewhat when the sparks start. It also changes with the surroundings of the coil, 
increasing as the coil is moved closer to a metal wall. This may have been one of the reasons 
that Golka’s coil worked better in Utah than in Colorado, because the metal walls were closer 
to the coil in Colorado. 


Cl) 
Va Ub < 03 Lo C2 


iron core air core 


Figure 1.1: The Classical Tesla Coil 


The symbol G represents a spark gap, a device which will arc over at a sufficiently high 
voltage. The simplest version is just two metal spheres in air, separated by a small air gap. 
It acts as a voltage controlled switch in this circuit. The open circuit impedance of the gap 
is very high. The impedance during conduction depends on the geometry of the gap and the 
type of gas (usually air), and is a nonlinear function of the current density. This impedance 
is not negligible. A considerable fraction of the total input power goes into the production of 
light, heat, and chemical products at the spark gap. In any complete analysis for efficiency, an 
equivalent gap resistance Ryaqp could be defined such that 1 Regis would represent the power loss 
in the gap. This would have rather limited usefulness because of the mathematical difficulty 
of describing the arc. 


The arc in the spark gap is similar to that of an electric arc welder in visual intensity. 
That is, one should not stare at the arc because of possible damage to the eyes. At most 
displays of classical Tesla coils, the spark gap makes more noise and produces more light than 
the electrical display at the top of the coil. 


When the gap is not conducting, the capacitor C) is being charged in the circuit shown 
in Fig. 1.2, where just the central part of Fig. 1.1 is shown. The inductive reactance is much 
smaller than the capacitive reactance at 60 Hz, so L; appears as a short at 60 Hz and the 
capacitor is being charged by the iron core transformer secondary. 
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Figure 1.2: C; Being Charged With The Gap Open 


A common type of iron core transformer used for small Tesla coils is the neon sign trans- 
former (NST). Secondary ratings are typically 9, 12, or 15 kV and 30 or 60 mA. An NST has 
a large number of turns on the secondary and a very high inductance. This inductance will 
limit the current into a short circuit at about the rated value. An operating neon sign has a 
low impedance, so current limiting is important to long transformer life. However, in Tesla 
coil use, the NST inductance will resonate with C;. The NST may supply two or three times 
its rated current in this application. Overloading the NST produces longer sparks, but may 
also cause premature failure. 


When the voltage across the capacitor and gap reaches a given value, the gap arcs over, 
resulting in the circuit in Fig. 1.3. We are not interested in efficiency in this introduction so 
we will model the arc as a short circuit. The shorted gap splits the circuit into two halves, 
with the iron core transformer operating at 60 Hz and the circuit to the right of the gap 
operating at a frequency (or frequencies) determined by C;, £1, L2, and C2. It should be 
noted that the output voltage of the iron core transformer drops to (approximately) zero 
while the input voltage remains the same, as long as the arc exists. The current through the 
transformer is limited by the transformer equivalent series impedance shown as R, + 7X, in 
Fig. 1.3. As mentioned, this operating mode is not a problem for the NST. However, the 
large Tesla coils use conventional transformers with per unit impedances in the range of 0.05 
to 0.1. A transformer with a per unit impedance of 0.1 will experience a current of ten times 
rated while the output is shorted. Most transformers do not survive very long under such 
conditions. Golka was not alone in burning out some of his transformers. The solution is to 
include additional reactance in the input circuit. 


The equivalent lumped circuit model of the Tesla coil while the gap is shorted is shown 
in Fig. 1.4. Ry and Rz are the effective resistances of the air cored transformer primary and 
secondary, respectively. The mutual inductance between the primary and secondary is shown 
by the symbol M. The coefficient of coupling is well under unity for an air cored transformer, 
so the ideal transformer model used for an iron cored transformer that electrical engineering 
students study in the first course on energy conversion does not apply here. 


At the time the gap arcs over, all the energy is stored in C,. As time increases, energy 
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Figure 1.3: Tesla Circuit With Gap Shorted. 
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Figure 1.4: Lumped Circuit Model Of A Tesla Coil, arc on. 


is shared among Cy, L1, Co, D2, and M. The total energy in the circuit decreases with 
time because of losses in the resistances R; and Ro. There are four energy storage devices 
so a fourth order differential equation must be solved. The initial conditions are some initial 
voltage vj, and 711 = 72 = vo = 0. If the arc starts again before all the energy from the previous 
arc has been dissipated, then the initial conditions must be changed appropriately. 


The Corums present the necessary solution technique in their manual [3] and also the 
computer code. The voltages and currents are not single frequency sinusoids. Rather there 
is a frequency spectrum with one hump for M small and two humps for M large. This 
is fascinating material for lovers of circuit theory, but is of somewhat limited usefulness in 
suggesting design changes for better performance. 


It appears to this author that the time domain solution is more useful than the frequency 
domain. We simply examine v1, v2, 41, and iz as time increases, either graphically or in some 
sort of tabular printout. We then change one or more of the energy storage device values 
and do it again. It is also helpful to calculate the energy stored in each device. If the total 
energy stored in the circuit is decreasing monotonically with time, at the rate power is being 
absorbed by R; and Rz, then one can be reasonably confident that the computer code is 
working correctly. 
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The time domain solution resembles a drunken walk in that it is difficult to predict what 
a given value will do next. Energy is moving among storage devices like cannon balls rolling 
around on the deck of an old sailing ship. Patterns can be changed readily by changing 
component values. We need a strategy for evaluating each solution for movement toward or 
away from some optimum. This strategy is developed by recognizing the following facts. After 
a small number of half cycles of 7;, the arc will dissipate and the spark gap will again become 
an open circuit. At this point we want as much energy as possible stored in the secondary, 
either as i3L2/2 or v3C2/2. Any energy stored in C; when the gap opens is not available to 
produce the desired high voltages on C). 


With proper design (proper values of C1, Ly, C2, L2, and M) it is possible to have all the 
energy in C transferred to the secondary at some time t,. That is, at t; there is no voltage 
across C; and no current through L,. If the gap can be opened at t;, then there is no way 
for energy to get back into the primary. No current can flow, so no energy can be stored 
in £,, and without current the capacitor cannot be charged. The secondary then becomes a 
separate RLC circuit with nonzero initial conditions for both Cy and L2, as shown in Fig. 1.5. 
This circuit will then oscillate or “ring” at a resonant frequency determined by Cy and Lg. 
With the gap open, the Tesla coil secondary is simply an RLC circuit, described in any text 
on circuit theory. The output voltage is a damped sinusoid. 


Lo Cy v2 


Figure 1.5: Lumped Circuit Model Of A Tesla Coil, arc off. 


Finding a peak value for v2 given some initial value for v, thus requires a two step solution 
process. We first solve a fourth order differential equation to find 72 and v2 as a function of 
time. At some time t, the circuit changes to the one shown in Fig. 1.5, which is described 
by a second order differential equation. The initial conditions are the values of 72 and v9 
determined from the previous solution at time t;. The resulting solution then gives the 
desired peak values for voltage and current. The process is tedious, but can readily be done 
on a computer. It yields some good insights as to the effects of parameter variation. It helps 
establish a benchmark for optimum performance and also helps identify parameter values that 
are at least of the correct order of magnitude. However, there are several limitations to the 
process which must be kept in mind. 


First, as we have mentioned, the arc is very difficult to characterize accurately in this 
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model. The equivalent R, will change, perhaps by an order of magnitude, with factors like 71, 
ambient humidity, and the condition, geometry, and temperature of the electrode materials. 
This introduces a very significant error into the results. 


Second, the arc is not readily turned off at a precise instant of time. The space between 
electrodes must be cleared of the hot conducting plasma (the current carrying ions and elec- 
trons) before the spark gap can return to its open circuit mode. Otherwise, when energy 
starts to bounce back from the secondary, a voltage will appear across the spark gap, and 
current will start to flow again, after the optimum time ¢t, has passed. With fixed electrodes, 
the plasma is dissipated by thermal and chemical processes that require tens of microseconds 
to function. When we consider that the optimum t; may be 2 ps, a problem is obvious. This 
dissipation time can be decreased significantly by putting a fan on the electrodes to blow the 
plasma away. This also has the benefit of cooling the electrodes. For more powerful systems, 
however, the most common method is a rotating spark gap. A circular disc with several elec- 
trodes mounted on it is driven by a motor. An arc is established when a moving electrode 
passes by a stationary electrode, but the arc is immediately stretched out by the movement 
of the disc. During the time around a current zero, the resistance of the arc can increase to 
where the arc cannot be reestablished by the following increase in voltage. 


The rotary spark gap still has limitations on the minimum arc time. Suppose we consider 
a disc with a radius of 0.2 m and a rotational speed of 400 rad/sec (slightly above 3600 rpm). 
The edge of the disc is moving at a linear velocity of rw = 80 m/s. Suppose also that an arc 
cannot be sustained with arc lengths above 2 cm. It requires 0.02/80 = 25 ys for the disc 
to turn this distance. This time can be shortened by making the disc larger or by turning 
it at a higher rate of speed, but in both cases we worry about the stress limits of the disc. 
Nobody wants fragments of a failed disc flying around the room. The practical lower limit 
of arc length seems to be about 10 ws. With larger coils this may be reasonably close to the 
optimum value. 


The third reason for concern about the above calculations is that the Tesla coil secondary 
has features that cannot be precisely modeled by a lumped circuit. One such feature is ringing 
at ‘harmonic’ frequencies. Neither the distributed or lumped models do a particularly good 
job of predicting these frequencies. Data will be presented later for a medium sized secondary 
(operated as an extra coil, explained in the next section), with a high Q resonance at about 160 
kHz. When applied power is switched off, the coil usually rings down at 160 kHz. Sometimes, 
however, it will ring down at 3.5(160) = 560 kHz. A third harmonic appears in many electrical 
circuits and has plausible explanations. A 3.5 ‘harmonic’ is another story entirely. 


These three reasons explain why we never see a paper giving a complete Tesla coil design 
with experimental data verifying the theoretical design. We get started with theory, but at 
some point have to move to an experimental optimization. The saying is, “Tune for most 
smoke”, which Harry Goldman attributes to Bill Wysock and Gary Legel. It is a tribute 
to the experimentalists that we have coils in existence with names like “Nemesis” that can 
produce sparks fifteen feet long [11]. 
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1.3. Magnifier 


As mentioned above, the classical Tesla coil uses two stages of voltage increase. Some coilers 
get a third stage of voltage increase by adding a magnifier coil, also called an extra coil, to 
their classical Tesla coil. This is illustrated in Fig. 1.6. 


Ch 
Va Ub : 3 Lo C2 


iron core air core magnifier 


Figure 1.6: The Classical Tesla Coil With Extra Coil 


The extra coil and the air core transformer are not magnetically coupled. The output 
(top) of the classical coil is electrically connected to the input (bottom) of the extra coil with 
a section of copper water pipe of large enough diameter that corona is not a major problem. 
A separation of 2 or 3 meters is typical. 


Voltage increase on the extra coil is by transmission line action, rather than the transformer 
action of the iron core transformer. Voltage increase on the air core transformer is partly by 
transformer action and partly by transmission line action. When optimized for extra coil 
operation, the air core transformer looks more like a transformer (greater coupling, shorter 
secondary) than when optimized for classical Tesla coil operation. 


The lumped circuit enthusiast would say that voltage rise is by RLC resonance. Both 
camps agree that voltage rise in the secondary and especially in the extra coil are not by 
transformer action. 


Although not shown in Fig. 1.6 the extra coil depends on ground for the return path of 
current flow. The capacitance from each turn of the extra coil and from the top terminal to 
ground is necessary for operation. Impedance matching from the Tesla coil secondary to the 
extra coil is necessary for proper operation. If the extra coil were fabricated with the same 
size coil form and wire size as the secondary, the secondary and extra coil tend to operate as 
a long secondary, probably with inferior performance to that of the secondary alone. There 
are guidelines for making the coil diameters and wire sizes different for the two coils, but 
optimization seems to require a significant amount of trial and error. 


In my quest for a better description of Tesla coil operation, I decided that the extra coil 
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was the appropriate place to start. It looks like a vertical antenna above a ground plane, so 
there is some prior art to draw from. While the classical Tesla coil makes an excellent driver 
to produce long sparks, it is not very good for instrumentation and measurement purposes. 
There are just too many variables. The spark gap may be the best high voltage switch available 
today, but inability to start and stop on command, plus heating effects, make it difficult to 
use when collecting data. 


I therefore decided to build a solid state driver. Vacuum tube drivers have been used for 
many years and several researchers have developed drivers using power MOSFETs, so this 
was not entirely new territory. It turned out to be a long term project. At the beginning, I 
had little idea about the input impedance of a coil above a ground plane, or how much power 
would be required to get significant sparks (say, half a meter in length or more). There have 
been many iterations, but I finally produced a design that would make sparks. Two major 
disadvantages are that it requires a digital oscilloscope with deep memory for tuning purposes, 
and one can make longer sparks using a standard spark gap. These disadvantages make it 
unlikely to sweep the Tesla coil community. There might be situations, however, where this 
approach would be useful. One is a museum installation, for example, where sparks of 0.5 to 
1 meter are acceptable, and long life and low maintenance are critical factors. 


The remainder of this document is a collection of my notes on this project, including some 
deadends. There are discussions on 

1. Capacitance 

2. Inductance and Transformers 

3. Gate Driver and Inverter 

4. Lumped Model 

5. Experimental Results 

Capacitance appears in many different places in the Tesla coil system, in the power supply, 


the controller, the driver, the coil body itself, and the top toroid or sphere. It therefore gets 
a lengthy treatment. Other items get a somewhat lesser treatment. 
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IDEAL CAPACITORS 


Capacitors are used in almost every activity of electrical engineering, yet information on 
capacitor characteristics is scattered through a variety of textbooks, databooks, and manu- 
facturers literature. The following is an attempt to organize some of this information. 


A typical capacitor is a two-terminal device consisting of two conductors separated by a 
dielectric. When a voltage difference V, is applied to the conductors, a charge of +Q will 
appear on one conductor and an equal and opposite charge —Q on the other conductor. The 
capacitance C’ is defined as the ratio of the charge on one conductor to the potential difference. 


Q 
=— 2.1 
= (2.1) 
where C is in farads, @ is in coulombs, and V, is in volts. Actually, one farad is a rather 
large capacitance, so capacitance values are usually expressed in terms of uF (10~°F) or pF 


(10> Fy, 
The total energy stored in a capacitor is 
1 1 1Q? 


1 
Wr=- E? — peas a 209: 
= sf dv = 5CV2 = 5QV. (2.2) 


where Wg is in joules, F is the electric field in V/m, and e€ is the permittivity. The integral 
expression shows that the energy stored in a capacitor with a fixed voltage difference across 
it increases as the permittivity of the material increases. 


The permittivity is usually expressed as the product of a relative permittivity ¢, and the 
permittivity of free space € 9. 


C= Gre5 (2.3) 


where 


€o = 8.854 x 107? F/m (2.4) 


The relative permittivity is unity for a vacuum and typically in the range of 2 to 6 for 
most dielectrics, as we shall discuss in more detail later. 


Capacitors are frequently used in series in a circuit, as shown in Fig. 2.1. There will be no 
actual charge transfer through the dielectric material. However, the electric fields will cause 
a movement of charge within the series string. The battery supplies a positive charge to the 
left plate of C;. This positive charge attracts an equivalent negative charge on the right plate 
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of Ci. The movement of this negative charge leaves behind a positive charge of the same 
amount, which the electric field will force onto the left plate of Co. The process continues 
until each capacitor has the same charge Q = Q, on its left plate. That is, 


Qs=Q1=Q2 =Q3 


(2.5) 
The total voltage across the series combination is 


V=H=V+K4+ V3 (2.6) 
and since 


ee 


a (2.7) 


OO Oe 


Gy Cy Ca C3 


(2.8) 
which can be solved for the series capacitance Cs. 


(2.9) 


Figure 2.1: Capacitors in Series 


A circuit of parallel capacitors is shown in Fig. 2.2. The voltage on each capacitor is the 
same and the amount of stored charge on each capacitor will be proportional to the individual 
capacitance values. It is not hard to show that the total parallel capacitance C, is given by 


Cp =C1+C.+C3 (2.10) 
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Ci C2 C3 


Figure 2.2: Capacitors in Parallel 
2.1 Capacitance of Common Geometries 


The ratio of Q to V, depends on the geometrical arrangement of the conductors and on the 
electrical characteristics of the dielectric. The capacitance of a parallel plate capacitor, as 
illustrated in Fig. 2.3, is 


cA 
C=— 2.11 
< (2.11) 


where A is the area in m? and d is the separation between plates. This formula is accurate 
only when fringing can be neglected, that is, when d is small in comparison with A. 


Area A 


Figure 2.3: A Parallel Plate Capacitor 


Later on, we will be interested in the capacitance of geometries where there are two 
different dielectrics. The simplest case is shown in Fig. 2.4. There is a layer of dielectric with 
relative permittivity €, > 1 of thickness x, and a second layer of air, with thickness y. The 
boundary between the two dielectrics can be considered a floating electrode. In fact, we can 
place a conducting plate on the boundary without changing the results at all. We basically 
have two capacitors in series. When we solve for the series capacitance, we get 


6A Ep 
¥y (€, + x/y) 


(2.12) 
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—{]:}— 


a oo 


Figure 2.4: Capacitor with Different Dielectrics 


Another important geometry is that of a coaxial cable of inner radius a, outer radius 8, 
and length @, which has capacitance 


Qrel 
C= in(b/a} (2.13) 


The geometry for the coaxial cable is shown in Fig. 2.5. 


(@) eee] + 


ee 9 - 1 —_J 


Figure 2.5: Coaxial Cable 


The common 50 2 coaxial cable 213/U (RG-8A/U) has a nominal capacitance of 29.5 
pF/ft (96.8 pF/m). The small 75 Q video cable 59B/U has a nominal capacitance of 21.0 
pF/ft (68.9 pF/m). Most other 50 and 75 Q cables will have capacitance values very close to 
these. Physically larger cables capable of carrying more power will have 6 and a increased in 
the same proportion so the ratio b/a and the capacitance will remain the same as that of a 
smaller cable. 


Another geometry of great practical interest is the twin conductor transmission line, shown 
in Fig. 2.6a. This is composed of two conductors of radius r, with separation 2h between 
conductor centers. The conductor-to-conductor capacitance C;,, is given by 


mel mel 
OBS Inf(h + Vh2 —12)/r] cosh 1(h/r) iy 


If the two conductors have a small radius and are located far apart, the expression for 
capacitance becomes 
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mel 
Cee = In(2h/r) (2:15) 


The error in the approximate expression is only 5.26% when h = 2r and 1.16% when 
h = 3r so the latter equation really has a wide range of usefulness. 


The conductor-to-plane capacitance C,, between a cylindrical conductor of radius r and 
a conducting plane a distance h from the cylinder, as shown in Fig. 2.6b, is twice the value 
given by the previous two equations. 


27el 27el 
(ey Ce 2.16 
> cosh-t(h/r) — In(2h/r) a) 


This equation can be used to find the capacitance between two unequal conductors. We 
find the capacitance of each conductor to an imaginary ground plane, and then combine the 
two values for C,, using the formula for capacitors in series. 


ai ——— 
gy 
2h Cee 
(a) (b) 
Figure 2.6: Twin Conductor Transmission Line 


Another geometry of interest is that of a spherical capacitor of two concentric spheres with 
radii a and b (b > a) as shown in Fig. 2.7. It is not practical to actually build capacitors this 
way, but the symmetry allows an exact formula for capacitance to be calculated easily. This 
is done in most introductory courses of electromagnetic theory. The capacitance is given by 
[3, Page 165] 


Are 


aT (2.17) 


If the outer sphere is made larger, the capacitance decreases, but does not go to zero. In 
the limit as b > oo, the isolated or isotropic capacitance of a sphere of radius a becomes 


Coo = 4r€a (2.18) 
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1OL 


Figure 2.7: Spherical Capacitor 


Coo gives us a lower bound for the capacitance of a spherical top loading element of a 
Tesla coil with respect to ground. One way of arriving at a reasonable estimate of the actual 
capacitance is to start with the isotropic capacitance C, and add a correction term, as we 
shall see in the next section. 


We will also be interested is isotropic capacitances of shapes other than spheres. Some will 
be difficult to impossible to calculate analytically, so we will use the isotropic capacitance of 
a sphere as a starting point for making an estimate. Rectangular boxes and short cylinders, 
for example, will have similar isotropic capacitances to a sphere. We just need to find an 
equivalent radius (or diameter) for these nonspherical shapes. Several different equivalents 
could be used, such as a geometrical mean equivalent, an arithmetic mean equivalent, or the 
radius of a sphere with the same surface area as the other shape. It turns out that the simplest 
method, the arithmetic mean, does quite well [1]. Consider a rectangular box with orthogonal 
edge dimensions a, b, and c. Define an equivalent sphere diameter ¢. where 


0. = avere (2.19) 


For other shapes we use the dimensions of the box in which the other shape can be placed. 
Making the change from radius to diameter, the isotropic capacitance is now 


Coo = 2r1ebe (2.20) 


This approach will give acceptable results in many cases. But, of course, there is no way 
of knowing the amount of error, or when some other approach would yield better results. It 
will get us in the right ballpark, however, and sometimes allow us to determine lower bounds 
of acceptable values obtained from other techniques. 


Suppose, for example, that we wanted the capacitance between two spheres separated 
by several diameters. We suspect that the parallel plate capacitor formula will not be very 
accurate and are unable to locate a better formula. What can we do? The lower bound for the 
capacitance between two spheres is just half of Co, the isotropic capacitance for one sphere, 
as can be argued from Fig. 2.8. 


The dark point at the center of the figure is obviously ‘the’ point at infinity that math- 
ematicians love to talk about. If the capacitance of each sphere with respect to this ‘point’ 
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Figure 2.8: Capacitance Between Two Spheres 


is C.., then the capacitance between spheres has a lower bound of C,,/2, from the formula 
for two capacitors in series. Bringing the spheres closer together will increase the capacitance 
but they cannot be separated far enough to reduce the capacitance below C,, /2. 


2.2 Toroid Capacitance 


An important emphasis of this book is the analysis of the Tesla coil. Among other things we 
will be interested in the capacitance of the top element (usually a ‘fat’ toroid) with respect 
to ground. We will also need the capacitance between adjacent turns (which look like ‘thin’ 
toroids) and finally the capacitance of a turn with respect to ground. As usual, we will use 
the published results as much as possible and leave the derivations to others. 


The dimensions of a toroid are shown in Fig. 2.9. 


top 


view 


k—— p ———| side 
O O fa 


Figure 2.9: Toroid Dimensions 


There are other coordinate systems besides rectangular, cylindrical, and spherical in which 
variables can be separated and Laplace’s equation solved. One of these is toroidal coordinates. 
Moon and Spencer use this coordinate system to solve for the capacitance of an isolated toroid 
(2, Page 375] as 
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Q_1/2(cosh mo)  Qn—1/2(cosh no) 
P_4/2(cosh No) Py-1/2(cosh no) 
where P and Q are Legendre functions of first and second order and a and 1 will be discussed 


later. The subscript ‘MS’ refers to Moon and Spencer, to distinguish the capacitance obtained 
from the value to be obtained from some empirical formulas later. 


Cys = 8ae (2:21) 


n=. 


Note that Eq. 2.21 is the corrected version. Godfrey Loudner [1] found that Moon and 
Spencer had an extra 7 in their expression (second equation from the top on p. 375), and also 
in Eq. 13.29 and Eq. 13.29a. There is also a typo on p. 373, second equation from the top, 
where V/2 should be replaced by V/z. Even the great ones can make a mistake! 


The Legendre functions can be written in many different forms, as integrals or infinite 
series, converging for arguments greater than unity or less than unity, and so on. The non- 
integer subscript (n—1/2) adds another layer of complexity. Many math books do not mention 
the non-integer case, so one must be diligent in finding the correct expressions. The order of 
difficulty is much greater than for the sphere. Our mothers warned us that there would be 
days like this, but let us proceed. 


Moon and Spencer give us an expression for Q,,_1/2(cosh 7) [2, Page 373] 


1 us cos n6 dé 


om h — mmm": 
Qn—1/2(cosh no) V2 Jo ./(cosh jo — cos 8) 


For some reason, they do not give a similar expression for P,_/2(cosh 9). However, they 
do give plots for both functions, which is convenient for checking computational results. 


(2.22) 


Smythe [7, Page 159] gives an expression for P in the form 


nPi (x) = (n' +1)(n' + 2)+--(n’ +m) [e+ Vx? — 1 cos 6)" cos mé dé (2.23) 


This expression is valid for x > 0 and for any n’, including 1/2, 3/2, etc. We are only 
interested in the case m = 0. If we interpret the product of terms ahead of the integral sign 
as a factorial with zero entries (that is with a value of unity), the expression becomes 


1 [7 / 
Py (a) = - | [cv + Vx? —1cos6|” dé (2.24) 
T JO 
These expressions for P and Q can be numerically integrated using any scientific program- 
ming language (QuickBasic, etc.). 


We now return to our discussion of a and 9. These are coordinate values in toroidal 
coordinates, hence need to be translated into a more familiar coordinate system. The major 
toroid radius w and the minor radius r are given by 
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w = acothn (2.25) 
a 

= Jos 2.26 

sinh ( ) 


Eliminating a between the two equations yields 


w =rsinhncothn = rcoshn (2:27) 


so ‘5 
— =coshn=2 (2.28) 
r 


We then solve for 7 as 


ar) 
eee — | (2.29) 
r 
and for a as 
a=rsinhyn = s(e" —e€") (2.30) 


This ‘exact’ formulation is of most interest to EM theorists and computer programmers. It 
will seem like overkill to most Tesla coil enthusiasts who just need to get in the right ballpark 
with a capacitance estimate. For this reason empirical formulas have been developed which 
yield an approximate value, adequate for most purposes, but obtained with much less effort. 
Empirical formulas for the capacitance of a toroid are given by [6] 


Cs = maa (d/D < 0.25) (2.31) 
Cs = 0.37D + 0.23d (d/D > 0.25) (2.32) 


where D is the toroid major diameter, outside to outside, in cm, dis the toroid minor diameter 
in cm, and the capacitance is given in pF. Table 2.1 gives some isotropic capacitance values, 
both from the Moon and Spencer numerical integration and the empirical formulas. The 
deviation or error of C's with respect to Cyyg is given in the last column in percent. 


We see that the empirical formulation agrees with Moon and Spencer to within 1% for the 
case of fat toroids, but gets progressively worse as the toroid gets thinner. The error is within 
5% for d down to about 0.5 cm (4 gauge wire). Toroids this thin do not have the mechanical 
strength necessary to serve as top loading elements of a Tesla coil, so we can conclude that 
the empirical formulas are quite adequate for most purposes. 
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Table 2.1: Isotropic Capacitance of Toroids 
w r D d Cys Cg error 
meters cm pF pF % 

3. «15 90 30 40.46 40.20 —0.63 
2 Ll 60 20 26.97 26.80 -0.63 
1.05 30 10 13.49 13.40 —0.63 
2 .08 56 16 24.55 24.40 —0.62 
2 .06 52. 12 22.02 21.93 —0.42 
2 .04 48 19.28 19.52 +1.23 
2. 402 44 16.00 16.43 +2.70 
2 OL 42 +3.42 
2 = .005 41 12.00 12.48 +4.05 
2 .0025 40.5 10.62 11.14 +4.97 
2 .001 40.2 9.09 9.76 +7.36 


Doar bw & 
e 
ed 
NI 
i) 
a 
i 
e 
© 


2.3 Solenoid Capacitance (Medhurst) 


The isotropic capacitance of a sphere was given above as a simple formula. We looked at 
the theoretical formulas for capacitance of a toroid, but basically gave up and went to a 
simpler empirical version. After that learning experience, we will not even try to write exact 
equations for the isotropic capacitance of a cylinder. We will immediately write the empirical 
equations as developed many years ago by a man named Medhurst. These will be expressed 
in several different versions, to meet different needs. The simplest expression for the isotropic 
capacitance of a cylindrical coil of wire, with diameter D and coil length @, is 


Cu =HD pF (2.33) 


where D is in cm, and H is a multiplying factor that equals 0.51 for 0/D = 2, 0.81 for 0/D 
= 5, and varies linearly between 0.51 and 0.81 for ¢/D between 2 and 5. Most coilers prefer 
values for £/D between 3.5 and 4.5, so this linear range is adequate for most purposes. 


An expression for H that works for ¢/D between 2 and 8 is 


£ 
H= 0.100976-—5 + 0.30963 (2.34) 


Another expression for H that works for 0/D between 1 and 8 is 
H=0 0005(—)4 —0 o097(£)3 +0 0648(—)? —0 0757() + 0.4723 (2.35) 
=), D : D D D : 
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2.4 Tesla Coil Capacitance 


We now have expressions for the isotropic capacitance Cg of a toroid and the isotropic ca- 
pacitance Cyy of a coil. The next step would be to set the toroid on top the coil and add 
the two capacitances to get an effective capacitance C;,. for the Tesla coil. Unfortunately, this 
only works when the toroid is far away from the solenoid. As the toroid is brought near the 
coil form, shielding occurs such that the effective capacitance is less than the sum of the two 
isotropic capacitances. The Tesla coil capacitance might be written as 


Cre = Cu + KC (2.36) 


where K <1. A value of K = 0.75 should result in a number for C;, within 20% of the correct 
value for most Tesla coils. The resonant frequency is related to the square root of Ci. so a 
20% error in capacitance results in only a 10% error in resonant frequency. 


Most readers probably feel disappointed here. We have gone to considerable effort and 
still come up short of an accurate formula for C;,.. Our effort is not entirely wasted because 
we can do ‘what if’ analyses relatively quickly. Questions about the effect of changing coil 
diameter, coil length, or toroid diameter can be answered with adequate accuracy. 


Someone might suggest using a modern digital capacitance meter to measure C;,. This 
method would probably have greater error than the above formula, because the leads of the 
capacitance meter have a similar capacitance value as C;,. Also the presence of the meter and 
a human body will change the capacitance. 


It is possible to calculate C;, numerically using Gauss’s Law. If one is careful about 
measuring and entering all the dimensions and the locations of grounded surfaces, one should 
get a value for C;. well within 5% of the correct value. There are programs available in the 
Tesla coil community that do this. 
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LOSSY CAPACITORS 


3.1 Dielectric Loss 


Capacitors are used for a wide variety of purposes and are made of many different materials 
in many different styles. For purposes of discussion we will consider three broad types, that 
is, capacitors made for ac, dc, and pulse applications. The ac case is the most general since 
ac capacitors will work (or at least survive) in de and pulse applications, where the reverse 
may not be true. 


It is important to consider the losses in ac capacitors. All dielectrics (except vacuum) have 
two types of losses. One is a conduction loss, representing the flow of actual charge through 
the dielectric. The other is a dielectric loss due to movement or rotation of the atoms or 
molecules in an alternating electric field. Dielectric losses in water are the reason for food and 
drink getting hot in a microwave oven. 


One way of describing dielectric losses is to consider the permittivity as a complex number, 
defined as 
e=e — je" =|ele” (31) 


where 
e’ = ac capacitivity 
é€” = dielectric loss factor 
6 = dielectric loss angle 


Capacitance is a complex number C™ in this definition, becoming the expected real number 
C’ as the losses go to zero. That is, we define 


C*=C-— jC" (32) 


One reason for defining a complex capacitance is that we can use the complex value in any 
equation derived for a real capacitance in a sinusoidal application, and get the correct phase 
shifts and power losses by applying the usual rules of circuit theory. This means that most of 
our analyses are already done, and we do not need to start over just because we now have a 
lossy capacitor. 


Equation 3.1 expresses the complex permittivity in two ways, as real and imaginary or 
as magnitude and phase. The magnitude and phase notation is rarely used. Instead, people 
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usually express the complex permittivity by ¢«’ and tan 6, where 


" 


€ 
tand = a (3.3) 
where tan is called either the loss tangent or the dissipation factor DF. 
The real part of the permittivity is defined as 
e. Seep (3.4) 


where e€, is the dielectric constant and €, is the permittivity of free space. 


Dielectric properties of several different materials are given in Table 3.1 [15, 5]. Some of 
these materials are used for capacitors, while others may be present in oscillators or other 
devices where dielectric losses may affect circuit performance. The dielectric constant and 
the dissipation factor are given at two frequencies, 60 Hz and 1 MHz. The righthand column 
of Table 3.1 gives the approximate breakdown voltage of the material in V/mil, where 1 mil 
= 0.001 inch. This would be for thin layers where voids and impurities in the dielectrics are 
not a factor. Breakdown usually destroys a capacitor, so capacitors must be designed with a 
substantial safety factor. 


It can be seen that most materials have dielectric constants between one and ten. One 
exception is barium titanate with a dielectric constant greater than 1000. It also has relatively 
high losses which keep it from being more widely used than it is. 


We see that polyethylene, polypropylene, and polystyrene all have small dissipation fac- 
tors. They also have other desirable properties and are widely used for capacitors. For high 
power, high voltage, and high frequency applications, such as an antenna capacitor in an AM 
broadcast station, the ruby mica seems to be the best. 


Each of the materials in Table 3.1 has its own advantages and disadvantages when used in 
a capacitor. The ideal dielectric would have a high dielectric constant, like barium titanate, a 
low dissipation factor, like polystyrene, a high breakdown voltage, like mylar, a low cost, like 
aluminum oxide, and be easily fabricated into capacitors. It would also be perfectly stable, 
so the capacitance would not vary with temperature or voltage. No such dielectric has been 
discovered so we must apply engineering judgment in each situation, and select the capacitor 
type that will meet all the requirements and at least cost. 


Capacitors used for ac must be unpolarized so they can handle full voltage reversals. They 
also need to have a lower dissipation factor than capacitors used as dc filter capacitors, for 
example. One important application of ac capacitors is in tuning electronic equipment. These 
capacitors must have high stability with time and temperature, so the tuned frequency does 
not drift beyond some specified amount. 


Another category of ac capacitor is the motor run or power factor correcting capacitor. 
These are used on motors and other devices operating at 60 Hz and at voltages up to 480 
V or more. They are usually much larger than capacitors used for tuning electronic circuits, 
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and are not sold by electronics supply houses. One has to ask for motor run capacitors at an 
electrical supply house like Graingers. These also work nicely as dc filter capacitors if voltages 
higher than allowed by conventional dc filter capacitors are required. 


The term power factor PF may also be defined for ac capacitors. It is given by the 
expression 


PF = cos (3.5) 
where @ is the angle between the current flowing through the capacitor and the voltage across 
it. 


The capacitive reactance for the sinusoidal case can be defined as 


Xo = — : 
C= TE (3.6) 


where w = 27f rad/sec, and f is in Hz. 


In a lossless capacitor, ¢” = 0, and the current leads the voltage by exactly 90°. If €” is 
greater than zero, then the current has a component in phase with the voltage. 


cos 6 = ia? (3.7) 


For a good dielectric, €’ >> €”, so 


" 


cos6 x — =tanéd (3.8) 


é 
Therefore, the term power factor is often used interchangeably with the terms loss tangent 
or dissipation factor, even though they are only approximately equal to each other. 


We can define the apparent power flow into a parallel plate capacitor as 


2 


A A 
62 VIS ey BO ay =e V9 erelj + DF) (3.9) 


By analogy, the apparent power flow into any arbitrary capacitor is 
S=P+jQ=V’wC(j + DF) (3.10) 
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Table 3.1: Dielectric Constant 
selected materials. 


Material 


Air 

Aluminum oxide 
Barium titanate 
Carbon tetrachloride 
Castor oil 

Glass, soda-borosilicate 
Heavy Soderon 

Lucite 

Mica, glass bonded 


Mica, glass, titanium dioxide 


Mica, ruby 
Mylar 

Nylon 

Paraffin 
Plexiglas 
Polycarbonate 
Polyethylene 
Polypropylene 
Polystyrene 
Polysulfone 


3-4 


€r, Dissipation Factor DF and Breakdown Strength V, of 


€r 
60 Hz 


1.000585 
1250 
2.17 

3.7 


Polytetrafluoroethylene(teflon) yal 


Polyvinyl chloride (PVC) 
Quartz 

Tantalum oxide 
Transformer oil 

Vaseline 
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10° Hz 


1.000585 


8.80 
1143 
2.17 
3.7 
4.84 
3.39 
3.3 
7.39 
9.0 
5.4 
2.5 
3.33 
2.25 
2.76 
2.7 
2.26 
2.25 
2.56 
3.1 


0.005 


0.014 


0.06 
<0.0002 
<0.0005 

<0.00005 

<0.0005 
0.0115 
0.0009 


0.0004 


DF Vb 
10° Hz = =V/mil 


- 79 
0.00033 300 
0.0105 50 

<0.00004 - 
- 300 
0.0036 - 
0.0283 - 
- 500 
0.0013 1600 
0.0026 - 
0.0003 - 
- 5000 
0.026 - 

- 250 
0.014 - 

- 7000 


<0.0002 4500 
<0.0005 9600 
0.00007 500 
- 8000 
<0.0002 1500 
0.016 


0.0001 500 
- 100 
- 250 

<0.0001 - 
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The power dissipated in the capacitor is 


P=V*wC" = V*wC(DF) (3.11) 


Example 


Find the real and reactive power into a ruby mica capacitor with area A = 0.03 m?, and a dielectric 
thickness d = 0.001 m, if the voltage is 2000 V (rms) at a frequency f = 1 MHz. 


f= v2 cptoli + DF) 


10°) (0.03) 


2 
= (2000)22% TOOL (O4) (8-854 x 10-7?)(7 + 0.0003) 


= 736040 + 10.8 


The capacitor is absorbing 36040 capacitive VARs (Volt Amperes Reactive) and 10.8 Watts. The 
real power of 10.8 W appears as heat and must be removed by appropriate heat sinks. 


The real power dissipation in a capacitor varies directly with frequency if the dissipation 
factor remains constant, and also with the square of the voltage. At low frequencies, the 
voltage limit is determined by the dielectric strength. At high frequencies, however, the 
voltage limit may be determined by the ability of the capacitor to dissipate heat. If the ruby 
mica capacitor in the previous example could safely dissipate only 10 W, but was to be used at 
5 MHz, the operating voltage must be reduced from 2000 V to keep the losses in an acceptable 
range. 


The dissipation factor varies significantly with frequency for some materials in Table 3.1. 
The actual variation must be determined experimentally. If interpolation or extrapolation 
seems necessary to find the loss at some frequency not given in the Table, the best assumption 
would be that DF varies linearly with log f. That is, if DF = 0.02 at f = 10? and 0.01 at f 
= 10°, a reasonable assumption at f = 104 would be that DF = 0.015. 


The basic circuit model for a capacitor is shown in Fig. 3.1. Any conductor, whether 
straight or wound in a coil, has inductance, and the capacitor inductance is represented by a 
series inductance L,. The effects of conductor resistance and dielectric losses are represented 
by a series resistance R,. Leakage current through the capacitor at de flows through a parallel 
resistance R,. In some manufacturer’s databooks, these are called ESL, ESR, and EPR, where 

L, = ESL = equivalent series inductance 


R, = ESR = equivalent series resistance 


R, = EPR = equivalent parallel resistance 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 3—Lossy Capacitors 3-6 


Ls Rs C 
YY 
ESL ESR 
Rp 
EPR 


Figure 3.1: Equivalent Circuit for a Capacitor 


This circuit indicates that every capacitor has a self-resonant frequency, above which it 
becomes an inductor. This is certain to puzzle a student making measurements on a capacitor 
above this frequency if the student is not aware of this fact. R, is readily measured by applying 
this frequency to a capacitor, measuring the voltage and current, and calculating the ratio. 
The capacitive and inductive reactances cancel at the resonant frequency, leaving only R, to 
limit the current. The resistance R, will always be much larger than the capacitive reactance 
at the resonant frequency, so this resistance can be neglected for this computation. 


The self-resonant frequency of a high capacitance unit is lower than that for a low ca- 
pacitance unit. Hence, in some circuits we will see two capacitors in parallel, say a 10 uF in 
parallel with a 0.001 uF capacitor as shown in Fig. 3.2. At first glance, this seems totally 
unnecessary. However, the larger capacitor is used to filter low frequencies, say in the audio 
range, while the small capacitor filters the high frequencies which are above the self-resonant 
frequency of the large capacitor. 


10 pF 0.001 pF 


Figure 3.2: Capacitors in Parallel to Filter Two Different Frequencies 


An example of the variation of R, and self-resonant frequency f,e, with capacitance value 
and rated voltage is given in Table 3.2. These are metallized polypropylene capacitors designed 
for switch-mode power supplies by the company Electronic Concepts, Inc. This application 
requires a low R, and a high f,-, so capacitors designed for other applications will tend to 
have higher values of R, and lower values of fires. 


The operating frequency range will obviously be less than the self-resonant frequency. This 
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Table 3.2: Capacitor Resistance and Self-resonant Frequency for Electronic Concepts Type 
5MP Metallized Polypropylene capacitors. 


VDC C Rs ares 
volts pF Q kHz 


100 1 0.015 1065 
100 2 0.012 703 
100 5 0.010 385 
100. =10 0.009 248 
200 1 0.020 861 
200 2 0.015 609 
200 5 0.011 323 
200 10 0.009 200 
400 1 0.019 784 
400 2 0.015 511 
400 5 0.010 283 
400 10 0.006 200 


could easily be as low as a few kHz for large motor run capacitors. 


Capacitors will always be rated for a working voltage as well as a specific value of capaci- 
tance. This voltage will always be well under the breakdown voltage of the dielectric. It may 
be specified either as a dc voltage or an ac voltage, depending on the application. Motor run 
capacitors are always operated on ac, so the voltage is specified as, say, 370 or 480 VAC. They 
can also be operated on dc, with a de rating of at least \/2 times the ac rating. Electrolytic ca- 
pacitors that can only be operated on dc will have their working voltage expressed as WVDC. 
This is the maximum dc voltage, plus the peak of the ac ripple voltage, that should be con- 
tinuously applied to a capacitor to prevent excessive deterioration and aging. Capacitors that 
are sometimes used on ac, and sometimes on dc, usually have working voltages expressed as 
WVDC, as was the case for the polypropylene capacitors in Table 3.2. 


We will now present some more detailed background information on several of the di- 
electrics shown in Table 3.1. 


Mica is a natural material that can be easily split into thin layers. It is very stable and 
does not deteriorate with age. The maximum capacitance is on the order of 0.03 wF. Mica 
capacitors tend to be quite expensive and the larger sizes are used only in critical applications 
like radio transmitters. 


Glass capacitors were first developed during World War II as a replacement for mica 
capacitors when supplies of mica were threatened. Glass capacitors exhibit excellent long- 
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term parametric stability, low losses, and can be used at high frequencies. They are used 
in aerospace applications where capacitance must not vary. The maximum size is limited to 
about 0.01 uF. 


Paper capacitors use kraft paper impregnated with non-ionized liquid electrolyte as the 
dielectric between aluminum foils. These capacitors can provide large capacitances and voltage 
ratings but tend to be physically large. The thickness of the paper layers and the particular 
electrolyte used can be varied to produce a wide range of electrical characteristics, which is 
the reason for not listing paper in the dielectrics of Table 3.1. 


Ceramic capacitors are made with one of a large number of ceramic materials, which 
include aluminum oxide, barium titanate, and porcelain. These are very widely used as bypass 
capacitors in electronic circuits. The older style is a single layer of dielectric separating two 
conducting plates and packaged in a small disc. The newer style is the monolithic, which 
appears in a rectangular package. It consists of alternating layers of ceramic material and 
printed electrodes which are sintered together to form the final package. The self resonant 
frequency is on the order of 15 MHz for the 0.01 uF size, for a total lead length of 0.5 inch, and 
on the order of 165 MHz for the 0.0001 uF = 100 pF size. Ceramic capacitors are classified 
into Class 1 (€, < 600) and Class 2 (€, > 600) by the Electronics Industries Association 
(EIA). Class 1 ceramic capacitors tend to be larger than their Class 2 counterparts, and have 
better stability of values with changes in temperature, voltage, or frequency. The best Class 1 
capacitor, with nearly constant characteristic, will be labeled ‘COG’ using EIA designators, but 
is often referred to as ‘NPO’ which stands for ‘negative-positive-zero’ (temperature coefficient). 


Plastic-film capacitors use extremely thin sheets of plastic film as the dielectric between 
capacitor plates, usually in a coil construction. They have low losses and good resistance to 
humidity. Four common materials are: 


— 


. Polycarbonate 


i) 


. Polypropylene 
3. Polystyrene 


4. Teflon 


Generally speaking, as one moves down in this list from polycarbonate to Teflon, the 
capacitors get larger, better, and more expensive. 


Historically, low-budget Tesla coilers have used either saltwater capacitors (as Tesla himself 
did) or homemade rolled polyethylene and foil capacitors. The saltwater capacitors are lossy 
and heavy. They are one of the first components to be replaced as a new coiler grows with 
his hobby. The homemade foil capacitors are limited by corona onset to about 7000 V. They 
are oil filled, so if the container leaks or is tipped over, one has a real mess. Serious coilers 
would look for commercial high voltage, high current, pulse rated capacitors made by specialty 
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companies like Maxwell. These worked fine, but were expensive if purchased new, and difficult 
to find on the used market. 


This all changed in 1999 when the Tesla coil community moved en masse to multi mini 
capacitors (MMC). These are small commercial capacitors that one buys by the sackful from 
Digi-Key and connects in series and parallel strings to get the required ratings. 


There are several different capacitor series, some of which are more suitable for Tesla coil 
work than others. One should look for capacitor types that are rated for “high voltage, high 
frequency, and high pulses”. A dissipation factor of 0.1% at 1 kHz is good, as is polypropylene 
for the dielectric. The ECWH(V) and ECQP(U) series of capacitors are good. The Panasonic 
ECQ-E capacitors utilize metallized polyester (Mylar) which has a higher dissipation factor 
(1% at 1 kHz) than the polypropylene capacitors. 


Capacitor ratings have received a great deal of attention by Tesla coilers. A given capacitor 
might be rated at 1600 VDC and 500 VAC. The AC rating is always for rms values, so the peak 
voltage would be about 700 V for a rating of 500 VAC. A capacitor in a Tesla coil primary 
experiences full voltage reversal, so it would seem wise to find the AC rating, multiply by 
V2, and divide that into the peak voltage available from the iron core transformer to find the 
number of series capacitors to use. 


Experience of coilers, however, has been that there are sufficient safety factors built into 
the capacitors that using the DC rating is quite acceptable. For example, a 15 kV transformer 
has a peak voltage of 15/2 = 21.2 kV. Divide 21.2 by 1.6 for 1600 VDC capacitors to get 
13.25. Use either 13 or 14 capacitors in series. 


The voltage rating is not determined by the dielectric breakdown voltage as much as it is 
by the Ionization or Corona Inception Level. This refers to the level at which ionization or 
partial discharges can begin to occur inside a bubble of entrapped air or within an air-filled 
void within the solid dielectric system. If one had perfect dielectrics and could always exclude 
any entrapped air, derating for this phenomenon would not be necessary. 


Corona inside a capacitor will chemically degrade the dielectric material over a period of 
time until the capacitor ultimately fails. Manufacturers of capacitors might select a voltage 
rating whereby their capacitors will last for one million hours before this occurs. In Tesla coil 
use, 100 hours of actual operation is a long time. This helps explain why coilers can exceed 
the manufacturer ratings without immediate problems. 


There is evidence that operating a capacitor above its ac voltage rating will shorten its 
life by a factor of the overvoltage ratio raised to the 15th power. Suppose we have a capacitor 
rated at one million hours at 500 VAC. For economic reasons we are thinking about operating 
it at 700 VAC. The life reduction factor is (700/500)'° = 155.57. Dividing one million by 
155.57 gives an expected life of 6428 hours or almost one year of continuous operation. This 
is more than adequate for most Tesla coil applications. 


One brand of capacitor is the WIMA. Terry Fritz comments about them: “The switching 
power supplies we build have WIMAs in pulsed duty similar to that seen in Tesla coil use. 
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The real indicator of how long they live is how warm they get. If they heat to about 5 
degrees C above the ambient, then they last forever. At about 8 degrees above ambient, we 
see occasional failures. At 10 degrees, they get to be a problem. In many situations, we go 
beyond the WIMA chart derating in pulsed applications with no problem at all, just being 
sure they don’t get hot.” 


WIMA experts say that partial discharges occur only above a certain voltage level and 
only on ac. Frequency does not seem to be a factor. They agree that the ac rating can 
be exceeded in Tesla coil applications, but without a serious investigation, were mentioning 
factors like 1.25 to 1.5. 


Some capacitors that are candidates for a MMC use actual metal foil for electrodes while 
others use metallized plastic. The metal foil devices are physically larger for a given rating 
and are much more robust in Tesla coil service. 


Dielectric material, electrode thickness, and lead size all limit the rate at which charge 
can be added to or removed from a capacitor. This limit is expressed as the dV/dt limit. The 
peak current Ipeqx can be determined from the dV/dt value by the equation 


Lae “C (3.12) 
where appropriate multipliers should be used to get everything in volts, seconds, and farads. 
By way of reference, a 1000 pF capacitor rated for 10,000 V/s will withstand 10 A surges. 
A capacitor with a dV/dt rating less than 1000 V per 5 ys probably does not have metal foil 
endplates and should not be used. 


The peak current that a Tesla coil primary capacitor must provide is determined by the ca- 
pacitor voltage and the surge impedance. If the secondary were removed, the surge impedance 
would be 


Vey ae (3.13) 


where C; is the capacitance in the Tesla coil primary and L, is the inductance of the primary. 
When the secondary is in place, the surge impedance will increase from the above value, so 
the peak current required from the capacitor will be less. Using the above expression for Z, 
should be a worst case calculation. If it tells you that there might be 100 A flowing after the 
gap shorts, and your individual capacitors are rated at 10 A each, then you would need ten 
parallel strings. Since this figure includes some factor of safety, you might be able to get by 
with somewhat fewer strings. One WIMA expert said that the maximum amperage rating is 
not as critical as the Corona Inception Level mentioned above. 


The polypropylene capacitors useful for Tesla coil work have a self-healing mode. A partial 
discharge will eventually cause a local failure. A burst of energy through this short will 
vaporize everything around the short, effectively removing the short. The capacitor continues 
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to operate, but with a slightly lower capacitance due to some of the electrodes being removed. 
It is a good idea to carefully measure the capacitance of the MMC periodically. A decrease of 
even 1% could indicate that the capacitors are being stressed, and that total failure is possible. 


Terry Fritz reports testing some WIMA MKS 4 (K4) capacitors rated at 1 uF and 400 
VDC on a large DC supply. These are metallized polypropylene dielectric capacitors that are 
encapsulated in epoxy filled 3 x 1 cm rectangles. Starting at about 900 V there are a few little 
snaps, indicating that the dielectric had punched through and the arc blasted the thin metal 
layer on the other side. Going to 1200 V gave a few more snaps. Around 1500 V there were 
many snaps and the capacitor shorted out. 


When a capacitor cleanly blows apart, that is a sign that they were destroyed by over 
voltage. When a capacitor puffs up and looks melted and burned, that is from too much 
current. It is always a good idea to think about possible failure modes, and put the capacitors 
in appropriate enclosures to protect bystanders. 


IT had a student build a small design project once, that illustrates this point. He used an 
electrolytic capacitor rated at about 25 VDC in some circuit for 120 VAC usage, and I did not 
notice it before testing. He was proudly demonstrating performance, using a variac to bring 
up the voltage. At about 90 V, the electrolytic capacitor exploded with a noise somewhere 
between a pistol and a small shotgun. The contents of the capacitor, about a cubic centimeter 
of plastic fibers, hit the student in the middle of the chest. He was not injured, but for a brief 
moment he thought he had been killed by this exploding capacitor. I suspect he still pays 
more attention to capacitor voltage ratings than most people. 


Many coilers recommend placing resistors across the MMC to bleed off the charge after 
power is removed. This can eliminate some very unpleasant surprises when making adjust- 
ments between runs. One can get true high voltage resistors but they are expensive. Most use 
a long chain of ordinary resistors therefore. Digi-Key sells 0.5 W carbon film resistors rated 
at 350 volts each for about two cents each in quantity. At this price, there is little point in 
not using adequate safety factors. 


Suppose we feel comfortable in operating a 0.5 W resistor at 0.1 W and at the rated peak 
voltage of 350 V each. We solve for the resistance as 


where we round off to the nearest standard value of 620 kQ. 


The number of resistors needed for a string across say a 15 kV transformer secondary 
would be 


15000V 2 
N- 150002 


350 = 60 resistors 


We then check the time constant to get a measure of how long it takes to discharge a 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 3—Lossy Capacitors 3-12 


Table 3.3: Dielectric Constant ¢, and Dissipation Factor DF of Water. 


Frequency in Hz 


Temperature 10° 10° 10 107° 
1.5°C é 87.0 87.0 87 38 
DF 0.1897 0.01897 0.00195 1.026 
25° é 78.2 78.2 78.2 55 
DF 0.3964 0.03964 0.00460 0.545 
85°C é 58 58 58 54 


DF 1.2413 0.12413 0.01259 0.259 


capacitor chain. If the capacitance happened to be 27 nF, the time constant would be 


T = RC = (60)(620000)(27 x 10~°) = 1 second 


A capacitor will be mostly discharged after five time constants or 5 seconds in this example. 


There are obviously many other design possibilities. If we assumed the resistors would 
accept a higher voltage before arcing over, then we would want to use higher ohm values to 
keep heating within bounds. Terry Fritz tested some Yageo 10 MQ 0.5 W carbon film resistors 
to see where they would actually fail. At 4000 V the resistor started to turn brown and smoke 
a little. This should be expected since it is now drawing 2 W, four times its power rating. At 
5300 V, it gave a very satisfying crack and arced along its outside surface. After that, the 2 
MQ resistor measured 11.1 MQ. So a resistor specified at 350 V actually failed at 5300 V, a 
safety factor of 15. Terry sees no problem with using these 0.5 W resistors up to 1000 V if 
the resulting wattage is not excessive. 


3.2 Dielectric Loss in Water 


Water did not appear in Table 3.1 because it is so much different from other substances that it 
needs special treatment. The relative permittivity ¢/, and the dissipation factor DF are shown 
in Table 3.3 for different values of temperature and frequency. 


The columns for 10° and 10° Hz (100 kHz and 1 MHz) would bracket the range of operation 
for most Tesla coils. Only the very large coils operate below 100 kHz and the very small coils 
operate above 1 MHz. The column for 10” Hz was included to show that relative permittivity 
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does not change between 10° and 10’ Hz and that the dissipation factor continues the same 
decline as between 10° and 10° Hz. As frequency increases to 101° Hz (the general range of 
microwave ovens) the relative permittivity decreases somewhat while the dissipation factor 
becomes quite large. It is obvious that heating food with radio waves is much more effective 
at 10!° Hz than at 107 Hz. 


At Tesla coil frequencies, the dissipation factor (DF), at some frequency fz can be found 
from the listed value (DF), at frequency f; by 


(DF), = (DF), 2 (3.14) 


3.3. Dielectric Losses of a Tesla Coil 


There are dielectric losses primarily in two places in a Tesla coil, the coil form and the 
insulation around the conductor. There may also be some loss in the air surrounding the coil, 
in the soil, and in insulating supports near the coil, but the ones of most interest seem to be 
the coil form and the wire insulation. We saw in Eq. 3.11 that the loss in a capacitor is equal 
to the square of the applied voltage times the radian frequency times the capacitance times 
the dissipation factor. We can make reasonable estimates for voltage, frequency, and DF, but 
what about the capacitance? 


A cross section of the coil is shown in Fig. 3.3. 


Vior ae Dy o= = 


Figure 3.3: Cross Section of Tesla Coil 


We show three turns of the several hundred that are typically used. This is a close wound 
coil where the conductors are immediately adjacent to one another. Space wound coils usually 
use a spacer of the same size as the conductor, which would result in a cross section showing 
half as many conductors in a given winding length. The coil diameter is D, the winding length 
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is ¢,, and the coil form length is ¢.¢. The wall thickness of the coil form is t. At the top of 
the coil is a toroid with voltage V;o,, while at the bottom the voltage is zero. To make life a 
little simpler, we assume the toroid can be represented by a flat sheet on top the coil form, 
and that the zero potential ground plane is directly under the coil. 


Most of the dielectric between toroid and ground plane is air, but a small portion is the 
coil form. We can think of the total capacitance between toroid and ground plane as a coil 
form capacitance C., in parallel with an air dielectric capacitor and a winding insulation 
capacitance C,,;. If the parallel plate capacitor formula is valid, we can quickly write an 
expression for the coil form capacitance. 


emt Def 


Cop = (3.15) 


boy 


As stated earlier, the expression is valid when fringing can be ignored. This condition is 
true when the separation is small compared with the area but here we have a small area and 
a large separation. So we go back and examine the situation for which fringing is negligible. 
In any capacitor, there are electric field lines directed from the positive to the negative plate, 
and equipotential lines (lines or surfaces of constant potential) which intersect each other at 
right angles. The areas enclosed by these lines are called curvilinear squares. The sides may 
be curved and of different lengths, but the four corners are all right angles. In a parallel plate 
capacitor, the equipotential lines are equally spaced straight lines, as are the electric field 
lines, so the curvilinear squares become true squares. It turns out that inside the coil form of 
Fig. 3.3 that electric field lines are mostly straight down and the equipotential lines are nearly 
equally spaced. 


For the equipotential lines to be exactly equally spaced, the voltage must increase linearly 
from bottom to top. The actual mathematical form of the increase is probably not quite 
linear, as shown when the turn-to-turn voltage of the top few turns, for example, exceeds 
the dielectric strength of the insulation and turn-to-turn sparking occurs. But it will not be 
vastly different from linear, either, so the parallel plate formula should give the proper order 
of magnitude. I would be surprised if it were in error by more than 20% or 30%. 


A circuit model for the equivalent resistances that represent losses in the coil form and in 
the winding insulation is shown in Fig. 3.4. The capacitances Cyr and Ci; are in parallel with 
other capacitances (not shown) such that the total capacitance is C;,. The power dissipated 
in the coil form is given by 


Pap = VeqeCop(DF) cy (3.16) 
and the equivalent resistance is given by 
Vie 1 


Ret = ———_ 3.17 
T° "Peg wCeg(DF)ey nt) 
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Example: What is the power loss and equivalent coil form resistance in a PVC coil form 
that is 1 m in length, 3 mm wall thickness, and diameter of 200 mm, at a frequency of 200 
kHz if the top voltage is 500,000 V? Assume e, = 3.0 and DF = 0.015. 


€o€rDept _ (8.854 x 107 '*)(3.0) (0.200) (0.003) 


=5.01x10-4 F 
lef 1 . 


Cop = 


Po = Vi2,wCes (DF) = (500, 000)27r(200, 000) (5.01 x 107 4(0.015) = 236 W 


1 1 


= ee Rie ee 
wCcs(DF).¢  2(200,000)(5.01 x 10-14) (0.015) = 


Rez 


To get the toroid to this voltage probably requires an input power of 10 kW or more, so 
236 W is not a large loss, percentage wise. On the other hand, the coil form would soon melt 
if this power were applied continuously. If maximum efficiency is desired, one can always go 
to a polyethylene coil form, where DF is less than 0.0002 


Vior 


Figure 3.4: Equivalent Resistances for Coil Form and Wire Insulation Losses. 


We turn now to the consideration of the dielectric loss in the winding insulation. Most Tesla 
coils are wound with magnet wire, a wire in wide use for winding motors and transformers. 
It has a dielectric coating that is tough both mechanically and thermally. It will withstand 
temperatures well above 100°C. Unfortunately, it is not particularly low loss at Tesla coil 
frequencies. Dissipation factors are similar to PVC or nylon. The coating called Heavy 
Soderon actually has a thin layer of nylon on the surface, to help withstand mechanical abuse. 


I have observed that the losses of tight wound magnet wire coils sometimes increase dra- 
matically with humidity, even at input voltages far below those necessary for breakout. It 
seems that some magnet wire coatings can absorb moisture, and also some coil forms. I 
bought two plastic barrels at a recycling place that I thought were polyethylene. Whatever 
they were, the input resistance of the coil could change by a factor of two or more as humidity 
changed from 25 to 100%. On the other hand, a coil of some magnet wire on a PVC form 
would see little change in input resistance, probably less than 10%. So moisture in or on the 
wire insulation and the coil form certainly has the capability of increasing losses. 
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Figure 3.5: Two Adjacent Turns of Magnet Wire 


The basic geometry that we will try to solve for power loss is shown in Fig. 3.5. 


This figure shows two adjacent turns of a space wound coil. The conductor has radius 
r, and is surrounded by a dielectric of permittivity «;. The turns are then contained in a 
dielectric of permittivity €2, where usually €2 = €,. However, if the winding has been coated 
with some substance like polyurethane, then €9 will be different from that of air. 


Like many problems in electromagnetics, there is no exact analytic formula for capacitance 
of the exact structure. We have a choice of an exact solution for an approximate structure, 
or an approximate solution for an exact structure, or combinations thereof. Fig. 3.5 shows 
the parameters for a structure that approximates our coil winding, and has an exact analytic 
solution. The advantage of an analytic solution is that we can easily see which factors are the 
most important. Questions about the effect of wire size, insulation thickness, and conductor 
spacing can probably be answered correctly even if the answers differ from the correct values 
by 20% or more. 


The geometry in Fig. 3.5 is for a parallel conductor transmission line. The conductor 
on the right is at potential V = V, and the one at the left at potential V = —V,. The 
vertical plane between the conductors is then at potential V = 0, by symmetry. The potential 
is V = V, (a constant) on the outside of the dielectric layer. If the dielectric is eccentric, 
then the surfaces lie on the surfaces of constant 7 of the bycylindrical coordinate system [2, 
pages 361-366]. There are three factors that make this an approximate geometry: First, the 
actual dielectric is centered rather than eccentric. Second, the magnet wire will lie against 
a coil form, which will affect the electric field lines and the capacitance. And third, there 
are additional conductors on either side of the two shown, which also affect field lines and 
capacitance. However, I think this geometry is the best that we can do, if we want an analytic 
solution. 


The capacitance between the conductor on the right and the outer surface of its dielectric 
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Qrel 
C= (3.18) 
No — Nea 


and the capacitance between the outer surface and the V = 0 plane is 


= 21 ent 
Ne 
where @ is the circumference of one turn of the winding. To get the power dissipation, we 


need the voltage V, — V, across €, and the voltage V, across €2. The solution process is fairly 
standard for Laplace’s Equation. We write general solutions for Laplace’s Equation in both 


C2 (3.19) 


dielectrics. 


Vi = Ay + Bin in € (3.20) 


Vo = Ao t+ Bon in €9 (3.21) 


The boundary equations are 


V=V at. =e (3.22) 
y=0 at 1=0 (3.23) 
Vi = V2 at: 7 = He (3.24) 


The fourth boundary equation is that the normal electric flux density is continuous on the 
7 = Nx surface. 


OV, OV2 
et ee 3.25 
€1 On €2 On ( ) 
After some algebra, we find 
Yj = Vata ee) eat (3.26) 
No€2 + Nx(€1 — €2) 
Voerm (3.27) 


2 ieee EEE 
No€2 + Nx(€1 — €2) 
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We can find V, by substituting 7 = 7, in either of the last two equations. The power 
dissipations in the two regions are 


Ve (No = Nx )w(27e1€)(DF)1 


P, = (Vo —Vx)?wC1 (DF), = 
1=( POOP E) [No€2 + Ne(€1 — €2)]? 


(3.28) 


_ Vee ngw(2rel)(DF)2 


Py = V;wO2(DF)2 = Wa Dnte soir (3.29) 


The voltage V, is half the voltage between turns. If there are N turns and the voltage at 
the top of the coil is Vio, then 


= Vior 
2N 


Vo (3.30) 


The total power P,,; dissipated in the wire insulation and in the gap between the turns is 
The relationships between the variables in Fig. 3.5 are 


w= \ar+r? (3.32) 


where i = 1 or 2. 
No = sinh7! (=) (3.33) 
Tr 


Ne = sinh! (=) (3.34) 
r 


s=w2.—w, = 1a? + re —\/a+rZ (3.35) 


We start with a table lookup to find the radius r; of the bare magnet wire (without 
insulation). Then we use w, = 71+ build + space (if any), where build refers to the thickness 
of the dielectric on the magnet wire, and solve for a and 7,. Then it gets a little tricky. Do 
we use an ro such that the minimum distance between the eccentric cylinders is the build, 
or the average distance, or something in between? My first impulse is to use the minimum 
distance, since most of the ‘action’ is between the two conductors, so we need the most accurate 
description of dielectric in that region. That is, I would write 
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rg —s =17, + build (3.36) 


This is a transendential equation in one unknown, rg. There is no simple closed form 
solution, but if one plugs in values for a and rj, it is not difficult to get a numeric value for ro. 
Note that if the dielectrics are touching, then ro, s, and we all go to oo. The model requires 
that P; = 0 for this case and P, is the dissipation in e€, if €, fills the entire space. 


Example: Find the dielectric losses in three different coils built on the same coil form used in the 
previous example. The winding length is one meter, the coil diameter is 200 mm, the top voltage is 
500,000 V, and the toroid size is adjusted as necessary to get a frequency of 200 kHz. We want to 
compare the cases: (a) a 16 gauge coil tight wound, (b) a 16 gauge coil with a small air gap, 2 mils, 
between windings, and (c) a 22 gauge coil space wound with the same number of turns as the 16 gauge 
tight wound coil. These cases all have dry air as the surrounding dielectric. Then we want to examine 
a fourth case of the same geometry as (b) but with water condensed on the magnet wire surface. We 
are a long way from having a good model for this situation, but we might get some interesting results 
if we assume that the 2 mil gap is filled with 10% water and 90% air and use interpolated values for 
permittivity and dissipation factor. The magnet wire dielectric is Heavy Soderon with thickness 1.65 
mils for either wire size. Results are shown in Table 3.4. 


Table 3.4: Dielectric Loss Example 


wire size 16 16 22 16 

r 1, mils 25.4 25.4 12.65 25.4 
space, mils 0 1 12.85 1 
wy, mils 27.05 28.05 27.15 28.05 
N, turns 725 700 725 700 
L, mH 19.0 17.7 19.0 17.7 
a, mils 9.303 11.901 24.023 11.901 
Ep] 3.39 3.39 3.39 3.39 
Ep 1 1 1 7.82 
(DF), 0.03 0.03 0.03 0.03 
(DF)»2 0 0 0 0.198 
r9 OO 70.32 16.03 70.32 
No 0.3585 0.4529 1.400 0.4529 
Ne 0 0.1684 1.194 0.1684 
Vo 345 357 345 357 
Pi 2.36 0.352 0.010 2.017 
Py 0 0 0 3.416 
2NP, 3423 493 14 2824 
2N P» 0 0 0 4782 


We see a loss of 3423 W for the first case. This seems high, but the input power to the coil at 
the moment the top voltage is 500 kV is probably at least 10 kW and perhaps well over 100 kW, so 
this is not impossible. We also note that small winding irregularities lower this number significantly. 
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If the average gap between windings is only 2 mils, the loss drops from 3423 to 493 W. It is difficult 
for amateurs to wind coils with no gap at all, so the effective loss is most likely under 1 kW. If we go 
to a true space wound coil, as for in the 22 gauge shown in the third column, the losses drop to only 
14 watts. If we planned to operate the coil in CW mode, or for very extended periods in disruptive 
mode, we may want to consider space winding just to keep the coil temperature rise low. 


When 10% water is added to the dielectrics for the case of a 2 mil gap, the losses increase dramat- 
ically. They jump from 493 W to 2824 + 4782 = 7606 W. These are not precise numbers, but indicate 
the trend that moisture increases coil losses substantially. 


3.4 DC Capacitors 


In this section we will examine the dc capacitors used to convert time varying voltages into 
steady dc, called filter capacitors. They usually only experience one polarity of voltage, hence 
can be of the polarized type. Aluminum oxide electrolytic and tantalum capacitors are of 
this type. The exact value of a capacitor used for filtering is not very critical. One may see 
tolerances as high as —20% to +80%. That is, a 100 uF capacitor would have an actual value 
between 80 and 180 wF with this tolerance. Similarly, capacitance variation with temperature 
is not critical. 


The most important goal is a high capacitance in a small volume. We accomplish this 
by using a material with a high dielectric constant in a structure with a large area and thin 
dielectrics. One ingenious solution to this goal is the aluminum electrolytic capacitor. This 
consists of a ribbon of aluminum foil on which a thin film of aluminum oxide has been formed 
electrochemically, and a water based electrolyte fluid which acts as the opposing plate. A 
second ribbon of aluminum foil is used to make electrical connection to the electrolyte fluid. 
The two ribbons of aluminum foil are mechanically separated by a ribbon of porous paper, as 
shown in Fig. 3.6. The three ribbons are rolled up and placed in a can filled with the liquid 
electrolyte. 


spacer and electrolyte (cathode) 


— + 


cathode film anode film 


oxide 


Figure 3.6: Cross Section of Aluminum Electrolytic Capacitor 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 3—Lossy Capacitors 3-21 


The oxide coated foil is the positive plate, called the anode. The aluminum oxide film is 
the dielectric, and the fluid electrolyte is the negative plate, called the cathode. The second 
aluminum ribbon serves only to make good electrical connection to the cathode. It is usually 
bonded to the aluminum can that houses the capacitor. 


The oxide dielectric has a thickness on the order of 0.01 um. The plate separation d has 
therefore been reduced to a very small value. Also, aluminum oxide has a dielectric constant 
in the range of 8 to 11, which is relatively high for a dielectric. This combination allows a 
very high capacitance for a small volume. 


The oxide is grown by placing the aluminum foil in an electrolyte solution and applying a dc 
voltage between the aluminum and the electrolyte. The resulting oxide layer has a thickness 
proportional to the applied voltage, about 0.0014 wm per volt at room temperature. The 
voltage at which the oxide is grown must be considerably higher than the proposed operating 
voltage to provide adequate dielectric strength. The leakage current increases rapidly as the 
operating voltage approaches the forming voltage value. 


The film of aluminum oxide has the property of diode action. That is, current can flow one 
way through the anode but not the other. The aluminum electrolytic is therefore limited to 
dc applications. The polarity marked on the can must be carefully observed since a negative 
voltage more than a volt or two will cause breakdown of the film and destruction of the 
capacitor. This limitation is overcome in non-polarized electrolytic capacitors intended for ac 
applications by simply using aluminum oxide layers on both aluminum ribbons and facing the 
ribbons in opposite directions. This produces the same effect as putting two diodes in series, 
but pointed in opposite directions. 


While the capacitance per unit volume is quite high for this capacitor type, there are 
several disadvantages. One is that the maximum operating voltage is limited to about 450 
V. This is obviously not a problem for most computer applications, but can be a limitation 
in situations requiring high power. Another limitation is that the dissipation factor is higher 
than for other capacitor types. A third limitation is that there may be significant leakage 
currents. That is, both R, and R, need to be considered in doing loss calculations. A fourth 
limitation is that the shelf life is somewhat limited. A capacitor with an acceptable leakage 
current when first made may have an unacceptable leakage current after not being used for a 
year or two. There is a process called reforming where the capacitors are heated to 85°C and 
a voltage is applied, which can restore acceptable performance when this problem occurs. 


Tantalum capacitors usually have solid electrolytes, are smaller, and have better char- 
acteristics than aluminum capacitors. The small size of tantalum capacitors is due to the 
high dielectric constant of the tantalum-oxide dielectric. In addition to being small, tantalum 
capacitors have relatively low leakage currents, good resistance to vibrations, minimal capaci- 
tance and ESR variations with changes in temperature and no major parametric changes over 
the life of the device. 


It is interesting to compare a high quality ac capacitor with a Tantalum. According to 
Table 3.2, a 1 uF 100 V polypropylene capacitor has a value of R, = 0.015 QO. A MEPCO 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 3—Lossy Capacitors 3-22 


Series 40YW Tantalum capacitor rated at 1 wF and 50 VDC has a typical R, of 16 Q, and 
is still within specification at 80 Q. This does not mean the Tantalum capacitor is inferior to 
the metallized polypropylene capacitor, but rather that the Tantalum capacitor should not 
be used in applications requiring large ripple currents. 


Example 


Calculate the power dissipation in a capacitor described by R, = 0.1 2, Rp = 10°Q, L, = 10uH, 
and C = 400uF if the applied voltage is v = 200 + 6sinwt, with w = 27(120) = 744 rad/sec. This 
would describe a 200 VDC supply with a 6 V ripple voltage. 


We need to apply superposition to solve for the dc power dissipation and then the ac power 
dissipation. At dc, C represents an infinite impedance, so we have R,, L,, and R, in series. This series 
impedance is essentially just R,. The average power dissipated in R, is 


v2, (200)? 
P. c= de GCC 0.04W 
a fp. TDS 


At 120 Hz, the capacitive reactance is a few ohms in parallel with R, = 10°Q, so the capacitor 
model essentially becomes R,, L,, and C in series. The series impedance of this circuit at 120 Hz is 


Vie Oe ee =a = 0.1 + j0.00744 — 73.316 
Ww 


with magnitude |Z,-| = 3.310. Obviously, the capacitive reactance would have been an excellent 
approximation to the total series impedance. The rms current is then 


6 
Tac = ——— = 1.28A 
V2(3.31) 


The average power due to this current flowing through R, is 
Pye = I2,.Rs = 0.16W 
The total power being dissipated in the capacitor is then the sum of Pg. and Pye. 
Prot = Pac + Pac = 0.04+ 0.16 = 0.2W 


We see that the ac ripple of only 6 volts causes four times the heating of the 200 volts of dc. Also 
the capacitor leads must be built to carry a current of 1.28 A. This would not be a problem, but if the 
capacitance were increased by a factor of ten, wire size would start to become important. 


3.5 Pulse Capacitors 


Applications for pulse capacitors include particle accelerators, metal forming, laser drivers, and 
X-ray generators (and Tesla coils, of course). In these situations, discharges are often highly 
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oscillatory, and voltage reversals up to 90% can be seen. Flash tubes and electrical impulse 
welders have relatively longer discharge times and smaller reversals. Pulse-discharge or energy- 
storage capacitors are usually charged over a relatively long time and discharged in a short time 
(microseconds or less). Voltages and currents may be very high. The inexpensive aluminum 
and tantalum capacitors do not meet these requirements, so special techniques are used to 
fabricate these capacitors. The traditional materials used include castor oil/Kraft paper, 
plastic film, ceramic, and mica, although the polypropylene capacitors considered earlier in 
the chapter would also work. This section is included to help coilers determine the suitability 
of surplus pulse capacitors for Tesla coil use. 


An important parameter is the pulse repetition rate. If this is faster than, say, 1 pulse/minute, 
heating effects become important. Overheating can shorten the lifetime of a pulse capacitor 
considerably. 


The lifetime is typically expressed in number of pulses or shots. Capacitors are readily 
available with a lifetime of greater than 100,000 shots, with voltages up to 100,000 V and 
peak currents up to 50,000 A, if the peak current of the first negative peak of the decaying 
sinusoid is no more than 30% of the first positive peak. 


There are some rather interesting limits to the speed of discharge. One is that the capacitor 
plates act like a parallel plate transmission line. A charged transmission can only discharge 
in a time equal to the two-way transit time, t, along the line where 


= ae (3.37) 
where @ is the greatest distance from the capacitor lead to the end of the plate, and c is the 
speed of light in vacuum. Light propagates at a speed of 0.3 m/ns in vacuum so it is quite 
possible to exceed 100 ns discharge time with a rolled parallel plate construction. A capacitor 
built with 0.1 wF of capacitance and 10 nH of inductance will not ring down in 100 ns because 
of this effect, even if that is the value predicted from the (low-frequency) circuit theory model. 
The only way to make a fast capacitor is to manufacture it with short sections such that this 
transmission line transient is much shorter than the desired lumped element discharge time. 


This transmission line discharge is a good way to get a single pulse of carefully controlled 
voltage and period. A coaxial cable is charged to a desired voltage and then discharged 
through the load, say with a mercury switch. This is an inexpensive method of getting a pulse 
of say 500 V and 40 ns. 
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INDUCTORS AND TRANSFORMERS 


We now consider the inductor. Like the capacitor, it appears in all sizes. There are signifi- 
cant differences, however. One is that the commercial availability of prefabricated inductors is 
much less than for capacitors. Engineers dealing with inductors will routinely buy coil forms 
and wire, and assemble their own device. With capacitors we can get by with only a vague 
notion of how they were designed and built. With inductors, we must do the design, which 
forces us to have a deeper understanding of the steps involved. 


A second difference is the amount of electromagnetic theory required. With capacitors, we 
could escape with a mention of electric field and permittivity. Inductors require that we jump 
right into some challenging concepts of magnetic fields and energy, and even set up some line 
or volume integrals. We will refer back to a first course in electromagnetic theory as needed, 
and even pull out some results from more advanced courses. 


4.1 Definitions 


Consider a coil of wire as shown in Fig. 4.1. The resistance of the wire can be modeled as 
a separate lumped device so we can think of the coil as being perfectly conducting. If i is a 
finite dc current, the voltage v will be zero, as would be expected across a perfect conductor. 
If a time varying current is applied, however, there will be a related voltage observed across 
the coil. From the circuit theory viewpoint, the relation is given by 


(4.1) 


Figure 4.1: Simple Inductor 


We will also observe a time varying voltage when we pass a magnet by the coil, even if the 
current is zero. (Faraday was the first to observe this, in 1831). Moving a magnet is philo- 
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sophically quite different from applying a current, but it turns out that we can mathematically 
describe both situations with the electromagnetic equivalent of Eqn. 4.1. 


v= N— (4.2) 


where N is the number of turns on the coil and ¢ is the magnetic flux passing through the 
coil. The direction of a flux ¢ that is produced by a current 7 is determined by the right hand 
rule. That is, if you curl the fingers of your right hand in the direction of the current flow, 
the thumb will point in the direction of the flux. 


Setting Eq. 4.1 and Eq. 4.2 equal to each other and integrating to remove the differential 
operator yields 


Li=N¢ (4.3) 


This equation has circuit quantities on the left and field quantities on the right, allowing us 
to move back and forth between the two ways of thinking. Solving for the inductance L gives 


op x 
[Sees 4.4 
ae (4.4) 
where we have introduced the flux linkage X. This term better describes the case where ¢ is 
not constant between adjacent turns. The flux linkage can be considered as the equivalent 


flux which gives all the correct results when passing through a single turn coil. 


The flux ¢ and the flux linkage ’ are proportional to the current i. The relationship is 
linear if there are no ferromagnetic materials in the vicinity, which then gives us a constant 
value for L independent of the actual value of 7. Thus for the air-cored coil, L is just a function 
of the geometry of the coil, much like our expression for capacitance that was calculated from 
area and separation of plates in the previous chapter. Unfortunately, inductance formulas 
tend to be much more complicated than the formula for a parallel plate capacitor. 


The electric power input to the inductor is 
de 


p=v= Ne (4.5) 


There are no losses in our perfectly conducting coil, so whatever power flows in at one time 
must flow out at another time. In the sinusoidal case, power flows in for half a cycle and back 
out the next half cycle. The power flow results in stored magnetic field energy in the coil. 
The differential energy input during the differential time dt is 


dW = pdt = Nidd (4.6) 


where W is the stored energy in the field. 
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We want to relate this stored energy to the field quantities B and H, where 
B= pippalt - (4.7) 


B is the magnetic flux density in tesla (T) or webers/m? (Wb/m7?), H is the magnetic intensity 
in A/m, p, is the relative permeability (= 1 for vacuum), and jg is the permeability of free 
space in henrys per meter. 


Mo =4n x 107" H/m (4.8) 


The relative permeability is very close to unity for all materials except for the ferromagnetic 
materials iron, cobalt, nickel, and a number of special alloys. For these materials, 4, may 
range from 10 to 10°. The relative permeability is also a function of magnetic intensity in 
ferromagnetic materials, making what would be a linear problem into a nonlinear one. 


The magnetic flux density B may also be expressed in gauss, where 10* gauss = 1 tesla. 
The earth’s magnetic flux density varies from 0.2 to 0.6 gauss, depending on location. The 
60 Hz magnetic flux density in a home or office is usually less than a few milligauss except 
near a source (electric heater, computer monitor, electric razor, blow dryer, etc.) where 
it may be a few tens of milligauss or even more. A modern well-designed 60 Hz power 
transformer will probably have a magnetic flux density between 1 and 2 T inside the core. It 
requires considerable effort and special designs to get much above 2 T. The necessary current 
density causes heating in the conductors, unless, of course, the conductors are cooled into the 
superconductor region. Fluxes as high as 8 to 16 T have been used in accelerators and energy 
storage systems. 


Other conversion factors which might be needed are: 


e 1 Oersted = 250/7 = 79.6 ampere-turns/meter 


e 1 Tesla = 10,000 gauss 


The flux ¢ passing through an area A is the integral of the magnetic flux density B over 
that area. 


o= | Baa Wb (4.9) 
which becomes simply ¢ = BA if B is constant over the area. 


We also need Ampere’s circuital law 
Ni= pHa (4.10) 
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This states that the current enclosed by any arbitrary path is given by the integral of the 
dot product of the vector H and a differential length dé along that path. If we extend our 
coil around into a toroid shape, H (the magnitude of H) will be essentially constant inside 
the toroid and Ampere’s circuital law becomes 


Ni=Hé (4.11) 


where @ is the length of a circle in the toroid. The energy stored in the magnetic field is now 


H 
iW Nrib= 2 age Wena (4.12) 
7 


where Vol = A@ is the volume where the magnetic energy is stored. The total energy can be 
found by integration. 


B 
Wes (Vol) [ HdB Joules (4.13) 
0 


If the permeability is constant (the magnetic circuit is linear) the integral can be quickly 
evaluated. 


(Vol) B? H? 

W= TE Ta (Vol) Hr o> J (4.14) 
This equation has much important information in it. Suppose that we have a magnetic circuit 
that is entirely ferromagnetic. H is determined by the current and is independent of the 
permeability. B = wH is large and the total energy stored is large. Suppose now that we 
cut a small air gap across the magnetic circuit. The flux @ drops substantially because of 
the increased reluctance of the magnetic circuit. B will have about the same (smaller) value 
in both the iron and the air gap, so H = B/w will be much larger in the air gap because 
of the lower permeability. The total integral of H-dé stays the same but a large fraction 
of the integral comes from the air gap portion. So the total energy stored decreases as the 
air gap length increases, and the fraction of the total energy stored in the air gap increases 
dramatically. 


Although not as instructive, the total energy can also be expressed in circuit quantities as 


1 
W= st J (4.15) 


We can also use the result from Ampere’s circuital law to determine the flux ¢. In the 
simple case of uniform flux density B and no air gap, this becomes 


[bp Lo Ni A 


b= BA = prtgH A = PPX 


(4.16) 
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4.2 Ferromagnetic Losses 


Two things happen when a time varying magnetic field exists inside a ferromagnetic material. 
Magnetic domains rotate in the material to align with the magnetic field and the Faraday 
voltage induced inside the material produces what are called eddy currents. The rotation of 
domains each cycle produces a frictional type loss called the hysteresis loss P,. Experimentally, 
we find that 


Pr = Kn Bias (4.17) 


where kK}, is an empirical property of the material, f is the frequency, Baz is the maximum 
flux density, and z is an empirically determined value, usually between 1.6 and 2.0 for power 
frequency transformer steels. 


The hysteresis loop of the material is obtained by plotting the magnetic flux density B 
against the magnetic intensity H as shown in Fig. 4.2. The area inside the hysteresis loop 
is the energy dissipated each cycle. Some ferromagnetic materials have very thin hysteresis 
loops, resulting in low losses, while others have relatively fat loops and correspondingly high 
losses. 


He aan 


Figure 4.2: Hysteresis Curve 


As the driving magnetic field increases, the resulting flux density increases at a slower 
rate as more domains become aligned with the magnetic field. This phenomenon is called 
saturation. If linearity is desired, then the transformer should be operated at low flux levels 
where the hysteresis loop is nearly linear. For power transfer, however, it is most cost effective 
to operate the device into its saturated region. The exact amount is a matter of engineering 
judgment. 
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We might define at least two definitions for permeability, from which we can get some 
guidelines for saturation. These are the dc and ac permeabilities 


B 
Mde = FF (4.18) 
AB 
2c — {aS 4.19 
Hac = Ry (4.19) 


The two permeabilities are identical for very low drive levels that are symmetric about zero. 
If a de bias exists, the ac permeability will always be smaller than jug.. The ac permeability 
is the main parameter of interest to filter choke designers. At least some engineers consider 
a material to be saturated when the dc permeability has dropped to half its initial value, or 
when the ac permeability has dropped to one-eighth of its initial value [12, page 26,28]. 


A magnetic circuit with eddy currents is shown in Fig. 4.38. The current is inversely 
proportional to the resistance seen by the induced Faraday voltage. In a large piece of steel 
the resistance can be very low, even though the resistivity of steel is not very low compared 
with a good conductor like copper. For this reason, magnetic circuits at power frequencies 
are usually made of thin sheets of steel, called laminations, which are on the order of 0.5 mm 
thick. 


eddy current 


Figure 4.3: Eddy Currents in Magnetic Circuit 


The equation for eddy current loss has the form 


PS kp B (4.20) 
The experimentally determined constant kK, depends on the resistivity and the dimensions of 


the material. A detailed analysis shows that K, is proportional to the square of the lamination 
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thickness, so it is important to keep the laminations as thin as possible. Eddy current losses 
can be kept acceptably low at 60 Hz with little difficulty, but become excessive at a few kHz, 
even with very thin laminations. Therefore, inductors or transformers built for operation 
above 1 kHz are rarely made of laminated material. Instead, they are made of even smaller 
pieces of ferromagnetic material, typically powdered iron or ferrites. Powdered iron suffers 
from low permeability and low resistivity compared with ferrites, so we shall concentrate on 
the latter. 


4.3 Ferrites 


Ferrites were developed during and after World War II. The chemical formula for ferrites 
is ZFeg O4, where Z stands for any of the divalent ions: zinc, copper, nickel, iron, cobalt, 
manganese, or magnesium, or a mixture of these ions. The bulk resistivities are in the range of 
10? to 10° ohm-cm, compared with 10~° ohm-cm for powdered iron. This very high resistivity 
reduces the eddy current losses so that ferrites can be used for frequencies up to 20 MHz or 
even more. 


Ferrites are basically a form of ceramic, made by mixing fine powders of appropriate oxides, 
compressing the mixture, and firing it in carefully controlled atmospheres at temperatures of 
about 1100°C to 1200°C. The most common ferrites are mixtures of two ferrite powders, either 
manganese-zinc or nickel-zinc. Magnetic properties can be varied over a significant range by 
changing the ratio of the two divalent ions and by changing the processing conditions. Each 
material is given a unique code number. Ferroxcube, for example, assigns a “3” as the first 
digit of its MnZn materials (3C85, 3B7, 3D3, etc.) anda “4” as the first digit of NiZn materials 
(4C4, 4A, 4A6, etc.) 


4.3.1 Ferrite Temperature Limits 


The Curie temperature (the temperature at which the ferrite becomes nonmagnetic) can be 
relatively low. The Ferroxcube 3E5 ferrite may have a Curie temperature as low as 120°C. 
Other Ferroxcube materials have Curie temperatures up to 300°C. If there is any possibility 
of the ferrite device operating at a temperature above 120°C due to internal losses, a ferrite 
material with an adequate Curie temperature must be selected. 


The copper wire used in winding inductors or transformers remains mechanically stable 
at temperatures far above the ferrite Curie temperature, so the wire itself is not of concern. 
However, the insulation on the wire may fail at relatively low temperatures. Polyvinyl chloride 
(PVC) insulated wire typically has a maximum temperature rating between 80°C and 105°C, 
for example. Magnet wire has a somewhat higher rated temperature, such as the Belden poly- 
thermaleze coating rated at 180°C. Belden also makes Teflon coated wires rated at 200°C and 
260°C. 
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Heat is dissipated from the surface of the inductor or transformer by a combination of ra- 
diation and convection. Heat radiated depends on the device surroundings, while convection 
depends on air flow over the device, so it is very difficult to accurately predict temperature 
rise in most installations. In any case, both radiation and convection will be directly propor- 
tional to the total exposed surface area of the core and windings. We can therefore describe 
the independent variable as total watts dissipated in core and copper per unit area of the 
device. The dependent variable, temperature rise, is directly proportional to the dissipation 
in W/cm?. One manufacturer (Magnetics) has calculated a temperature rise of 10°C for a 
surface dissipation of 0.01 W/cm? and a rise of 100°C for a surface dissipation of 0.1 W/cm?, 
given some reasonable assumptions. Each increment of 0.01 W/cm? results in a temperature 
rise of 10°C. 


Example. 


An inductor has a total heat dissipation of 0.06 W/cm? and is in an ambient temperature of 50°C. 
What is a reasonable estimate of inductor temperature? 


Based on the Magnetics guideline, 0.06 W/cm? should yield a 60°C temperature rise above the 
ambient, so the inductor temperature will be 60 + 50 = 110°C. PVC insulated wire should not be used 
in this situation. 


4.4 Mutual Inductance 


Consider two inductively coupled coils as shown in Fig. 4.4. The current 7; produces a flux 
1, that links with 71. Part of ¢11 is lost as leakage flux ¢y¢ and part of it, flux ¢2) links both 
currents 71 and ig. Current 72 likewise produces a flux ¢22 with part of it, flux $12, that links 
both currents. The relationship among these fluxes is 


bir = gat bie (4.21) 
b22 = P12 + G20 (4.22) 
The self-inductance of circuit 1 is 
N- 
Ly = Nout (4.23) 
a4 


and similarly for £29. The mutual inductance of circuit 1 with respect to circuit 2 is based 
on the flux in circuit 1 that is produced by the current in circuit 2. 


- Nid¢12 
19 


Dis (4.24) 


and similarly for D9. 
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Figure 4.4: Mutual Inductance 


It can be shown that Dy = L2; in a homogeneous medium of constant permeability. To 
emphasize this fact we define a new symbol M for the mutual inductance. 


M = [yo = Lo (4.25) 
The maximum value for M is VL4,L22. We define the coefficient of coupling k as 


ee! (4.26) 
VE Lop 


The coefficient of coupling can reach values as high as 0.998 in iron-core transformers. It 
is difficult to make k much above 0.5 in air-core transformers. 


The voltage v2 produced by the primary current 7; is given by 


vo = M— (4.27) 


If the two inductors of Fig. 4.4 are connected in series, the total inductance is 


DL= 14, + Log 42M (4.28) 


where the plus or minus is determined by whether the mutual flux tends to reinforce or 
cancel the fluxes of the individual coils. This is a convenient method to measure the mutual 
inductance. Just measure the series inductance twice, once with one coil reversed, subtract 
one result from the other, and solve the resulting expression for M. 
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4.5 Inductance Formulas 


Let us now examine the inductance formulas for some simple geometries. First we will look at 
the inductance of a nonmagnetic coaxial transmission line. The radius of the inner conductor 
is a, and the inside radius of the outer conductor is b. From Ampere’s circuital law, it is easy 
to show that, fora<r <b, 


——— A/m (4.29) 


and therefore 


I 
B= UH = ue 
27r 


T (4.30) 


We cannot use ¢ = BA since B varies from inner to outer conductor. Instead, we integrate 
to find the flux crossing any radial plane extending from r = a to r = b and from, say, z = 0 
to z=. 


a) 
o=fBas= | | Hom ar dz = Ho m* Wb (4.31) 
0 Ja £. Tv 


20 a 


The flux links the current once, so N = 1. From Eqn. 4.4 the inductance for this length ¢ 


0) Mol, 0b 
. 
LD= += 5 In H (4.32) 


Suppose now that we try to use this expression to find the inductance of a segment of 
isolated straight conductor. As the radius of the outer conductor b —> oo, the corresponding 
inductance also becomes infinite. What this result tells us is that we never actually have a 
isolated straight conductor carrying a current without some return path. We must always 
consider the return path for current if we expect to get inductance values that have any 
relationship to reality. 


In the situation of a toroidal coil of N turns and a current J, as shown in Fig. 4.5, the 
magnetic flux density is 


NI 
gaia: a (4.33) 
2rr 
For a toroid cross section that is rectangular, as shown in Fig. 4.6, the integration is 


straightforward. 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Inductors and Transformers 4-11 


b PT ppptg NI NIT, b 
erro dr dz = PRO y= (4.34) 


r=<aJz=-0 277 20 a 


o= 


We then multiply the flux by N to get the total flux linkages, and divide by I to get the 
inductance. For the rectangular cross section case, this is 


N°T 6b 
| eg sca aa Fr H (4.35) 
27 a 


I 
~_— 
= \a 
_ 

I 


Figure 4.5: Toroidal Coil 


qi 


Figure 4.6: Toroidal Coil Cross Section 


4.5.1 Tesla Coil Inductance (Wheeler) 


Many times we just want to make a quick estimate of the inductance of some simple structure 
without making extensive calculations. Many approximate formulas have been developed in 
the days before hand calculators and computers, of which one will be given here. Watch out 
for the fact that the above formulas are given in the standard SI units, but the following 
formula is in the English system. 


The low-frequency inductance of a single-layer solenoid is approximately [8, p. 55]. 


r2N2 


ig a ay 4.36 
~~ Or+ 100 © eo) 
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where r is the radius of the coil and ¢ is its length in inches. This formula is accurate to 
within one percent for @ > 0.8r, that is, if the coil is not too short. It is known in the Tesla 
coil community as the Wheeler formula. The structure of a single-layer solenoid is almost 
universally used for Tesla coils, so this formula is very important. In normal conditions (no 
other coils and no significant amounts of ferromagnetic materials nearby) it is quite adequate 
for calculating resonant frequency. 


Example 


What is the approximate inductance of an air-cored solenoid with r = 8 inches, @ = 30 inches, and 
N = 175 turns? 


(8)2(175)2 1960000 
a IE en = 5270 pH 
9(8) + 10(30) 372 . 


If one needs the inductance for other geometries, Terman [8] has a number of expressions. A 
somewhat more recent paper by Fawzi and Burke [2] gives formulas for calculating the self and mutual 
inductances of circular coils in a form suitable for computer calculation. 


4.6 High Frequency Transformers 


Conventional low frequency transformers consist of coils of wire wound around steel lamina- 
tions. As mentioned earlier, the losses, especially the eddy current losses, become excessive 
at frequencies above a few kHz with this technology. Transformers built for operation at 
frequencies above a kHz or so are built around ferrite cores or air cores. Ferrite cores yield 
very compact, efficient transformers. Saturation limits operation to moderate power levels, 
however. If extremely large currents or powers are involved, then the air core transformer 
may be the logical choice. We shall discuss both types. 


The symbolic construction of a transformer is shown in Fig. 4.7. A voltage v,; produces 
a current 7; which in turn produces a flux ¢;. Part of ¢, links the secondary winding and 
produces a voltage vg by Faraday’s Law. If some load is connected, a current 72 will flow. 
This current produces a flux which opposes the original ¢; according to Lenz’s Law. This 
reduces the voltage induced in the primary so that more primary current will flow for a given 
source voltage. A complete description of transformer action in terms of field quantities has 
been developed only rather recently [4, 6, 9, 10, 16]. 


4.7 The Ideal Transformer 


It will be convenient to describe the actual transformer in terms of an ideal transformer. This is 
a transformer with no copper losses, no hysteresis or eddy current losses, and perfect magnetic 
coupling between primary and secondary. For such a device, the relationships between input 
and output voltages and currents are 
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Figure 4.7: A Two Winding Transformer 


uy NM 
—=— 4.37 
oe (4.37) 
i; No 
-=— 4.38 
a ie (4.38) 
The relationship between the input and output apparent power is 
N. 
Vy4t}1 = wa ey = 212 (4.39) 
The input impedance is 
N,/N: N? 
Fe IN 259 ae (4.40) 
al tgNo/N1 Ng 


The ideal transformer thus changes the level of voltage, current, and impedance between 
primary and secondary. 


4.8 The Actual Transformer 


A complete circuit model of the actual transformer is shown in Fig. 4.8. 


In Fig. 4.8, R, and Ro are the resistances of the primary and secondary windings, D; and 
Ly are the leakage inductances, R,, is an equivalent resistance representing the hysteresis and 
eddy current losses, and L,,, is the magnetizing inductance. The circuit is usually simplified 
by eliminating the ideal transformer and replacing all the impedances, voltages, and currents 
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ideal 


a Ty Ry Ro Lo 
— 
| 


Figure 4.8: Transformer Model 


on the secondary side with their equivalent values as seen by the primary. If we define the 
turns ratio @ as 


Perce (4.41) 


the simplified circuit is as shown in Fig. 4.9. 


i j(X1 + a?X2) Ri +a7Ry In/a 


Figure 4.9: Actual Transformer Referred to Primary 


We have also shifted to the phasor notation in Fig. 4.9, replacing inductances by their 
equivalent reactances (X1 = wl}, etc.), and the instantaneous voltages and currents by the 
phasor quantities. The notation can be shortened even more by defining an equivalent resis- 
tance and reactance 


Neg =R,+ aR (4.42) 
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Xeq = X1 +07 Xe (4.43) 


The load current is given by 


L 
== en (4.44) 
a Req + jXeq ta ZL 


The copper losses are given by 


Popper = (Ip/a)?Req (4.45) 
and the core losses are given by 
V2 
Poore = z (4.46) 


Other conducting materials can be used, such as aluminum, but it is tradition to refer to 
these series losses as copper losses regardless of the conducting material. It is good practice 
to keep the copper losses and core losses within a factor of two of each other, at least on large 
power transformers. The two losses tend to work against each other in a design. The core 
losses are reduced by reducing the maximum magnetic flux density in the core, which requires 
either a larger core cross-sectional area or more turns on each winding. Either approach 
requires more wire, which increases the copper loss. There may be instances where the core 
losses in a ferrite core used at high frequencies are much higher than the copper losses in 
reasonably sized wire. In such cases, one should go ahead with proper sizes rather than try 
to reduce the wire size and increase the copper losses to attempt to maintain some arbitrary 
parity. 


The air core transformer has no core losses, of course. R,,, can be removed from Fig. 4.9 
in such cases. 


4.9 Transformer Design 


Now we are ready to design a simple transformer. We want to select wire sizes, core size, core 
material, and number of turns on primary and secondary so that the transformer will meet 
the requirements without overheating, but without being so large that it is more expensive 
than necessary. 


The most important rating is the required voltage of operation. This determines core size 
and material, and the number of turns. The wire size is then selected to handle the transformer 
current rating without excessive copper losses. As we have seen before, the voltage is related 
to the flux by Faraday’s Law. 
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vg = No— (4.47) 
To a good approximation, the input and output voltages are sinusoidal. That is, 
vg = V2V2 cos wt (4.48) 


where V2 is the rms output voltage. 


The flux is then found by integration. 


1 2V: 
go= No [vvat = vale sin wt (4.49) 


The flux density B is given by B = $/A, where A is the cross- sectional area of the core. 
The maximum flux density, Bmaz, is found when sinwt = 1. 


V2V2 


Bmax = wN>A 


(4.50) 


Bmax is determined either by published data for a particular type of magnetic material, or 
by measurement. It is typically in the range of 1 T for low frequency laminated transformer 
steel, and in the range of 0.1 to 0.3 T for ferrite cores. The most efficient use of the magnetic 
material occurs when Baz is slightly above the knee of the magnetization curve. As the 
material saturates, the permeability 4: = B/H starts to decrease. The inductance is directly 
proportional to permeability, so the inductance starts to decrease also. But the inductance is 
defined as L = N@/i. The flux ¢ is proportional to the sinusoidal voltage so it does not satu- 
rate. Therefore, as the permeability and the inductance decrease, the current 7 must increase. 
The peak of the magnetizing current increases rapidly above the knee of the magnetization 
curve, and can exceed the peak of the rated current if the transformer voltage is increased too 
far. The magnetizing current becomes very nonsinusoidal, with a high harmonic content, at 
higher voltages. As a rough guideline, the peak of the magnetizing current should not exceed 
perhaps 10 % of the peak of the rated current. That is, if the rated current were 5 A rms, 
with a peak of /2(5) = 7.07 A, then the peak of the magnetizing current should not exceed 
about 0.7 A. 


Once we know Bygz, the minimum number of turns can be determined from the above 
equation. 


Example. 
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A Ferroxcube 204XT250-3F3 ferrite core is to be used for a transformer at 50 kHz. The rms input 
and output voltages are to be 10 V. Determine the proper number of turns on each winding. 


We first examine a published hysteresis curve for this material, which indicates that saturation is 
acceptable up to about 0.3 T, at least for low frequency operation. We then check a chart of core 
loss versus flux density which shows a recommended operating range of 100 to 300 mW/cm? for this 
material. If the frequency is above about 25 kHz, then B,,q, must be reduced to maintain the total 
heating in this range. At 50 kHz, Braz = 0.2 T causes a heating of slightly under 200 mW/cm?, which 
is deemed acceptable. From the published mechanical data, the area A is 0.148 cm?. The minimum 
number of turns is then 


v2(10) 


SEN A AY oe ts 
27 (50 x 103)(0.2)(0.148 x 10-4) urns 


No min = 


We would normally round up to the next higher integer, 16 turns. 


Suppose we were interested in using the same ferrite core for the same 10 V transformer at 60 Hz. 
Assuming Baz = 0.3 T, the minimum number of turns is 


v2(10) 


N. Mineo? BAN ACONT AT AO, A AL 
a. 2n(60)(0.3)(0.148 x 10-4) 


= 8450 turns 


Any attempt to wind this many turns through a toroid opening of only 0.312 inches inside diameter 
would be frustrating at best. This points out the fact that low frequency transformers must be relatively 
large. 


The power rating of a transformer can be determined from Ampere’s circuital law and 
Faraday’s law, which state, for the sinusoidal case, 


Hl=Ni (4.52) 


d 
v= no = Ndmagw Cos wt (4.53) 


Converting these equations to rms values yields the apparent power 


H 
S=VI= WNBA) = wBH Al = wBH(Vol) (4.54) 


where B and H are rms values and Vol is the volume of the magnetic material. We see 
that the apparent power is directly proportional to the frequency and to the volume of the 
transformer. This helps explain why aircraft use 400 Hz rather than 60 Hz. The transformer 
volume and mass are reduced by the same ratio, thus increasing the aircraft payload. 


Example. 


The Ferroxcube core of the previous example has a volume of 0.462 cm? and a relative permeability 
of 1800. What is the power rating at 50 kHz and a By, = 0.2 T? 
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= 1.283 W 


Boot Bas __ 2n(50,000)(0.2)2(0.462 x 10-8) 
(Fs) (Vol) = ay a800) (a x 10-7) 
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LUMPED RLC MODEL 


We are now at the point of looking at models for the Extra coil. I have already discussed 
the conflict between the Corum brothers and most of the Tesla coil community regarding 
distributed versus lumped models. My background is similar to that of James Corum, so I 
started this research with the assumption that the distributed approach would be superior. I 
spent considerable time learning some computer codes the Corums wrote, and even supplied 
them with codes to calculate inductance. But I was unable to predict all the interesting fea- 
tures with the distributed approach, and it appeared that a massive effort would be required 
if one wanted to do so. So I turned to the lumped RLC model. The inductance and capaci- 
tance values seem to be reasonably well understood, as detailed in the last two chapters, but 
resistance is a big challenge. I do not consider the following to be a definitive work on Tesla 
coil resistance, but rather a starting point for discussion. 


The first circuit we want to model is shown in Fig. 5.1. The IGBT driver (discussed in 
the next chapter) is driving a vertically mounted coil with a toroid on top. The driver and 
coil are physically separated for reasons of health and safety. The wires between the two have 
some inductance L, and capacitance C,. The obvious choice for these conductors would be 
a coaxial cable of appropriate voltage and current rating. However, the capacitance C; must 
be charged and discharged by the leading edge of the applied voltage pulse each half cycle. 
The associated current shows up as switching noise in the current sense resistors. One can 
reduce C, by about half by using an open wire transmission line rather than coax, which is 
what I did. Open wire lines are often operated as balanced lines (neither wire grounded). In 
this case, the open line is unbalanced (one wire ‘hot’ and one wire grounded). 


Damping of the high frequency components is improved by greater inductance and greater 
resistance in the line between driver and coil. Both inductance and resistance increase with a 
wire of smaller diameter. I tried both 14 ga and 8 ga wire, and chose the 14 ga PVC coated 
solid copper wire. 


The hot conductor must be protected from flashover from the top of the coil. Total 
protection is probably like total protection from lightning or tornadoes, probably not possible 
and certainly expensive. What I did will provide at least some protection. I took scrap 
sections of two inch thick styrofoam about one foot wide and cut a small groove down the 
center. The groove was oriented up and the 14 ga wire was laid in the groove. This would 
help prevent the 14 ga wire from arcing to the floor. Another section of styrofoam was laid on 
top the first. Then a length of coax was laid on top of the sandwich, directly over the 14 ga 
wire. The coax was connected to earth ground. A flashover to ground would normally hit the 
coax. A flashover to the hot conductor would have to miss the coax, go through two inches 
of styrofoam and then a layer of PVC to get to the 14 ga wire. This gets the probability of a 
damaging strike down to an acceptable value. 


The driver applies a square wave of voltage between the left end of L; and the ground 
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Driver 


Figure 5.1: Drive for Tesla Coil 


plane of the Tesla coil. This square wave is composed of an infinite series of cosinusoids, the 
fundamental and all odd harmonics. The peak value of the fundamental is 4/7 times the dc 
supply voltage. That is, if Vz. = 500 V, then the instantaneous voltage of the fundamental is 
uv; = 637coswt. The first three terms of the Fourier series are 


4 1 1 
f(t) = —[cos wt — 3 cos 3wt + 5 cos ut 4+++5] (5.1) 
us 


The harmonics of the exciting wave will drive higher order resonances of the Tesla coil, if 
these resonances are harmonically related. It appears, at least for the coils I have built, that 
any higher order resonances are not exact multiples of the fundamental frequency. That is, I 
apply a square wave of voltage to the feed point, and observe a current that looks sinusoidal 
at the fundamental frequency. The lumped RLC model automatically excludes higher order 
resonances, so if they are of significance, we must use a distributed model to describe them. 
My feeling at the time of this writing is that higher order resonances are not a problem, at 
least not enough of a problem to exclude the use of the lumped model. 


5.1 The RLC Model 


The simplest model that can be proposed for the Extra coil is the series RLC circuit shown in 
Fig. 5.2. The inductance Ly is given by Wheeler’s formula. It can also be calculated from first 
principles [2] using elliptic integrals and a computer. The two methods typically agree within 
one percent. It can be measured to the same accuracy with a hand-held inductance meter. 
The only thing that changes its value is significant amounts of ferromagnetic materials inside 
or near the coil. These are avoided because of high losses. The inductance does not change 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 5—Lumped RLC Model 5-3 


with frequency or when sparks occur. If everything about a Tesla coil were as well behaved 
as the inductance, modeling would be easy. 


Vio Vo2 Voase 


Figure 5.2: Tesla Coil with Series Resonant LC Circuit 


The capacitance C2 is the modified sum of the capacitance of the coil itself, as given by 
Medhurst, and the isotropic capacitance of the toroid. The toroid and the coil shield each 
other, so the effective capacitance is always less than the sum of the Medhurst capacitance 
and the toroid capacitance. There should be a simple empirical formula to calculate Co, given 
the Medhurst and toroid capacitances, but I do not know what it is. Cy can be calculated 
numerically using Gauss’s Law. Terry Fritz has written a program to do this. If one is careful 
about measuring and entering all the dimensions and the locations of grounded surfaces, one 
should get a value for C well within 5% of the correct value. Of course, Cy increases when 
the spark occurs due to the size of the toroid becoming effectively larger. This lowers the 
frequency of operation while the spark is present. 


There is really no way to directly measure Cj. C2 is in the range of tens of pF, the same 
range as the leads of a capacitance meter. The presence of a meter and operator will change 
the capacitance of the coil. One can measure Ly and the resonant frequency and then calculate 
Cy from these measured values. 


If only a crude approximation is desired, one can use the Medhurst capacitance plus 75% 
of the isotropic toroid capacitance for Cp. Since resonant frequency is related to the square 
root of capacitance, a 20% error in capacitance results in only a 10% error in frequency, which 
is sometimes close enough. 


Now we come to the third element of the model, Rec. This is the input impedance of the 
coil at resonance. We have seen that we can readily get a value for Lo to within 1%, and a 
value for Cy to within 5% with a little more work. Or we can measure frequency and calculate 
Cy to within one or two percent. But Rrc is another story, as I have mentioned. For now, 
let us assume we have an appropriate value for Rrc and explore the features of the model. 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 5—Lumped RLC Model 5-4 


RLC Circuit With Sinusoidal Source 


We have seen that a square wave can be decomposed into a sine wave (or cosine wave) of 
fundamental frequency plus an infinite series of harmonics. If the harmonics do not have 
much effect upon operation, we can safely model the square wave as a sine wave of the same 
frequency, at least for some purposes. For brevity, we will call the three circuit elements just 
R, L, and C, and redraw the circuit slightly as shown in Fig. 5.3. If only a single frequency 
sine wave is present, we can use phasor analysis. The input phasor voltage V, drives a current 
I, through the circuit. By Kirchhoff’s Law, 


Vi=VL_t+Vct+VR (5.2) 


Figure 5.3: Series Resonant LC Circuit 


Using Ohm’s Law in phasor form gives 


1 
Vi =0hZ=hjwlL+ — +R) (5.3) 
jw 


The voltage across the inductor leads the current through the inductor by 90°, while the 


voltage across the capacitor lags the current by 90°. The two voltages tend to cancel each 
other, and actually do cancel at the resonant frequency wo, as defined by 


1 


b= 5.4 
0) uC (5.4) 
We then solve for the resonant frequency in radians per second as 
1 
Wo = = (5.5) 


The resonant frequency in Hertz is 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 5—Lumped RLC Model 5-5 


1 
= 5.6 
fo= sie 68) 
The series impedance, expressed as a phasor, is 
Z = 2lo=R+ 5 (wh =) (5.7) 
7 _ oN wC , 
where 
1 2 
Z=4/R?4 L-— 5.8 
(wr-—) (5.8) 
and 
L-(1 
@ = arctan eee) (5.9) 


At resonance, Z = FR, a real number, so the input current I, is in phase with the input 
voltage V1. This means that the voltage Ve will lag the input voltage Vi by 90°. We will 
say more about this in a moment. 


The bandwidth of the series circuit is defined as the range of frequencies in which the 
amplitude of the current is equal to or greater than 1/\/2 times the maximum amplitude. 
This current produces a heating effect of half the maximum value, so the frequencies at the 
bandwidth limits are called the half-power frequencies w, and wy. With a little algebra, the 
bandwidth is shown to be 


R 
B=wo- wy =i (5.10) 
The quality factor Q is defined as the ratio of the resonant frequency to the bandwidth, 
so for the series circuit we have 


_ wo wol | 1 1 L 


@ BR wR RVC 


(5.11) 


The amplitude of the voltage across either the inductor or the capacitor at the resonant 
frequency wo is Q times the amplitude of the source voltage. 


The above is standard Circuit Theory I. We see that to get a large voltage in our Tesla coil 
that it needs to be high Q, having large L and/or small R. With this model, a space wound 
coil of 14 gauge wire (called 14S) that I use for testing would have a Q of about 500. The 
inductance is about 17 mH, the resonant frequency is 160 kHz, and the input impedance Rro 
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is about 25 Q. If I apply a sinusoidal source voltage with an rms value of 1000 V, the voltage 
across the inductance should rise to 500,000 V. The voltage does not actually reach this value 
because a spark occurs first, but if the toroid were large enough and smooth enough, it should 
reach this voltage. 


The RLC model predicts such a voltage, but how do we check the model for validity? Does 
the Tesla coil toroid actually reach 500 kV? If we try to measure the toroid voltage directly, 
we have to attach some sort of probe, which must have some capacitance and resistance. The 
probe capacitance lowers the resonant frequency. At least part of Rc is frequency dependent. 
The probe may have frequency dependence as well. We can measure the toroid voltage with 
the new Q and new Reco but we cannot be positive about the toroid voltage when the probe 
is removed. 


I have spent considerable effort to measure the toroid voltage with a battery operated, 
fiber optic coupled, electric field probe. Results are very repeatable and show some interesting 
features, but there is always the question about absolute calibration. To calibrate the electric 
field probe, I have to connect a voltmeter to the toroid at some point. In addition to changing 
the resonant frequency and Rc, the voltmeter probe will also change the electric field pattern. 


So we agree that we cannot measure the toroid voltage to within 1%. But what is a 
reasonable estimate of the accuracy of measurement? I would estimate that the uncertainty 
of the absolute value of toroid voltage is about 10%, almost certainly no worse than 20%. I 
used to say that such an accuracy was close enough for government work, but such a statement 
would be politically incorrect today. Actually, I quit using the phrase after my uncle, a long 
time employee of the Corps of Engineers, became angry when I said it. 


Sparks are a statistical phenomena, anyhow. Even when the toroid voltage is accurately 
known, say in a Marx generator, there is uncertainty about whether the spark will occur and 
about how long it will be. One voltage may produce a spark one time out of ten, and a higher 
voltage may see a spark nine times out of ten. These two voltages may easily differ by more 
than 10%. So if I know the toroid voltage to within 10%, it should be quite adequate. 


Another way of estimating the toroid voltage is to observe the spark length and multiply 
by some constant. Two polished brass spheres separated by one cm will break down in dry 
air at 30 kV. So, the assumption goes, if we have a 10 cm spark from toroid to air, the toroid 
voltage must have been 300 kV. This is one of those great ideas that just do not work in 
practice. One of the main problems is that spark length is a nonlinear function of voltage. 
Once a spark is started, it does not take as much voltage to keep it going. 


An example of the error involved is given by the following: I was applying a square wave of 
voltage with an rms value of about 1200 V to the Tesla coil. With a Q of 500, the maximum 
possible toroid voltage would be 600 kV. Sparks were occurring at lower voltages, probably 
not over 500 kV. I observed a spark 143 cm long to a ground rod. Using 30 kv/cm would 
predict a toroid voltage of 4350 kV, about a factor of nine too high. This technique certainly 
gives one the satisfaction of claiming big voltages, but doesn’t have much relation to reality. 
One could scale the constant down to say 3 kV/cm and be closer to the true value, but the 
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nonlinearity of the spark keeps this from being very good either. Measuring the electric field 
with a fiber optic coupled probe should get much better results for the toroid voltage than 
the technique of multiplying by a constant. 


Anyhow, I believe the RLC model does a reasonable job of predicting the toroid voltage 
before breakout, to within perhaps 10 to 20% of the ‘correct’ value. The streamer adds 
resistance to the system, which lowers Q and toroid voltage. The model and experimental 
results seem to be in reasonable agreement on these features. 


RLC Circuit With Square Wave Source 


The sinusoidal steady state analysis of the previous section does not tell us how the input 
current builds up when voltage is applied. For that we need to go back to the time domain. 
We will use the actual square wave input voltage and ignore the resistance for the first few 
cycles. Until the current builds up, the voltage across the resistor will be much smaller than 
the voltage across either the inductor or capacitor, so this is not a bad approximation. 


The input voltage v; is equal to Vg the first half cycle and —Vg the second half cycle. The 
current everywhere in the circuit is 7 and the capacitor voltage is v,. The initial current is Ip 
and the initial capacitor voltage is V9. Kirchhoff’s voltage law is 


di 
L— e— 12 
dt +v Va (5 ) 


The component equation for the capacitor is 


Ve 
2 dt 


=i (5.13) 


The solution of this set of equations for t > to is as follows: 


Va — Vi 
i(t) = Ip cos wo(t — to) + ———® sin w(t — to) (5.14) 
0 
and 
v¢(t) =Va- (Va = Vio) cos wo(t = to) + Zolo sin wo(t — to) (5.15) 
where the angular resonance frequency wo is 
1 
Wo = 27 fo ar — (5.16) 


VLC 


and the characteristic impedance (sometimes called the surge impedance) Zp is 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 5—Lumped RLC Model 5-8 


262 12 Q (5.17) 


We assume we are operating at resonance, wo of the applied square wave = wo of the RLC 
circuit. The time axis is selected so tj) = 0 and t = 0 at the start of the square wave. The 
initial current through LF and the initial voltage on C are both zero. The equations for 7 and 
V.. become 


and 
ue(t) = Va — Va cos wot 


The current starts and ends at zero, with a peak value of Vg/Zp. The initial current for 
the next half cycle is zero. Capacitor voltage starts at zero but builds to a value of 2Vq, which 
is the initial value for the next half cycle. Waveforms are shown in Fig 5.4. 


For the second half cycle, the applied voltage is —Vg. The starting time to is half the 
period of the square wave. The current is 


EVE = 2Va 


i(t) Fi 


sin Wo (t a to) 


and the capacitor voltage is 
u-(t) = —Va — (—Va — 2Va) cos wo(t — to) 


The current starts and ends at zero, but this half cycle it reaches a peak value of —3Vq/Zo. 
The factor t — to starts at zero at the beginning of the second half cycle so the cos function 
starts at +1 and goes to —1. The capacitor voltage starts at 2Vz, as it must, and ends at 
—AV 4, which is the initial capacitor voltage for the third half cycle. 


One can continue on with this solution process as long as desired. The peaks for current 
are Vy/Zo, —3Va/Zo, +5Va/Zo0, —TVa/Zo, etc., occurring at the midpoints of each half cycle. 
The peaks for capacitor voltage are 2Vqz, —4Vq, 6Vq, —8Va, etc., occurring at the ends of each 
half cycle. The capacitor voltage lags the current waveform by 90°, the standard relationship 
between voltage across and current through a capacitor. 


This suggests one way of measuring Zp without measuring the resonant frequency. We 
apply a half cycle of voltage (actually somewhere between one third and two thirds of a cycle), 
and measure the peak current. The problem is that the peak current during the first half cycle 
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Va 


Vi 


Figure 5.4: First Two Cycles of Current and Capacitor Voltage in RLC Circuit 


may be much smaller than the charging current of the transmission line to the Tesla coil, which 
would mask the desired effect. Since the toroid voltage is only double the input voltage after 
the first half cycle, and electric fields are small, we might be able to put a fiber optic coupled 
pulse generator right at the input terminals of the Tesla coil, and minimize the effects of stray 
capacitance. On the other hand we could very well discover that the distributed capacitance 
effects will always mask the desired measurement. 


As mentioned earlier, these waves will not increase forever, but will asymptotically ap- 
proach values controlled by the circuit resistance, (if a spark does not occur first). At steady 
state we will see the relationships found in the previous section. 


Phase Shift, Input to Toroid 


One interesting issue is the amount of phase shift from bottom to top of a Tesla coil. There 
are strongly stated opinions among Tesla coilers as to what this shift is, or should be. Some 
believe the phase shift is zero. I believe the toroid voltage (with respect to ground) lags the 
input voltage (also measured with respect to ground) by 90°. On one occasion I calibrated 
the various phase shifts in my instruments (described in the following) and arrived at a phase 
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shift of 94° from base to toroid. This fits the RLC model well within experimental error. 


There are several different ways of thinking about the phase shift, some quite detailed and 
prone to error. Any of us can get sloppy on notation and make such an error. So I will try to 
be precise. 


Earlier in the chapter I showed that current is in phase with the input voltage for the 
lumped RLC model at resonance. The model predicts that the toroid voltage will lag the 
input voltage by 90°. One concern here is that part of the total capacitance is internal to 
the coil, which confuses the computation of the toroid voltage magnitude. However, if the 
current entering the toroid is in phase with the current entering the base of the Tesla coil, as 
the RLC model requires, then this current flows through a pure capacitance from the toroid 
to ground. So the phase of the toroid voltage would be 90° behind the base voltage, even if 
the magnitude is in error. 


We get the same result if we think of a Tesla coil as a quarter-wave antenna above a 
ground plane. If the antenna impedance matches the transmission line impedance (quite 
possible) then only a traveling wave exists on the line. The output of a quarter-wave section 
of transmission line always lags the input by 90°, and likewise the top of a quarter-wave 
vertical antenna with respect to the base. The Tesla coil is not the same as a quarter-wave 
whip antenna because of the mutual inductance and mutual capacitance, but I think this is 
not adequate to change the phase shift from 90° to zero. 


But how does one actually measure the phase shift? A well done measurement can be 
more convincing that a basket of theoretical arguments. As discussed earlier, the presence of 
a probe presents some serious problems with loading the coil and affecting both magnitude 
and phase. It appears that the next best thing to having a physical probe attached to the 
toroid is to measure the electric field. With adequate calibration, this should answer most of 
our questions. 


Suppose we have our Tesla coil inside a grounded metal building. Some flux lines go up 
and hit the roof. Some go sideways to the walls. Some go down to the earth. Each flux line 
(or tube) can be considered to be between a portion of a capacitor plate at Tesla coil voltage 
and a similar portion at ground potential. That is, the capacitance of the Tesla coil can be 
thought of as hundreds of sub-capacitors going in all directions, and connected in parallel. 
If we take one of these sub-capacitors and put an electrode in it, (not at ground or toroid 
potential), we have just split the capacitor into two series capacitors. Voltages divide across 
capacitors inversely as the capacitances. The voltage between the electrode and the adjacent 
surface can be readily measured by a standard electronics circuit if its capacitance is much 
larger than the value of the sub-capacitor between toroid and ground. The concept for an 
electrode near ground is shown in Fig. 5.5. 


The electrode can also be mounted on the toroid, with the transmitter inside the toroid. 
If the electric field is measured at the toroid, it is essential that the electronics be battery 
powered, and the signal transmitted away from the coil by an optical fiber. It is not as essential 
if the voltage is measured near ground potential, but is still a good idea. The alternative, 
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voltage optical 
follower, 
transmitter 


portion 


of toroid 


Figure 5.5: Concept of Electric Field Measurement 


a coaxial cable to the electronics, will pick up some of the signal from the Tesla coil on its 
ground braid. RF on a ground lead can cause erroneous results. 


Actually there are four distinct locations where the electric field might be sensed. 


1. Next to toroid 


i) 


. At bottom of coil 
3. At grounded surface above toroid 


4. At toroid height by grounded wall 


Each location has its own advantages and disadvantages. I started experimentation with 
fiber optic systems by putting the transmitter in the toroid. This has the advantage of being in 
a high electric field. The electrode can be physically small and/or located close to the toroid. 
One disadvantage is that a spark at that location could destroy your electronics. Another 
disadvantage is that the smaller toroids are difficult to modify and do not have enough room 
for the electronics. A half-spun toroid worked quite well, but the classic toroids did not. 


My next attempt was to place the electronics at the bottom of the coil. I put a metal disk 
at the top of the coil form, about two thirds of the diameter of the secondary, and electrically 
connected to the toroid. I then put a polyethylene bracket across the bottom of the coil 
form with a smaller disk facing up. A wire from this disk ran through the bracket and into 
an aluminum box holding the electronics, which was mounted to the bottom of the bracket. 
This worked quite well for testing at lower powers. The electric field was strong and was not 
affected by movement of people or things external to the coil. One concern was that the sensor 
would ‘see’ the sides of the coil. The voltage would be lower on the sides but the distance 
would be smaller, so the influence would be similar. This posed the most problem for looking 
at the toroid voltage during discharge. The sensor would show a voltage decay lasting for five 
to ten cycles at discharge. The toroid voltage may have dropped much faster than that, but 
one could not be sure because the coil sides took a longer time to discharge. 


The main problem, however, was the occasional spark from top electrode to the grounded 
aluminum box through the inside of the coil form. This spark was about twice as long as the 
sparks to air at that time, hence was a bit of a surprise. On at least one occasion, that spark 
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to ground upset something in my driver that caused my IGBTs to give up. I made a number 
of changes to my system after that, including removing the grounded aluminum box. I also 
added some polyethylene disks to block the inside of the coil. Now a toroid to ground strike 
has to go through several layers of polyethylene plus two or three feet further to get to a real 
ground, if it tries to travel on the inside of the coil. 


Philosophically, the best place seems to be above the toroid, at ceiling level. Directly above 
the coil, the toroid appears as a large disk. The influence of the toroid would be maximum 
while the influence of the coil sides would be minimum. It is also the most inconvenient 
location regarding the changing of batteries, switching ranges, and the like. For this to work, 
the ceiling must be conductive, and immune to corona. That is, a large area of the ceiling 
must be covered with a good conductor (like copper sheeting), or a large grounded toroid 
must be used, with the sensor at its center. I will consider doing this if I ever build another 
lab. For now I will place the sensor at about toroid height, and located next to the grounded 
metal wall of the building. 


It has been suggested [1] that the grounded toroid placed directly above the coil could be 
on a rope and pulley to make it easier to change batteries. In my case I already have a rope 
and pulley in place to raise and lower a large and awkward toroid into position on top the 
coil. If the grounded toroid was not too large and heavy, it could be clipped on the rope and 
raised into position after the high voltage toroid was lowered into place. I may try that. 


It was also suggested [1] that the system could be calibrated by charging the high voltage 
toroid to some known voltage and discharging it directly to earth through a low impedance 
path. The coil would be physically in place but disconnected at top or bottom or both. The 
capacitance of the toroid would resonate with the inductance of the shorting path at some 
frequency well into the MHz range. If the electronics were still linear at this high frequency, 
this might work very nicely. 


In summary, the simple RLC model appears to do reasonably well in predicting 


1. Resonant frequency 
2. Toroid voltage magnitude, given the resistance 


3. Toroid voltage phase 
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RESISTANCE OF COIL 


The resistance Rpg in the RLC model is an effective or equivalent resistance which rep- 
resents all the losses in the Tesla coil. It includes 


1. Ohmic or copper losses 

2. Dielectric losses, coil form and conductor insulation 
3. Eddy current losses in toroid, strike ring, and soil 
4. Radiation losses 


5. Losses in the spark 


It is surprising how difficult it is to calculate these various losses. Making meaningful 
measurements can also be challenging. If we operate at low input voltage so we are below 
spark breakout, then we can ignore the last term for the moment. I will ramble through some 
considerations for the other losses. 


6.1 Temperature Effects 


Almost all coils are wound with copper wire. It is moderately priced and widely available. 
There might be an occasional aluminum coil, usually from some ‘bargain’ at a surplus auction. 
Aluminum has higher resistivity than copper, so to get a given resistance the wire must be 
physically larger. We saw earlier that to get a high toroid voltage we needed a coil with large 
L and/or small R. If we use larger wire to keep the resistance the same, the coil must be 
physically larger and the inductance will decrease. We would expect therefore that aluminum 
coils would always be inferior to copper coils. 


Example 


You are given a choice between two spools of magnet wire, each 1000 feet in length. The copper is 
24 gauge, with nominal resistance 25.67 Q, while the aluminum is 22 gauge, with nominal resistance 
26.46 Q. You have a piece of 5 inch diameter PVC pipe and are considering using the entire 1000 ft to 
wind a coil. You wish to compare the inductances of the two prospective coils. 


You assume a build (thickness of dielectric) of 1.65 mils. The diameter of the 22 gauge wire is then 
23.35 + 2(1.65) = 28.65 mils and the diameter of the 24 gauge wire is 20.1 + 2(1.65) = 23.4 mils. The 
nominal number of turns would be 
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The winding length for the 22 gauge wire would be 


Ly = 764(28.65) = 21,890 mils = 21.89 inches 
and 17.19 inches for the 24 gauge wire. 


By Wheeler’s formula, the inductance for the 22 gauge coil would be 


2N? 2.5)?(764)? 
pe eS n= eid 
9r+10€,  9(2.5)+ 10(21.89) 
The 24 gauge coil has an inductance of 18,120 wH. For this particular example where r and N are 
held fixed, ¢,, is proportional to wire diameter and since it appears in the denominator, a smaller wire 
diameter results in a larger inductance. 


If we hold the winding length fixed, say at 17.88 inches, then the 1000 ft of aluminum wire will fit 
on a coil of 624 turns and 6.121 inches diameter. The inductance is now 17,670 wH, larger but still 
less than that of the copper coil. The aluminum coil has more resistance and less inductance than the 
copper coil, both undesirable. As a side issue, if we hold the winding length fixed, the coil gets to be 
relatively short and fat. Standard wisdom is that the ratio of length to diameter needs to be on the 
order of four or five, so this is another negative factor. You should buy the copper wire. 


Aluminum forms a coating of aluminum oxide when exposed to air. This oxide is not a 
conductor (copper oxide is) so there is always a problem in making electrically solid connec- 
tions. This is not an unsolvable problem since electric utilities use aluminum wire almost 
exclusively, for cost and weight reasons. The trained individual with the right tools might be 
able to get acceptable performance from an aluminum coil. The rest of us should stick with 
copper. 


Silver is slightly more conductive than copper, but the price is much higher. Copper losses 
are not considered that big a factor in coil performance, so there is little incentive to go to 
silver wire. 


The dc resistance of copper wire is determined by table look-up or by measurement with an 
ohmmeter. The tables typically give the resistance of 1000 ft of wire at a specified temperature, 
say 20°C or 68°F. The tables give the resistance to four significant places, but using more 
than three is somewhat of a joke because of the resistance variation with temperature. If the 
table resistance is Ry; at a temperature T,, then the resistance R2 at some other temperature 
T» is 


(6.1) 


where T; is the inferred absolute zero temperature, —234.5°C for copper. If the table resistance 
is given for 20°C, then the resistance of a copper wire at temperature T> is 


Tp 40345 
Rei 6.2 
2= fis + 934.5 Oi?) 
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The resistance has dropped to 90% of the tabulated value at —5.45°C (not impossible for 
an unheated shop in some parts of the world) and is at 110% of the tabulated value at 45.45°C. 
The wire could certainly rise to this temperature during operation even on a pleasant day. 


The inferred absolute zero temperature does not have much to do with absolute zero 
(—273°C). Pure iron has an inferred absolute zero of —162°C and Manganin has an inferred 
absolute zero of —167,000°C. The fact that the inferred absolute zero of copper is not dras- 
tically different from absolute zero makes it useful as a temperature sensor. If one is curious 
about how warm a Tesla coil got during operation, the best way of measuring the temperature 
is to measure the dc resistance and work backwards through the above equations. 


6.2 Skin Effect 


At dc, current flows uniformly across the entire cross section of the conductor. As frequency 
increases, however, a phenomenon called skin effect causes less of the total current to flow in 
the center of the wire. Having less conductor available causes the resistance to increase. 


An expression for skin depth can be derived as 


hoses (6.3) 
Va f uo 
where f is the frequency in Hz, pz is the permeability of the conductor (4m x 10~" for nonfer- 
romagnetic materials), and o is the conductivity. 


The skin depth for copper at 20°C is 
Pa (6.4) 


Most introductory electromagnetic theory books derive the expression for ac resistance as 


b 
Rac = Raex 6.5 
ac dc 26 ( ) 
where 0 is the radius of the wire and 6 is the skin depth, in consistent units. This equation is 
only valid for 6 < b. As might be expected, this excludes most Tesla coils, so we must find 
other expressions. 


A number of more advanced electromagnetic theory books derive an expression for the ac 
resistance. Typical is the treatment by Ramo, Whinnery, and Van Duzer [6]. They start with 
Maxwell’s Equations, and write a differential equation for the current density J,(r) inside an 
isolated straight cylindrical conductor centered on the z axis. 
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ddy. Vay 
74+ -—+T* J, =0 6.6 
dr? . r dr = 6:8) 

where 
2 . : _2 
TY = —jwpo = —j(2n fur) = ig (6.7) 
The solution for current density is 

J, = AJo(Tr) (6.8) 


where Jo(Tr) is the Bessel function of the first kind and zero order. 


The infinite series for this Bessel function is 


ore) (—1)” DENA OO (_1)M(T2)mp2m 
2) Sumer) 


Io(Tr) = S> — (6.9) 
m=0 ( m=0 2 (ml)? 
When we insert the above expression for T?, the series becomes 
ym (2\™ 2m 
oo (-1)™(-A)"(B)or 
= ) 
Note that 
ly apy Sg (6.11) 
2 1 
p2m r 2m 
With these substitutions, we can write 
2m 
oo (9) (5) 
Jo(Tr) = >° Pml)z (6.14) 


The series can now be split into its real and imaginary parts. 
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r\4 
Real|Jo(Tr)] = 1— Gs fe NOT 


a2(2N2 * 24412 
r\2 r\6 ryi0 
mae lo(tr)| = Gy 7 cm " ve ue 


These are related to the ber and bei functions as follows. 


Real|[Jo(T'r)| = ber(v2r /6) 


Imag[Jo(Tr)] = bei(v2r/6) 


6-5 


(6.15) 


(6.16) 


(6.17) 


(6.18) 


There are several ways to proceed to find the ac resistance. One is to use one of Maxwell’s 
curl equations to find the magnetic field inside the wire, a constant times the derivative of the 
electric field. This magnetic field is used to find the total current in the wire. The electric 
field times the wire length @ is the total voltage. The ac resistance is just the real part of the 


ratio of voltage to current. 


To my knowledge, every method of finding the ac resistance is tedious. I will proceed with 
a basic Circuit Theory I type approach. Consider the straight wire as formed of N concentric 


cylinders, as shown in Fig. 6.1. 


Figure 6.1: Wire Formed of Concentric Cylinders 


The wire radius is 6. The radius to each boundary between cylinders is designated by rz, 


where ry = b, r; = b/N, ro = 0, etc. The cross sectional area of cylindrical shell k is 
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Area = n(r? — rz_,) (6.19) 


The dc resistance of cylindrical shell & is 
pe 


m(rZ a 7-4) 


where p is the resistivity and @ is the length of the wire. The dc resistance of the entire wire 


Rack = (6.20) 


is 


(6.21) 


The average current density in a shell is approximately the current density at the midpoint 
of a shell. That is, 


Joe = AJo(Tr) for = (6.22) 


The current in shell k is 


Th = Jenn (rk — TR-1) (6.23) 


Note that J; is a phasor. Both the magnitude and the phase of the current and current 
density vary from the center to the surface of the wire. 


The total current is given by 


N 
T= S- I, (6.24) 
k=1 


where all the real parts of the various J, are added together, and all the imaginary parts, to 
form the phasor J. The power dissipated in shell k is 


Py = [ply Rack (6.25) 


where J;,J; indicates that the current has been multiplied by its complex conjugate to form 
|I;,|?. The total power dissipated in the wire is 


P= > Pr (6.26) 
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But by definition 


P=lIP Rez (6.27) 
Therefore, 
Pe ee 
Roc = a = i) 
me uP wey 
and 
Rac yl Pr 
— = a 6.29 
Rae Rae|L|? ( ) 


We need to do one more step to get this formula into computational form. The ratio of 
resistances is independent of a specific current flow so we can set A = 1 in Eq. 6.22. For 
shorthand we will write the real part of Jz, as REJ; and the imaginary part as IMJ;. After 
pages of algebra, all the terms like pé@ cancel out and we have the final result. 


Pe _ SRL REV? + 1M] (48) — (B=4)”) ian 


Me feseamea (a) —())] + etme ((9)'- CY) 


The algebra basically normalizes itself to use a cylinder of unit radius. The variable is 
skin depth, which appears in the argument of the Bessel function as r/6. If the wire radius is 
one skin depth, the ratio Rac/Rac = 1.020, so the de resistance can be used for any wire for 
which the radius is less than a skin depth. Other values are given in Table 6.1. 


Table 6.1: Rac/Rac for various values of b/6. 
b/d Rael faa 
1.020 
1.263 
1.763 
2.261 
2.743 
3.221 
3.693 
4.154 


CONDO KBR WN EH 
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This analytic technique of finding the ac resistance has been around for a long time and is 
described in many books. Before the advent of computers, it was difficult to use. One would 
use tables of Bessel functions which listed real and imaginary components separately, such as 
Jahnke and Emde [4]. The notation among different authors differed, making it difficult to 
compare results. Of course, computers have now made the task of calculating these infinite 
series much more manageable. Actually, if the wire radius is not more than about 8 skin 
depths, the first 12 terms of Eq. 6.14 are more than adequate. 


Because of the difficulties of using this analytic technique before computers, people devel- 
oped lookup tables and empirical models to find Ra-/Rac. For historical interest and for the 
benefit of any reader adverse to writing computer code, we will present the approximations 
developed by Terman [8], who has a detailed discussion of this topic. He defines R,, in terms 
of a parameter x, where 


2f 
oe 6.31 
plo”) tee 
for nonmagnetic materials. Here, d is the conductor diameter in meters, f is the frequency in 
Hz, and p is the resistivity in ohm meters. As x gets very small, due to either low frequency 
or small wire, the ac resistance approaches the dc resistance. Above about « = 3, Rac/ Rac 


varies essentially linearly with x according to the expression 


Rae 
Rae 


= 0.3535a + 0.264 (eS 3) (6.32) 
Terman gives the following tabular values of Ra-/Rac for x between 0 and 3. 


Table 6.2: Rac/Rac for various values of x. 
0 1.0000 1.5 1.026 
0.5 1.0003 1.6 1.033 
0.6 1.0007 1.7 1.042 
0.7 1.0012 1.8 1.052 
0.8 1.0021 1.9 1.064 


0.9 1.0034 2.0 1.078 
1.0 1.005 2.2. 1.111 
1.1 1.008 2.4 1.152 
1.2 1.011 2.6 1.201 
1.3 1.015 2.8 1.256 
1.5 1.020 3.0 1.318 


For copper, p = 1.724 x 10~° ohm meters. For those of us still using wire tables in English 
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units, where wire diameters are given in mils (1 mil = 0.001 inch), Terman [8] has reduced 
the expression for x to 


& =0.271dmV/Funz (6.33) 


where d,, is the wire diameter in mils and fyyyz is the frequency in MHz. 


For example, the 14 ga coil being examined here has a dc resistance of 3.99 Q at 20°C. 
For a large toroid, the resonant frequency is 160 kHz. The nominal diameter of 14 ga wire is 
64.08 mils. We calculate x as 


x = 0.271(64.08) V0.16 = 6.946 


The ac resistance of the wire in the coil (assuming the wire is uncoiled and is supported 
in one straight line) is then 


Rac = 3.99(0.3535(6.964) + 0.264) = 10.85 © 


We see that skin effect makes a significant difference in resistance, especially where larger 
wire sizes or higher frequencies are used. 


The two methods presented here (Terman and the one using Bessel functions) should yield 
very nearly the same results if we use the relationship 


b 
r= 1412- (6.34) 
For example, if b/d = 3, then x = 4.236 and 
Rae _ = 
Foe = 0-3535 (4.236) + 0.264 = 1.761 (6.35) 
‘dc 


which is reasonably close to the value of 1.763 given in Table 6.1. 


6.3. Proximity Effect 


The effect of adjacent turns in the coil causes the current density to be even more nonuniform 
than for the straight wire, which raises the resistance even more. This effect is called the 
proximity effect. Terman [8] has a curve for two straight parallel cylinders carrying current 
the same direction that shows an increase of about 33% for the wires touching physically (but 
not electrically), and an increase of about 10% for the case when the two wires are separated 
by a gap equal to the wire diameter. 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 6—Resistance of Coil 6-10 


Medhurst 


But our problem is not that of two parallel wires going to infinity, but that of hundreds of 
turns of wire in a finite space. Perhaps the first person to examine this problem experimentally 
was Medhurst [5]. We saw his results for the self-capacitance of a coil back in Chapter 2. The 
same paper gives tables for ¢yy, the increase in resistance over Ra, as a function of coil length 
over coil diameter and wire diameter over wire spacing. Since such tests are critical to our 
Tesla coil model, it seems appropriate to discuss his paper in some detail. 


Medhurst used between 30 and 50 turns of bare copper wire wound in grooves of a low 
loss dielectric former. The former material was Distrene, which he claimed to have a power 
factor of about 0.0003. The coil and former were dried carefully before testing. This means 
that dielectric losses in the coil form, winding insulation, and humidity should be negligible. 
He did not mention the possibility of eddy current losses, but I suspect these were negligible 
as well. Testing was done at low power, so corona or spark losses would have been zero. His 
test method did not separate ohmic (heating) losses from radiation losses. Standard wisdom 
is that radiation from small coils at relatively low frequencies is negligible, and I have seen no 
contrary evidence, so it will be assumed that his results apply to ohmic heating. 


He used two different wire sizes, 18 and 20 s.w.g. I assume these refer to the British 
Standard Wire Gauge, where 18 s.w.g. has a diameter of 48 mils (1.219 mm) and 20 s.w.g. 
has a diameter of 36 mils (0.9144 mm). For the tight wound case, he used double-silk-covered 
(dsc) or single-silk-covered conductors. Such coatings are not readily available today, so it 
would be hard to replicate his findings. 


A confusing aspect of his paper is that the symbol H is used in four different ways. It is 
used for magnetic field and for the unit of inductance, both of which are immediately obvious 
and not a problem. It is used for a function of coil length over coil diameter as we saw in 
Chapter 2, where Cyy = HD for the self-capacitance (or Medhurst capacitance) of a coil. 
Then the symbol is also used for a function of wire diameter and skin depth in determining 
the ac resistance. Some other symbol would have been more appropriate. 


The coils Medhurst tested are not ‘typical’ Tesla coils. Tesla coils usually have far more 
than his 30 - 50 turns. He tested at frequencies of around 1 MHz, somewhat above most 
Tesla coil operating frequencies. He restricted himself to cases where the skin depth is a small 
fraction of the wire diameter. 


We might define a Medhurst resistance Ry, similar to the Medhurst capacitance Cy, 
where 


Ryu = Rac(l tr ks(om ow 1)) (6.36) 


where ky is a monotonic function of frequency that is zero for very low frequencies and unity 
for very high frequencies. That is, 
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Ru = 6m Rac f large (6.37) 


and 


Ru = Rac © Rae f small (6.38) 


Medhurst makes no attempt to find kr. Therefore we will have to restrict ourselves to 
the case of high frequencies, where the proximity effect is fully ‘saturated’, or where the wire 
radius is several skin depths thick. $j, is given in his Table VIII. A portion of that table is 
shown in Table 6.3. 


Table 6.3: Experimental values of $7, the ratio of high-frequency coil resistance to the resis- 
tance at the same frequency of the same length of straight wire. 
LetD) 
d/z, 1 2 4 6 8 10 ow 
1 5.00 4.10 3.54 3.31 3.20 3.23 3.41 
0.9 4.10 3.36 3.05 2.92 2.90 2.93 3.11 
0.8 3.17 2.74 2.60 2.60 2.62 2.65 2.81 
0.7 2.47 2.32 2.27 2.29 2.34 2.37 2.51 
0.6 1.94 1.98 2.01 2.03 2.08 2.10 2.22 
0.56 1.67 1.74 1.78 1.80 1.81 1.83 1.938 
0.4 1.45 1.50 1.54 1.56 1.57 1.58 1.65 
0.38 1.24 1.28 1.82 1.34 1.34 1.385 1.40 


In this table, ¢,, is the coil winding length, D is the coil diameter, d is the diameter of the 
copper wire, and z 1 is the center-to-center spacing between adjacent turns, all in consistent 
units. 


This table indicates that the proximity effect can easily double or triple the measured 
input resistance over that predicted by R,, for a straight wire of the same length. In the 
previous subsection, a 14 gauge coil was mentioned which had a dc resistance of 3.99 2, and 
an ac resistance of 10.85 2 at 160 kHz. By interpolation in Table 6.3, ¢jy is found to be 
1.85. Therefore, the predicted ohmic resistance of the coil would be (1.85)(10.85) = 20 Q 
(at sufficiently high frequencies). The measured input resistance is about 23 Q at resonance. 
This measured resistance includes dielectric losses and eddy current losses in addition to ohmic 
losses. The difference of 3 2 would represent eddy current, dielectric, and transmission line 
losses, which is probably not far from reality. If the proximity effect is not fully ‘saturated’, 
the predicted resistance would be less, and the difference between predicted and measured 
would be greater. 
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Poynting 


I then started thinking about methods to explain and perhaps even calculate ¢),. Note that 
we are still operating in the circuit theory mode here, such that the current is the same in 
every turn of the coil. The effect of distributed capacitance, which causes the current to be 
different in different turns will be discussed later. I worked on two different methods that I 
had not seen in the literature. One was to use Poynting’s vector and calculate power flow into 
the copper of the coil from the total magnetic and electric fields. The other was to split the 
current into filaments and require a distribution of filament currents that would minimize the 
magnetic field in the center of the wire. Both methods worked to some degree if ¢,,/D was 
not too short. After discussing these two methods, I will present some results from a recent 
paper [2]. This paper is computationally superior to my two methods, but still is not all that 
great for short coils. I have come to realize that finding the resistance of a short coil is a tough 
problem. I am tempted to take up golf or something less frustrating! 


The Poynting vector is defined as 


S=ExH W/m’? (6.39) 


where x refers to the cross product of two vectors. 


A student in Circuit Theory I will hear the instructor say that power flows down a wire 
inside the wire (electrons bumping along). The student will walk next door and hear the EM 
Theory instructor say that obviously power flows down the outside of a wire in the form of 
electric and magnetic fields. The speed of propagation is determined by the dielectric constant 
of the insulating material, not by any property of the conductor, hence the action must be 
happening outside the copper. The wire gets hot due to some leakage of the fields from outside 
to inside through the mechanism of Poynting’s vector. 


Most students compartmentalize this information, so think power flows one place in one 
class and the opposite place in another class, and rarely ask where the power really flows. The 
problem is that no one really knows. Correct answers are obtained with either approach. 


There is a problem, in that Poynting’s vector does not appear to work in some cases. 
Consider the following example. 


Example 


On a clear day the average electric field is 130 V/m at the surface of the earth, directed down 
(the earth is negatively charged with respect to space). At the earth’s magnetic equator, the earth’s 
magnetic flux density B = wH is horizontal, with magnitude about half a gauss, or 0.5 x 1074 T. If 
E is directed ‘down’ and H is directed ‘south’, then the Poynting vector S is directed ‘west’. The 
magnitude is 


EB 1 : 10-4 
S=FH = _ 1800.5 x 10™) _ gong W/m? 
Lo dn x 10-7 
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A power density of 5200 W/m? would be very noticeable, but our senses tell us that nothing is 
actually flowing. No one has built a machine to extract and sell this power. A student who asks about 
this problem is likely to get one of two responses. The first is a raised voice and some subtle ridicule for 
asking such a dumb question. The second is a variant of “Trust me. Sometimes the electron acts like 
a wave, sometimes like a particle. We who have already been initiated into the scientific priesthood 
will train you to know when it acts like what.” But I digress. 


Standard wisdom among the priesthood is that Poynting’s vector works when E and H 
are the cause and effect of each other, or are both produced by the same source. If current 
is flowing in a wire, it produces a magnetic field around the wire and a voltage drop (electric 
field) along the wire, so this is a proper application. 


For a straight conductor centered on the z axis, the magnetic field H is given by 


ay (6.40) 


~~ Onr 
where J is the current flowing in the conductor, r is the distance from the z axis, and ag is 
the unit vector in the ¢ direction around the conductor. The current J produces a voltage 
drop along the conductor and the electric field is E = J R’a, where R’ is the resistance per 
unit length. The cross product a, x ag = —a,, so the Poynting vector is directed into the 
conductor. 


The magnitude of the Poynting vector entering a straight wire of radius 0 is 


I 
S, = EH =IR' 


a (6.41) 


The total power entering the conductor is determined by integrating the Poynting vector 
over its surface. If the wire length is @ and the radius is b, the power is 


27 
P= fi [i IR')(s—)bdedz = PRE (6.42) 
which is exactly what is predicted by nae: circuit theory. 


The magnetic field of a filamentary loop, carrying a current J, centered on the z axis and 
located in the zy plane is [7] 


I a 
Ly = 7 (6.43) 
Qn pat p +z =p) 2 
I 1 i pF =e 
A, = — ——————_ | Kk +. — SE 6.44 
“On Jat p+? (a—p)? + 2? ee 


where K and E are the complete elliptic integrals of the first and second kind. The general 
structure is shown in Fig. 6.2, which shows three turns out of an N-turn Tesla coil. Coil 
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diameter is D and the radius of the individual wires is b. The distance between adjacent turns 
is z,; and the overall winding length is @,,. Current is flowing into the paper on the right, and 
out of the paper on the left. The electric field in each conductor is in the same direction as 
the current. 


Fi Hy; Hei 
ei 
1 


oS 


— 


pli 

© 
es 

Figure 6.2: Three Turns of Tesla Coil 


The elliptic integrals K and E are functions of a parameter k, where 


k2 = dap 


=Ge (6.45) 


One can either look up the values for K and E in a math table, or can evaluate their 
infinite series representation. The series converges relatively slowly, so computation of the 
first one or two thousand terms is not unreasonable, using double precision numbers. 


As might be expected, people have developed techniques for calculating these elliptic 
integrals with less computational effort. Gauss’s method using the arithmetic-geometric mean 
method converges very fast, usually in less than 20 terms for 10-digit accuracy [1] . The 
incremental H field on the inside surface of turn 7 that is produced by current J in turn 7 is 
given by 


2 2. 7 V2 2 
Aw; = ee + at ~ (a—b)* ~~ 9)" 5 (6.46) 
2m ./(2a — b)? + (i— j)222 
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and similarly for the field on the outside of the coil, AH,2;. The total field H,1; is found by 


summation. 


N 
Hai = >> AA ai 
j=l 


(6.47) 


Similar expressions hold for H,;, Hp1i, and H,2;. These four fields are tabulated in 
Table 6.4 for a 21 turn coil with ¢,,/D = 1 and 2b = z;. This corresponds to the upper left 


corner of Table 6.3. 


Table 6.4: Predicted values of H.4;, Hz2i, Hpi, and Hp; in A/m for a 21 turn coil with 
ly/D=1 and d/z; =1. D = 8 inches and d = 2b = 0.4 inches. 


Ay; 

65.16 
73.92 
77.62 
79.92 
81.58 
82.82 
83.76 
84.46 
84.93 
85.21 
85.31 
85.21 
84.93 
84.46 
83.76 
82.82 
81.58 
79.92 
77.62 
73.92 
65.16 


moEoewmntantkKRwWwnNe|+ 


NON FPRPRP re eRe re 
FOO ON MD OB W WD 


F,24 
-19.79 
-23.69 
-23.76 
-23.19 
-22.52 
-21.90 
-21.36 
-20.94 
-20.64 
-20.46 
-20.40 
-20.46 
-20.64 
-20.94 
-21.36 
-21.90 
-22.52 
-23.19 
-23.76 
-23.69 
-19.79 


Api 
-66.47 
-35.09 
-24.72 
-18.60 
-14.32 
-11.06 

-8.42 

-6.20 

-4.26 

-2.49 

-0.82 

0.82 

2.49 

4,26 

6.20 

8.42 

11.06 

14.32 

18.60 

24.72 

35.09 


Api 
-35.09 
-24.72 
-18.60 
-14.32 
-11.06 
-8.42 
-6.20 
-4.26 
-2.49 
-0.82 
0.82 
2.49 
4,26 
6.20 
8.42 
11.06 
14.32 
18.60 
24.72 
35.09 
66.47 


We see that the magnetic field is up on the inside of the coil and down on the outside. 
Likewise, it is directed radially inward over the bottom half of the coil (turns 1-10) and 
outward over the top half of the coil. The z component is largest in the middle of the coil, 
while the p component is largest at the ends. 


The total power flowing into the hollow cylinder bounding the coil is then 
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P= J Salarea) = (Area)(E)(5H) 


N 
€,20(IR’)(Hyt — Hyp) + &aa(TR')(9 "(Hai — Aeai)) (6.48) 
i=1 


(6.49) 


where ¢; is the length of one turn and I have added the magnetic fields H,; and H,» to include 
the end effects at the top and bottom of the coil. The minus sign in front of the H,9; term is 
necessary to indicate that the Poynting vector is into the copper on the outside of the coil as 
well as on the inside, and likewise for the H,» term. 


The power flowing into the same length of straight wire would be 


Pay = (Avea)(B)(H) = (N)(C,)(2n0)(TR\(5>) (6.50) 


Power is proportional to resistance so the ratio of P to P. is another ¢, call it ¢, to 
distinguish it from the Medhurst ¢. The ¢p is 


_ 42b(LR’)( Apt — Hos) + fz, (IR')(M, (Hai — Hi) 


6.51 
b Né,2nbI R'(I/27b) Ger 
Canceling the common terms gives 
yt — Hop + (Aa =A) 
SS Eo 6.52 
For the case shown in Table 6.4, 
2273.92 
dp = Rca e = 3.456 (6.53) 


(21) (31.33) 


We see that this approach yields an increase in resistance of the coil as compared with 
the same amount of straight wire of 3.456, as compared with the Medhurst prediction of 
5.55. This is within a factor of two, as mentioned earlier, but is not really close enough 
for any definitive computations. After spending considerable time looking at the problem, 
I decided that one significant flaw in the analysis is that it assumes the actual current flow 
throughout the conductor cross section can be modeled by a filamentary current (of the same 
total magnitude) flowing exactly in the center of each conductor. The actual current density 
is higher near the surface of the conductor due to skin effect, but will not be symmetric inside 
to outside, or top to bottom, hence the equivalent current filament will not be in the exact 
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center. The computations are sensitive to the position or location of the current, and if this 
is not accurately known, it will not be possible to obtain highly accurate results with this 
Poynting vector approach. 


Boundary Conditions 


My next analytic attempt was to split the current flow into four filaments, located on the 
surface at the top, bottom, inside, and outside of each turn. The magnetic field components 
H, and H, were calculated at the center of the wire cross section, using Eqns. 6.43 and 6.44. 
The computations are very similar to those of the previous Poynting vector effort. Instead of 
finding the fields on the wire surface from a current at the wire center, I find the fields at the 
wire center from current filaments on the surface. 


In the high frequency limit, the magnetic field at the center of a conductor will be zero. 
Currents will distribute themselves to meet this boundary condition. The problem then be- 
comes one of iteration to find a current distribution that will cause the magnetic field to be 
zero at the center of the conductor. For example, suppose we start with a total current flow 
of 4 A, or four filaments each carrying 1 A. We calculate the sum of H, and H, at the center 
of each turn. We then move a portion of the current in one filament, say 0.1 A, to another 
filament, and recalculate the sum of the fields. If the sum is smaller, we are headed the right 
direction, and try moving some more current from the first filament to the second. When we 
reach the zero field condition, we have found a solution to Maxwell’s Equations. 


One obvious problem with this simple algorithm is that the limit is found when all the 
current flows in one filament. If we think of the wire as four quadrants in parallel, but no 
current is flowing in three of the quadrants, it is like removing three out of four parallel 
resistors. The resistance of one quadrant is four times the resistance of the entire wire. This 
approach will yield a maximum ¢y, of 4.00, not the 5.55 appearing in the upper left corner of 
Table 6.3. One can get past this limit by splitting the conductor into more than four filaments, 
say 8 or 16, or by allowing negative current in some of the filaments. 


I tried 8 filaments, but this did not immediately fix the problem. It did not appear to be 
possible to get to the case of zero magnetic field in the center of the wire while requiring all 
currents to be non negative. If we allow currents to be negative, such that current is flowing 
in one direction on one side of the conductor and in the opposite direction on the other side, 
there must be some lateral or azimuthal flow of current. This violates our original assumption 
of current flowing only in the direction of the conductors. My mother told me there would be 
days like this! 


While this is certainly an interesting problem in its own right, it has only limited appli- 
cation to Tesla coils which usually have a length/diameter ratio on the order of four. The 
maximum @yy is now 3.54, from Table 6.3, which can be analyzed with four filaments and non 
negative currents. Results of the analysis are shown in Table 6.5. The column labeled dj, 
is from Table 6.3 while the column labeled ¢ 7 is for the predicted resistance ratio from this 
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iterated computer program. 


Table 6.5: Comparison of @,, with ¢, for a coil with ¢,,/D = 4 
d/zi om bs 
0.9 3.05 3.11 
0.8 2.60 2.63 
0.7 2.27 2.30 
0.6 2.01 1.99 
0.5 1.78 1.80 
0.4 1.54 1.60 
0.3 1.382 1.40 


It can be seen that agreement is quite good between the Medhurst measured values and 
the calculated values using the four-filament approach. I feel that the Medhurst table has been 
theoretically validated. I am sure that there are other methods of theoretically determining 
om, perhaps easier and more accurate, but this approach convinced me that proximity effect 
does indeed have a classical explanation, if one wants to spend the time and effort necessary 
to find it. In the meantime, Table 6.3 should be quite adequate for anyone looking for an 
estimate of the high frequency resistance of a coil of wire. 


Fraga 


After preparing the above material, a reviewer [1] informed me of a paper by Fraga, Prados, 
and Chen on this topic [2]. An examination of the paper indicated that the authors seemed 
to be restricting themselves to the following case. 

1. Long solenoids 

2. Tight wound solenoids 

3. Solenoids with negligible distributed capacitance 

4. Low frequency (b < 6) 

5. Multilayered coils 

A typical Tesla coil has a length/diameter ratio of about four (not necessarily long). Some 
coils are space wound. Distributed capacitance is always a problem. It will be discussed in 
the next section. Tesla coils are usually operated at frequencies where the wire radius will be 


between one and five skin depths, not less than a single skin depth. And Tesla coils are almost 
never multilayered. So it was not obvious that the paper would be of much use to the Tesla 
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coil community. On the other hand, Tesla coilers are infamous for using things outside their 
normal operating range (e.g. pulling 60 mA from a neon sign transformer rated at 30 mA), so 
I looked some more. The authors have developed a relatively simple closed-form expression 
for the resistance of a coil, not requiring summations or interpolations, so it would be very 
nice if it worked for the typical Tesla coil. 


Their Eq. 35, expressed in my notation for the single-layer coil is 


R.- 27 N2pegalsinh 20 + sin 26] 
* 8 s£y [cosh 20 — cos 26] 


where N is the number of turns, a is the radius of the coil, and ¢,, is the winding length of 
the coil. The other terms are described in the following. 


(6.54) 


One of the difficulties in the analysis of a coil is that the conductor surfaces are not at a 
fixed distance from the axis of the coil. A fingernail pressed against the coil moves in and out 
as the hand moves down the coil. Boundary conditions cannot be simply expressed in this 
geometry. These authors circumvent this problem by converting round wire to square wire with 
the same cross-sectional area. They deal with the gap between the wires, due to insulation, 
by elongating the square conductor toward its neighbor until it touches mechanically (but 
not electrically). The resistivity is increased a proportional amount such that the resistance 
remains the same. 


Let the radius of a cylindrical wire be b, covered by a dielectric coating of thickness s. 
Define a square conductor with side y with the same area. Then 


y? = 1b? (6.55) 


or 


y=V/nb (6.56) 


The length of one turn is 2)+ 2s. The effective resistivity for this rectangular wire of 
thickness y is 


J. 2OA23 "2p 
Peff = P y mvs 


where p is the usual resistivity of the conductor (1.724 x 10~§ ohm meters for copper at 20°C). 


Ss 


(a+ b? (6.57) 


They then define an effective skin depth 


= {2a = [Ae 8/0) _, [A174 x 1051+ 9/8) _ 4 ray, ft a/b 
oe in V VamoQaf) \ vata x10 )amp OMY 5 (6.58) 
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for copper. The last variable, 6, is then defined as 


= (6.59) 


One big advantage of this approach is that both the skin depth and the proximity effect are 
included in one equation that would fit on most programmable calculators. Input quantities 
are wire radius, coil radius and length, number of turns, and frequency, all readily available. 
Using the Medhurst approach outlined earlier requires one to first find the ac resistance and 
then multiply this by a factor @,4 found by interpolation in Table 6.3. 


A disadvantage of the Fraga method is that the resistance is equivalent to the last column 
in Table 6.3. That is, we find the resistance of a section of an infinitely long coil. This method 
is not capable of properly dealing with the short coil. A comparison of dj;yg and what I call 
or = Rs/Rac is given in Table 6.6. The Fraga formula was evaluated for 18 gauge wire at 1 
MHz to make it as comparable as possible to the Medhurst experiment. 


Table 6.6: Comparison of ¢)y with ¢p for a coil with ¢,,/D = oo 
d/z1 ou oF 
10 3.41 3.19 
0.9 3.11 3.03 
0.8 2.81 2.86 
0.7 2.51 2.67 
0.6 2.22 2.48 
0.5 1.93 2.26 
0.4 1.65 2.03 
0.3 1.40 1.75 


We see that the Fraga formula agrees best with Medhurst for a wire diameter to wire 
spacing ratio of between 0.8 and 0.9, which happens to be the typical ratio for tight wound 
magnet wire. Therefore, there is hope for the Fraga formula for typical Tesla coils (tight 
wound magnet wire with a length/diameter ratio of at least four). 


I then proceeded to compare the coil resistance R, predicted by Fraga with the coil resis- 
tance Ry predicted by Medhurst and with the measured resistance Rro for my coils. Results 
are given in Table 6.7. 


This table gives the wire diameter d = 2b in both mils and mm, the coil radius a, the 
winding length @,,, the total length of wire used in the coil, the number of turns N, the 
Wheeler inductance L, the Medhurst capacitance Cy, the Medhurst factor ¢,y interpolated 
from Table 6.3, and the dc resistance Rg. at 20°C. Measured inductance is always close to 
the calculated value from Wheeler’s formula. Measured resistance agreed with the calculated 
value to within 1% or so, when corrected for temperature. 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 6—Resistance of Coil 6-21 


Table 6.7: Predicted and Measured Coil Resistance for Several Coils 


12T 145 144T 16B 18B~ 18T 20T 22T = 22B 


2b, mils 80.81 64.08 64.08 50.82 40.30 40.30 31.96 25.35 25.35 
2b, mm 2.053 1.628 1.628 1.291 1.024 1.024 .8118 .6439 .6439 
a, meters 231 =.198 = .107 .24 235.107 107 107 .24 


fw, meters 1.613 1.166 1.392 .7 AT 881 .952 945 ok 
wire, meters 509 482 623 675 617 534 702 424 675 


turns NV 301 387 797 445 418 794 1052 631 445 

L, mH 14.2. 17.2 19.34 49.4 55.02 29.1 47.3 17.2 49.4 

Cm, pF 30.58 23.95 20.70 22.66 21.22 15.45 16.14 16.08 22.66 
om 1.68 1.85 3.03 420 4.25 3.15 3.02 1.62 1.54 
Rae 2.66 3.99 445 889 12.9 11.2 23.4 22.4 35.7 
fo, kHz 246.4 247.6 251 153.5 147.3 236.9 181.1 301.5 153.5 
Rs 22.83 29.99 44.11 51.13 62.49 67.55 97.33 65.94 69.00 
Ru 17.96 24.7 46.35 80.93 95.61 76.50 111.39 58.38 66.09 
Rro 18.7 245 43.5 93.1 118.0 70.5 94.2 47.0 73.0 
f 217.38 211.3) 128.4 123.6 176.1 1385.9 227.9 129.5 
Rs 28.10 40.48 46.76 54.04 58.23 84.47 56.76 62.23 
Ru 23.13 41.67 74.97 88.92 66.52 98.37 51.57 55.00 
Rro 25.6 42.3 87.1 103.0 65.9 88.0 46.7 69.5 
f 158.8 151.8 91.9 86.2 122.55 94.2 158.1 93.1 
Rs 24.02 34.31 39.56 46.90 48.57 70.40 45.76 51.74 
Ru 20.06 35.86 64.76 76.65 57.03 85.04 44.56 55.00 
Rro 24.5 39.6 75.0 78.9 58.1 78.7 42.4 66.8 


The suffix S refers to space wound, while T refers to a tight wound coil. The suffix B 
refers to a barrel whose sides are not perfectly straight, for coils 16B and 22B. The winding 
is tight wound on the flat portions of the barrel and space wound on the transition portions. 
The barrels were assumed to be made of polyethylene when I purchased them at the local 
recycling plant. The barrels used for 12T and 18B are straight sided and have thinner walls 
than the 16B and 22B barrels. They were once used as tanks for water softeners. Coil 145 is 
on a coil form built from a 0.125 inch sheet of polyethylene. Coils 14T, 18T, 20T, and 22T 
are on PVC forms. Coil 20T has 3 layers of polyurethane on it, the others have no coating. 
Coil 12T is made in two sections for ease of handling. 


The 12 ga wire is Essex type USE-2 (or type RHH or type RHW-2), which are types 
specified in the National Electrical Code. This type has a relatively thick insulation, which 
spaces the conductors farther apart than the thin insulation of magnet wire. It has a nominal 
insulation thickness of 45 mils as compared to a 15 mil thickness for the more common Type 
THHN. The 14, 16, 18, and 20 ga wires are magnet wires coated with Heavy Soderon. This 
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Essex coating has a top layer of nylon. Magnet wires are available with one, two, three, and 
four layers of insulation. Heavy Soderon is equivalent to a two layer coating. The 22 gauge 
wire is not magnet wire but has a yellow insulating jacket, that I assume is PVC. The wire 
was acquired surplus and the numbers are difficult to trace. 


The Medhurst capacitance and the Wheeler inductance are used to calculate the resonant 
frequency fo. The actual resonant frequency was measured with a HP54645 digital scope. 
Agreement was within 5%. The operating resonant frequency will always be lower than the 
unloaded value due to the toroid on top, so this number is of mostly academic interest anyhow. 


The observed resistances include displacement current effects and any dielectric losses, 
which might increase the actual resistance well above that due to copper losses alone. The 
highest frequency in the table refers to the case of no top load, while lower frequencies were 
obtained with different sizes of toroids mounted on top the coil. 


Coil 14T has the greatest length/diameter ratio (6.47) of the group, so we would expect 
Fraga’s formula to work the best for this case. Indeed, it predicts a resistance only about 4% 
less than that predicted by Medhurst, both close to the experimental values. 


As we move to coil 18B, which has the smallest length/diameter ratio (1.0), Fraga’s formula 
significantly under predicts the resistance, being about 55% of observed and about 63% of the 
Medhurst values. 


On the other hand, coil 22T has a reasonable length/diameter ratio (4.42), but greater 
dielectric thickness s, and Fraga’s formula over predicts the resistance, being about 124% of 
observed and about 109% of Medhurst. 


For the seven coils 14T—22B, Fraga’s formula predicts resistance values about 90% of what 
Medhurst predicts. All things considered, this is pretty good. If we stick with coils close 
wound with magnet wire, and a length/diameter ratio of 4 or more, Fraga’s formula should 
be quite acceptable. 


The next step was to determine the character of the transition of resistance values as 
frequency is raised from dc to those frequencies where the proximity effect is fully implimented. 
That is, there is no proximity effect at dc. The measured resistance of a coil is the same as the 
same length of straight wire. As frequency increases, however, the skin effect causes resistance 
to increase, and the proximity effect causes resistance to increase even more. Presumeably, 
at a high enough frequency where the skin depth is a small fraction of the wire radius, the 
proximity effect saturates. The measured resistance continues to increase, but only due to the 
skin effect. 


I measured the input resistance as a function of frequency for seven of my coils, 14T, 
16B, 18B, 18T, 20T, 22T, and 22B. This includes all five coils with a Medhurst factor jy = 
Rrc/ Rac of 3.0 or more. I used a standard bench function generator with sine wave output, 
followed by a linear amplifier. The amplifier is a simple single-stage inverting op-amp, using 
the Apex PA-19, rated at 4 A, +36 V, slew rate 900 V/s. Voltage was measured with 
a standard scope probe, current with a Philips PM9355 current probe. The output of the 
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current probe was fed to the second channel of the scope, a HP54645D digital oscilloscope 
capable of calculating rms values of measured waveforms. I applied a voltage between about 
2.5 and 10 V (rms) to the base of the coil, tuned the function generator for resonance, observed 
the current value, and calculated the resistance as the ratio of voltage to current. 


The test location was inside a 54 by 90 ft metal building (manufactured by Morton) that 
is typically used for livestock or storage of ag equipment. One 15 by 30 ft corner was framed 
in, insulated, and equipped with furnace and air conditioner. Inside this instrument room 
was a double copper wall Faraday cage with footprint 8 by 12 ft. Electricity was provided 
to this screen room through an isolation transformer and power line filters. The computer, 
oscilloscope, and other sensitive equipment were located in this screen room. It is interesting 
to note that reception on a transistor AM radio was not affected by moving it from outside to 
inside the Morton building, but reception was impossible inside the screen room, even with 
the massive copper door open. The AM band starts at 550 kHz, slightly above most of my 
testing in the 100-300 kHz range, so it would appear that the Faraday cage was effective in 
the necessary frequency range. I never observed any failures in electronic equipment inside 
the cage that I could link to transient high fields outside the cage. 


This Morton building has a wood frame and wood trusses supporting the roof. There is 
no ceiling so the wood trusses and roof are visible inside the building. The bottoms of the 
trusses are about 16 ft above the floor. Depending on the location, the roof will be about 17 
to 25 ft above the floor. 


Except for the screen room, all electrical outlets in the building were wired in the con- 
ventional manner for North America, with the third wire connected to utility ground. This 
ground wire is connected to earth at every power pole on the utility system. I installed my 
own ground system under the dirt floor of the Morton building, consisting of three lengths of 
copper tubing buried a foot or so below floor level. The soil is very hydrophilic, so when I 
water the grass on the outside of the building, the copper tubing is located in wet earth. The 
measured resistance between this local ground and utility system ground was on the order 
of 1 2, a quite acceptable value. Open circuit voltage between the two grounds may be as 
high as several volts, and current flow between grounds may be as high as several amps. The 
utility ground is a major source of noise in sensitive measurements, so all measurements were 
made using only the local ground. The Faraday cage was connected to this local ground, as 
well as the metal skin of the Morton building. 


A damp earth floor made the interior of the building too humid, so I covered the floor 
with polyethylene sheets, and covered those with about four inches of milled asphalt. This is 
a product obtained when asphalt roads are recycled. If carefully packed, it forms almost as 
nice a surface as concrete, and is much less expensive. It is also much easier to penetrate if 
one wanted to install something in the earth. 


Even though the local ground is very satisfactory for making sparks with Tesla coils, I was 
concerned that variations in earth moisture would affect my proximity effect measurements. 
I therefore installed a metal ground plane on top of the asphalt millings, about 8 by 16 ft, 
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consisting of sheets of aluminum siding bonded together with copper flashing and sheet metal 
screws. This would obviously not be suitable for spark tests when hundreds of kV and tens of 
amps are present, but for low level measurements of a few volts and less than one amp should 
be ok. 


The coil under test was placed on a piece of 2 inch thick blue styrofoam on this metal 
ground plane. Location was about 10 ft from the wall of the Morton building and about 10 
ft from the wall of the instrument room. A length of 50 2 coaxial cable was run under the 
ground plane, through the wall of the instrument room, through the wall of the Faraday cage, 
and to the function generator. The shield of the coax was connected to the ground plane 
close to the base of the coil. The center conductor of the coax was connected to the base of 
the coil, so the coil was driven like a vertical helix above a ground plane. With no toroid or 
other connection to the top of the coil, the inductance of the coil would resonate with the 
Medhurst capacitance. This would be the maximum frequency at which the impedance could 
be measured. 


Capacitance could be increased and frequency lowered by placing larger and larger toroids 
on top the coil. However, it is not feasible to have toroids large enough to lower the frequency 
to a few kHz, as required for this test. Different size toroids also change the local field 
distribution, which might affect the results. So I decided to connect a capacitor to the top of 
the coil, and the other lead of the capacitor to a heavy ground wire that went up from the 
top of the coil to about 10 ft above the floor, over to the wall of the instrument room, and 
down to the floor where it was connected to both the metal ground plane and to the copper 
tubing of the local ground. Electrically, this forms a simple series RLC circuit. If large value 
capacitors are used, the resonant frequency can be under 1 kHz. 


Like many other things about Tesla coils, using the wrong type of capacitor can lead to 
surprises. The voltage rating is important, as a little thought will reveal. The Q of these 
coils is high, perhaps as high as 500. Driving the coil with 10 V would then put as much as 
5000 V across the capacitor. My first effort was to use variable capacitors here. I had a ham 
type variable capacitor that would handle the voltage, but receiving type variable capacitors 
would arc over between the plates. Even when the voltage was acceptable, the resistance of 
the sliding contacts in the capacitors was too high (and too erratic) at the required current 
levels. 


The final solution was to place an aluminum cake pan, flat bottom up, on top of the coil 
and electrically connected to it. A half spun aluminum toroid was connected to the hanging 
ground wire. Polyethylene sheets were placed on the cake pan and the half spun toroid would 
be placed on the sheets, flat side down, to form a simple parallel-plate capacitor. Different 
thicknesses would allow for the resonant frequency to be lowered to about one fourth of the 
maximum value. For even lower frequencies, lumped capacitors were used. 


Figure 6.3 shows the variation of Rrc/ Rac versus frequency for the coil 14T, plotted solid 
with diamond symbols. It also shows the Fraga resistance R, divided by Rac, plotted dashed. 
Theory and experiment match extremely well. Again, the Fraga formula works well for long 
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coils that are tight wound with magnet wire. 
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Figure 6.3: Rrc/Rac for coil 14T 


I then performed the same measurements of Ryo for coil 18B, the shortest coil in my 
collection. Results are shown in Fig. 6.4. The Medhurst @jy is 4.25 for this coil where ¢p 
is about 3. The predictions by Fraga are below the measured values over the entire range of 
frequencies. Even worse, the experimental and theoretical start to diverge at about 50 kHz. 


I have performed other tests which show that this particular magnet wire absorbs moisture, 
and this moisture will cause increased losses. Dielectric losses are entirely separate from the 
proximity effect. It appears that the dielectric losses become significant above 50 kHz for coil 
18B, causing the ratio Rrc/ Rac to go well above the value ¢jy predicted by Medhurst for the 
case with negligible dielectric losses. 


Figs. 6.3 and 6.4 also show another effect, a very interesting concept that is otherwise 
difficult to explain. This concept is that there is little penalty in performance if one uses a 
smaller wire in a coil. That is, the effect on spark length is not as strongly related to the wire 
resistance as one would expect. 


If one looks closely at these figures, it is evident that the transition region moves higher in 
frequency as the wire gets smaller. At a given frequency, the larger wire will always be closer 
to saturation. Consider an extreme example where two coils are each wound with 1000 ft of 
magnet wire, one with 14 ga and the other with 28 ga. Assume the coils are resonant at 70 
kHz and that $j, = 3. The dc resistance of the 14 ga coil is 2.525 2 while the dc resistance 
of the 28 ga coil is 64.9 2, a factor of 25.7 greater. The ac resistances at 70 kHz are, from 
Egs. 6.30 or 6.32, (1.888)(2.525) = 4.768 Q for the 14 ga coil and (1.0034)(64.9) = 65.1 Q for 
the 28 ga coil. The 14 ga coil has reached full saturation from the proximity effect at 70 kHz, 
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Figure 6.4: Rro/Rac for coil 18B 


so the effective Rrc is (4.768)(3) = 14.3 Q. The 28 ga coil has not started into the proximity 
effect yet at 70 kHz, so its resistance is still just 65.1 Q. The ratio of resistances at 70 kHz is 
65.1/14.3 = 4.55, a considerable reduction from 25.7. 


It would have been nice to finish this treatment of copper losses with a formula that was 
accurate to within 5% for any frequency, any length/diameter ratio, and any wire diameter 
to wire spacing ratio. But that remains for someone else. I hopefully have described the 
problem so that the reader can get within perhaps 20% of the correct value, and sometimes 
even better. 


6.4 Displacement Current Effect 


We have examined two methods for empirically or theoretically determining the copper loss 
in a coil, the methods of Medhurst [5] and Fraga [2]. These both assume that the conduction 
current is the same at all points in the coil, which is the usual case for circuit theory type 
RLC models. But this is not necessarily the case for a Tesla coil. This is one place where 
the lumped circuit models just cannot go. We finally have to use a distributed approach. 
If the conduction current is less than the input current in part of the coil, we would expect 
the effective resistance to also be less. Likewise, if the conduction current is greater than the 
input current, then we would expect the effective resistance to increase. Conduction current 
can be greater than the input current in one part of a coil, and less in another part, due 
to displacement currents. We will try to illustrate this concept with Fig. 6.5. A Tesla coil 
is connected to a toroid with a switch S;. The input current at the base is the conduction 
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current Zn. 


Every part of the coil has a capacitance to every other part. We show four capacitors in 
this four turn coil, Co; from turn 2 to turn 1, C42 from turn 4 to turn 2, Cy2 from the toroid 
to turn 2, and Cy, from turn 2 to ground. Each capacitor has a current flow in it, with the 
same subscripts. 


Figure 6.5: Conduction Current 


The incremental current entering turn 2 is 


Ai = igo + ty2 — 121 129 (6.60) 


If the two currents 749 and 2,2 are greater than the other two currents, then A? is positive. 
If the current in turn 1 is 7;,, and Az is added in turn 2, then the current in turn 3 is greater 
than 2;,. This situation is quite possible for the lower part of the coil where many turns above 
the turn in question are adding currents and only a few turns below are subtracting currents. 
It is more likely to occur if Cy, is small, that is, if the coil is mounted well above a ground 
plane. 


On the other hand, as we get toward the top of the coil, there are many turns below the 
turn in question that are subtracting currents from Az and only a few turns above that are 
adding currents. The minimum current occurs in the top turn of the coil. If switch $j is 
open (no toroid), this minimum current is zero. So we have the situation where the current 
increases in the lower part of the coil, will hit a maximum probably somewhere in the middle 
third of the coil, and then start to decrease toward the minimum current value at the top of 
the coil. 


Possible variations of conduction current in the coil winding as a function of position y is 
shown in Fig. 6.6. The current marked “.S; closed” is a possible current with toroid while the 
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current marked “S; open” is a possible current without a toroid. 


$1 closed 


0 top 


Figure 6.6: Conduction Current in a Coil 


There are several qualitative factors that can be deduced from these curves, without going 
to the effort to do the full distributed analysis. We started on this quest by asking “Is the 
lumped RLC model for the Tesla coil a useful concept, or must we proceed immediately to 
the distributed model?” If this distributed effect is not too big, such that any fudge factor is 
no more than a few percent, then we can still use the lumped model. If the effect is too big, 
then we will be inclined to use the distributed model in order to get decent results. 


We note that the current in the coil is greater than the input current in the lower part of 
the coil and less over the upper part of the coil. These tend to offset each other so that the 
effective resistance may not be much different from the resistance calculated for the uniform 
current case. 


As larger toroids are added, the current in the coil increases, and likewise the effective 
resistance. However, the resonant frequency decreases with larger toroids, which lowers the 
effective resistance for the uniform current case. This predicts a smaller change in resistance 
with frequency than predicted by either Medhurst or Fraga. 


The final qualitative factor is that if a coil has a geometry that causes the peak current 
in the coil to be well above the input current (say 30% or more), the effective resistance will 
also be well above that for the uniform current case. We can now go back to Table 6.7 and 
check out these predictions. 


We note that for coil 14T, the measured resistance is below both R, and Ry for a coil 
without a toroid and above both R, and Ry, for a coil with the largest toroid. The mea- 
sured resistance holds more nearly constant with variation in frequency than is predicted by 
Medhurst or Fraga, a very consistent observation throughout all my testing. 


It appears that as coils get shorter and fatter, the interior current in the coil gets larger 
and the effective resistance increases as compared with the predictions of Medhurst and Fraga. 
Consider coils 22T and 22B. Coil 22T (relatively long and thin) has a measured resistance 
below the predictions, while coil 22B (short and fat) has a measured resistance above the 
predictions. The other two short coils (16B and 18B) also have measured resistances above 
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the predictions. 


The average ratio of Rc over R, for the coils 14T, 18T, 20T, and 22T was 1.004, while the 
average ratio of Rrc over Ry for these four coils was 0.927. It appears to me that for normal 
Tesla coil geometries (length/diameter = 4 or more) that the uniform current assumption is 
not too bad. It appears to yield accuracies within +10%, which should be acceptable in most 
applications. We conclude that the displacement current effect is very real and easily observed 
in data sets like Table 6.7, but the errors involved in ignoring it are not so severe that we 
cannot use the lumped model. 


6.5 Dielectric Losses 


We turn now to losses in the coil form and in the wire insulation. Both the coil form and the 
wire insulation form a part of the coil capacitance. By Gauss’s Law the coil capacitance is the 
sum of electric flux lines leaving the coil, divided by the coil voltage. Some flux lines go from 
toroid to earth through the coil form, some from turn to turn through the wire insulation, and 
some through air. These can be considered as three capacitors in parallel. Since the volumes 
of the coil form and the wire insulation are much smaller than the volume of air, and the 
relative permittivity is only two or three times that of air, their capacitances will be a small 
part of the total. The losses may still be significant, of course. 


If we assume a linear increase of voltage along the coil, the flux lines in the coil form will 
be uniform from top to bottom (no fringing) just like the parallel plate capacitor. We should 
be able to use the formula for the parallel plate capacitor without great error. 


Dielectric losses are usually modeled by a resistor in parallel with the capacitance, rather 
than in series. They are related to the capacitor voltage rather than to the capacitor current. 
We have two dielectrics, the coil form and the wire insulation, so we have two resistors in 
parallel with the Tesla coil capacitance. We saw equations for the power dissipation in the 
coil form, Pf, and in the wire insulation, P,,;, back in Chapter 3. If the toroid voltage is Vior, 
then the parallel resistances can be defined as 


— lor (6.61) 


a (6.62) 


A somewhat more detailed RLC model is shown in Fig. 6.7. The resistances directly 
related to current are shown as Ry, Reddy, Rspark, and Rad, where the last three items are 
the equivalent resistances representing losses to eddy currents in toroid and ground plane, the 
spark itself (when present), and any losses due to radiation. We considered the Medhurst 
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resistance Ryy earlier in the chapter. You can replace Rj by R, from Fraga’s formula if you 
prefer. The resistances related to voltage are Re and Ryj. 


Figure 6.7: Detailed Lumped Model of Tesla Coil 


For single-frequency, steady-state operation, the parallel combination of a capacitor and 
two resistors can be modeled as a series capacitor and resistor, call it Ra. This is straightfor- 
ward Circuit Theory I, but a bit tedious. We write an expression for the parallel impedance 
of Ref, Rwi, and the capacitance, rationalize it, and simplify the real term. We assume that 
the parallel resistances are much larger than the capacitive reactance, as they will be for any 
coil with acceptable losses. The algebra goes as follows: 


ge 1 _ 1 _ G — jwC re (6 63) 
- ly Rg + 1/Ref + JWC'¢ 7 G + JwCre = (G + jwCre)(G Ts jwCre) ; 
We define Rae as the real part of Z. 
es G Pot + Pwi 
Radice = lZ= l | —=— > | = SS FF Ss 6.64 
te = Real Z = Real ( See) ~ ag = Vesa ee 


To make things even more complicated, we will define Rgje as the series combination of 
three resistors, the series coil form resistance R, fs, the series wire insulation resistance Ry»jis1, 
and the series coating resistance Ry»js2. If there is no coating on the coil (polyurethane or 
equivalent) and if the wire and atmosphere are dry then Ryjs2 = 0. 


Réie = Refs + Rwist + Rwis2 (6.65) 


The revised circuit model for the Tesla coil is shown in Fig. 6.8. 


We now proceed to get specific equations for these three series resistors, referring back to 
Chapter 3 for the power dissipated in terms of V;o,.. 


Rope = agit = VenwOs(DF ey _ Cop(DFhep oo 
cfs = = — 
Vie, CE, V2 ,W2Cz, WCF 
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Figure 6.8: Lumped Model of Tesla Coil with Series Rgj- 
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It would appear that Rae decreases as 1/w or as 1/f. This is certainly contrary to 
our intuition. What is even more surprising is that the dissipation factor of water is also 
proportional to 1/w over the typical Tesla coil frequency range. (DF) water = 0.396 at f = 10° 
Hz and 0.0396 at f = 10° Hz [3]. This would make Rgje vary as 1/w?. I have not convinced 
myself that this is correct, but can see no flaw in the above analysis. I will terminate the 
theoretical discussion of dielectric losses on that note. 


There are many things going on that make it difficult to be precise about frequency 
variation of the other losses, but generally speaking, Ra increases as \/f, and Reddy increases 
as f?. Rraq Will increase at a rate somewhere between f and f?. Rspark Can be ignored below 
the spark inception voltage. Depending on which terms are dominant loss terms, we may not 
see a pronounced change in input impedance with frequency. That has been my experience. 
Input impedance will drift from day to day, (mostly with humidity), but there is no obvious 
frequency dependence. Of course, other things are happening. We know that R,, increases 
with temperature, while Rg. increases with humidity. If these were the only factors, we 
would expect a cold, dry winter day to have the lowest impedance, and a hot, muggy day to 
have the highest impedance and the worst performance. In cases where moisture is a factor, 
performance might improve after a period of operation which caused the coil form to heat up 
and dry out. 


Moisture has been a very frustrating factor in my testing. Coils are located inside a metal 
skin building with no climate control, in eastern Kansas. Temperatures can vary from below 
freezing to 40°C (104°F) or more. Relative humidities vary from 25% to 100%. Typical of 
my measurement problems are two sets of input impedance data in Table 6.8 for 3/17/01 
and 4/6/01. The 3/17/01 data were collected when the bay temperature was about 9°C and 
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the relative humidity was about 28%. On 4/6/01 the temperature was about 17°C and the 
relative humidity was 100%. It had been damp all week with heavy fog the day before. 


Table 6.8: Rro Measured on Two Different Days 


Coil 1445 14T 16B 18T 18B 20T 22T = 22B 
frequency, kHz 249.1 266.1 145.5 242.6 145.6 183.6 307.0 148.3 
Rrc 3/17/01 24.55 43.5 93.1 70.5 94.2 47.0 73.0 


Rro 4/6/01 27.6 45.5 175.8 75.7 127.7 100.0 53.4 126.7 


We see that two of the coils experienced large changes in Ryo, coils 16B and 22B. Both 
coils used a plastic barrel as a coil form that I thought was polyethylene. I got the barrels at 
the local recycling plant. Coil 22B used the same type of wire as coil 22T which was wound 
on a piece of PVC, so the difference in Rrc between these two coils had to be the coil form. 
These results indicate that some coil forms are worse than others. These barrels evidently 
soak up water in amounts sufficient to raise the input impedance by a factor of two. 


Coils 14T, 18T, 20T, and 22T were wound on PVC while coils 145 and 18B were wound on 
polyethylene. Only coil 20T had any type of coating put on top the winding (polyurethane). 
Both PVC and polyethylene appear to have about the same increase in Rrco with humidity, 
so it is hard to argue that one should spend more money on the more expensive polyethylene. 


As long as one stays with good quality coil forms, it appears that high humidity will raise 
Rrc by 5-10% from the low humidity case. One could reverse engineer Eq. 6.67 and find an 
effective dissipation factor (DF); that would be an appropriate function of humidity, but I 
am not sure it would be worth the trouble. 


The eddy current loss will be a strong function of how near the conducting material is 
located to the coil. A coil sitting on a ground plane would have a much larger Reddy than 
one sitting on a one meter high stack of Styrofoam blocks. If soil moisture affects the eddy 
current loss in the earth beneath the coil, then this term could vary widely from day to day. 
If tests are being done inside a metal building, then the walls and roof of the building would 
contribute to the eddy current loss. 


Many coils have a strike ring located around their base, to intercept sparks before hitting 
the feed line or other components. There is general agreement in the Tesla coil community 
that this ring should be open rather than shorted. There have been observations where a 
shorted copper ring has significantly degraded spark length. A spun aluminum toroid also 
presents a shorted path for eddy currents, but there is less agreement that this represents 
a significant loss. It is argued that conduction currents are smaller at the top of the coil, 
therefore induced currents must be less. 


I built a toroid of 0.25 inch copper tubing pieces on insulating disks, connected together at 
one point by a conducting disk. The ends were placed into heat shrink tubing, which was then 
shrunk to hold the ends a fixed small distance apart. This toroid was then compared with a 
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spun aluminum toroid of similar capacitance, and also with a smaller toroid made of one inch 
copper tubing with diameter slightly greater than that of the coil form. The smaller toroid 
was an attempt to get a shorted turn as near to the coil as possible. It lacked the capacitance 
to be an effective toroid for long sparks, of course. I could not find any significant difference 
in input impedance between the insulated toroid and the spun aluminum toroid. The shorted 
copper ring, however, had about 10% higher input impedance than the toroid that was not 
a shorted turn. This suggests that you would not notice any improvement if you cut your 
beautiful spun aluminum toroid into pieces to eliminate eddy currents. The effect is there, 
and can be measured if one really works at it, but is not that significant in most situations. 


Overall, my tests indicated that Reaay is no more than a few percent of Rrc. If a little 
thought is given to separation of conducting materials from the immediate vicinity of the coil, 
eddy current losses can be ignored. Likewise in all my tests, R;aq is very close to zero. I was 
unable to detect a signal from the coil more than perhaps 100 m away. At worst, it would be 
a number like 0.01 2, which is a negligible portion of a typical measured resistance of 25 to 
50 2. 


6.6 Conclusion 


I believe that Fig. 6.8 is a reasonable model for a Tesla coil. There is scientific basis for 
calculating (or estimating) Ry, Refs, Rwisi, and Rwisg, and for ignoring Reddy and Rrad- 
It has the proper indication for changes in the model when a spark occurs. The model 
indicates that when a spark occurs, the equivalent resistance Rspar, increases from zero to 
some finite value, so the input resistance increases during a spark. This is exactly what 
happens experimentally. 


Unfortunately, great precision is difficult to impossible to obtain. Rjy or Rs, can be cal- 
culated to within a few ohms given the techniques in this chapter. I consider this a vast 
improvement over my state of knowledge when I started this project. Trying to get more 
accuracy is probably not warranted because of the strong influence of moisture on coil resis- 


tance. If Ry is 50+5Q, and Rage might vary from 0 to 5 Q or more as humidity goes from 0 
to 100%, there is little point in reducing the uncertainty on Ryy. 


In my opinion, a complete distributed model will not be any better in dealing with skin 
effect, proximity effect, and dielectric losses, and would certainly be more of a programming 
problem. The one thing that this lumped model cannot deal with directly is the displacement 
current effect. A distributed model can determine the actual current distribution, which can 
then be used to find a predicted effective resistance of the coil. 


Both approaches (lumped and distributed) have advantages. I believe the lumped ap- 
proach is better at determining resistance. However, the lumped approach will not show 
anything about resonances at harmonic frequencies, and cannot deal with things like the cur- 
rent distribution. Hopefully there will be peaceful coexistence, where each method will be 
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used to its full advantage. 
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TESLA COIL DRIVER 


In this chapter, we deal with the design of the Tesla coil driver. We will discuss several 
components or subsystems, including the power supply, the controller, the gate driver and 
inverter, and the current sense resistors. There are a limited number of interconnections 
between the subsystems. We will try to be clear as to what the inputs and outputs are for 
each subsystem. 


The previous six chapters were written in 2001 and deal with the theory of Tesla coils. 
Some experimental results of a number of coils, built with different wire sizes, diameters, and 
heights, are included. The original Chapter 7 described my Tesla coil driver as it existed in 
2001. It worked, but I thought some improvements were needed. The humidity in my Morton 
building near Manhattan was just too high to do serious Tesla coil research, so I moved my 
lab to a small commercial building in Canon City, Colorado in 2008. The building is adequate 
in most ways, but the 10 foot ceiling does not allow for long sparks. My 30 inch aluminum 
sphere mounted directly on top of coil 14S is only 36 inches from the ceiling. I will probably 
limit the power so the sparks do not hit the ceiling. 


7.1 Power Supply 


The power circuit is shown in Fig. 7.1. We start with 240 VAC at 60 Hz in the upper left of 
the power circuit and end with the series RLC circuit model for the Tesla coil in the lower 
right. The idea is to apply a square wave voltage v;. at the resonant frequency to the base 
of the Tesla coil. The coil is basically operated as a quarter-wave antenna above a ground 
plane. The Q of coil 145 is quite high, perhaps on the order of 700 at resonant frequencies 
around 150 kHz. A square wave voltage input results in a sine wave of current that grows 
with time. Each cycle of the square wave input adds energy to the electric and magnetic 
fields of the Tesla coil, and increases the voltage of the coil top load. When the voltage gets 
high enough, the top load emits a spark or streamer. Continued driving at the new resonant 
frequency will make the streamer appear thicker, richer, and whiter, but will not increase its 
length. The spark length is proportional to the square root of the input power at the time of 
spark initiation. 


The classical Tesla coil that operates with a capacitor charged to 10 kV or more that 
dumps its energy into a few turns of a primary coil may be able to supply 1000 A or more 
until spark initiation. The order of magnitude of power input is then (10 kV)(1000A) or 10 
MW. This puts out a much longer spark than my solid state system, which is capable of input 
power in the range of 20 kW. Matching performance of the classical system with a solid state 
system would be very expensive, if even technically possible, so solid state driven coils will 
not have spark lengths significantly longer than 2 or 3 feet. 
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Figure 7.1: Tesla Coil Power Supply 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 7—Tesla Coil Driver 7-3 


Readily available solid state switches (IGBTs) are rated at 1200 V. In an H-bridge config- 
uration, they can then supply a square wave of voltage at +600 V. It is prudent to not operate 
above perhaps +500 V. This fact determines the ratings of the components in the top portion 
of Fig. 1. But a square wave voltage of +500 V is not adequate to reliably produce sparks. 
To step up the voltage I went to a system of pulse transformers with primaries connected in 
parallel and secondaries connected in series. These pulse transformers were hand wound on 
plastic coil forms which would accept shaped ferrite cores to improve the magnetic circuit. 
The input had 17 turns and the output had 18 turns of 14 gauge magnet wire, so the trans- 
formers are basically unity isolation transformers. They are functional to perhaps 1 MHz. 
The outputs of four such transformers in series then allows u%:. to be as high as +2000 V. A 
voltage of +1000 V is adequate to produce sparks of up to 24 inches long, which is probably 
enough for a low ceiling. 


As current 24. into the Tesla coil base builds up, the voltage v,. will droop below that 
predicted for ideal unity isolation transformers due to resistance losses, eddy current losses, 
and hysteresis losses. Instead of 4 times the voltage on C4 and Cs, it may be as low as 
3.5 times this voltage. This effect is difficult to precisely quantify, so will be ignored in our 
discussion of ratings of the power supply circuit. The Tesla coil input voltage and current are 
measured directly so the coil input power can be calculated without any need to know droop 
or losses in the power supply circuit. 


On the other hand, the combination of a square wave voltage and a sinusoidal current is 
different from what we are used to, so we need to go back to circuit theory to make sure we 
have all the correct multiplying factors. 


A voltage square wave of value +V, has the same heating capability when applied to a 
non-inductive resistor as a dc voltage V,, hence is said to have an rms value of Vac = V,. The 
corresponding dc current in a resistor would have an rms value of I,, = Ip. The average power 
delivered to the resistor is the product of the rms voltage and the rms current, 


Pave = Vol = Vaclac (7.1) 


Suppose now we apply a sinusoidal voltage V,sinwt to a non-inductive resistor. The 
resulting current is I, sinwt. The average power is 


Pave = - | Vly sin? 6dé = eA = Volp = Vaclac (7.2) 
t Jo V2 V2 2 


When the voltage is square wave and the current is in phase with the voltage but is 
sinusoidal, the integral for average power becomes 


1 pt 2 2 
Po : Vplp sin 848 = =Vply = =Vp(V2Tac) = 0.9Vaclac (7.3) 
TT JO 
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For this case, the average power is no longer the simple product of rms voltage and rms 
current (as for dc and single frequency sinusoids), but has a 0.9 multiplying factor. The 
difference is due to the fact that the square wave voltage is composed of an infinite series of 
harmonics (fundamental, third, fifth, etc.). Each harmonic contributes to the rms value of the 
square wave. The current has no harmonics, so the higher voltage harmonics do not produce 
any contribution to the average power. 


Actually the product 0.9V,-lgc is the apparent power S for these two waveforms, which 
happens to be equal to the average power when the waves are in phase. If there is a phase 
shift, Paye will be less than S by a factor cos @, where @ is the phase angle between the two 
wave. I can readily determine S from the product of rms voltage and rms current (and a 
0.9 multiplying constant). However, P,y,e requires either a wattmeter for measurement or an 
additional measurement of phase angle. Building an accurate wattmeter in the hundreds of 
kHz range is a challenge, and was not done for this research. My scope will estimate the 
phase angle between waveforms but the quality of the estimate is not very good when small 
amplitude signals have significant ripple, as is often my situation. Since I always try to operate 
in the in-phase condition, I will approximate or estimate the average power as the apparent 
power. It should be remembered that the actual average power will be slightly less than the 
numbers I calculate. 


If vc = 1000 V and i,. = 20 A, then the estimate for average power is 


Pave = 0.9(1000) (20) = 18000W (7.4) 


The input current 7;. at the time of spark varies with surface conditions on the top load 
(dust and scratches), humidity in the room, air pressure, cosmic rays producing an ionized 
particle at just the right place and time, etc. But 20 A is within the range, perhaps at the 
high end of the range, so is useful for discussing ratings. This current passing through the 
reactances C’3 and L3 of the Tesla coil model produces the top load voltage necessary for 
spark breakout. Increasing the power supply input voltage just means that the 20 A level 
will be reached sooner. There are some time delays with spark initiation, so sometimes a 
greater input voltage will allow some overshoot to occur, such that the spark starts with an 
input current of 25 or even 30 A rather than 20 A. If I remove a spark initiation point (a one 
inch diameter sphere on the surface of the 30 inch sphere), the average current required for 
spark initiation will rise a little more. If the top load would hold off any spark until the input 
current reached 40 A, at a maximum %, of 2000 V, the nominal input average power would 
be 


Pave = 0.9(2000) (40) = 72000W (7.5) 


In Kansas, I estimated that the spark length was 
£, =0.17V/P inches (7.6) 
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If that formula holds in Colorado, an input power of 72 kW would yield a spark length of 
about 45 inches. 


So for this particular Tesla coil, I need a power supply and inverter that will supply a 
current 24. of 20 A, and perhaps up to 30 or 40 A, in short bursts without failure. A current 
of 20 A in the secondary of a unity isolation transformer will require a 20 A current in the 
primary. Each IGBT (Insulated Gate Bipolar Transistor) must carry 20 A for half the time 
(and be off the other half of the time). The Harris HGTG18N120BND IGBTs are nominally 
rated at 1200 V and 18 A. The 1200 V is a rating that we never want to exceed. However, the 
device will switch up to 100 A until losses cause it to overheat. Operating at 20 A and 50% 
duty cycle should not be a problem if the heat sinks are adequate. In a Tesla coil application 
where we operate in single-shot mode or in a low duty cycle, say, 10 ms on and 1 sec off, the 
device should function with rms currents up to 50 or 60 A while operating. 


A current of 20 A in each of four IGBTs means that the +V capacitors are delivering 80 
A half the time, and the —V capacitors are delivering 80 A the other half of the time. I used 
12 gauge wire which is rated at 20 A continuous, but will survive 80 A or more at very low 
duty cycles. The resistance of the 12 gauge wires will act to damp some of the high frequency 
noise in the power supply, so there is no strong incentive to go to larger wires. 


The current i, rises from zero to the value at spark initiation in a time period in the range 
of 0.5 to 2 ms, then drops to a nearly constant value in the range of 3 to 4 A for the remainder 
of the applied voltage waveform. A total time of 10 ms gives a reasonably rich, fat spark. One 
spark per second helps keep interest in the audience. The average current during the pulse is 
something under 10 A. A duty cycle of 0.01 yields an average current to the Tesla coil of less 
than 0.1 A. Multiply by a factor of four, and the average current out (and in) the capacitors 
C; and C2 is less than 0.4 A. I would not expect a heating problem in the variac or the two 
transformers. They are selected for voltage rather than current or power considerations. 


The 240 V variac is rated at 9 A. It is supplied by the building’s single-phase 120/240 
service. We need separate transformers and full-wave bridge rectifiers to Cy and C2 for voltage 
isolation. We also need about 500 V out of the rectifiers when the variac is full on. At no 
load the capacitor voltage is the input rms voltage times 2. We could use 240:240 isolation 
transformers where 240./2 = 339 V, which is too low, or 240:480 transformers, which yield 
a capacitor voltage that is too high. I needed 240:360 transformers where 3602 = 509 V. I 
was unable to quickly identify a commercial source for such transformers, so I had a friend, 


Floyd Energren, fabricate the two transformers for me, at a cost of $390.00. They are rated 
at 1 kVA each. 


The inductors Ly, and L2 were made from Amidon T520-26 toroids wound with nominally 
200 turns of 18 gauge magnet wire. Two of these toroids were paralleled to form L, and two 
more for Ly. One toroid had a measured inductance of 5.1 mH with an ELC-120 inductance 
meter, so two in parallel should have an inductance of 2.55 mH. I evidently measured the 
parallel combination as closer to 2.4 mH. Actually the permeability increases substantially 
with current flow for these toroids. Where I measured 5.1 mH at no current, I measured 13.5 
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mH at 2 A. The inductance decreases with both higher and lower currents. 


The purpose of these inductors is to smooth out the current flow into the capacitors, 
making the peak current smaller and time of current flow greater. See Mohan for a detailed 
discussion. 


The bridge rectifiers are rated at 1200 V and 35 A. 


Capacitors C, and C2 were each formed of 10 electrolytic capacitors rated at 1400 pF 
and 450 V, connected in parallel. The nominal capacitance is then 1400(10) or 14 mF, but 
the measured capacitance is closer to 15 mF. I would guess that the capacitors would survive 
short periods of operation at 500 V. 


Resistors R; and R2 are rated at 100 W and 2 kf. They are hollow cylinders, mounted 
vertically so air can flow readily both inside and outside. They are connected to the capacitors 
by a Normally Closed contact which is opened when 240 VAC is applied to the variac. The 
initial power dissipated when a 2 kQ) resistor is connected across 500 V is 125 W. The RC 
time constant is 30 seconds, so at worst case conditions the resistors are slightly overloaded 
for only a few seconds. One should allow at least 2 to 3 minutes after shutoff before working 
on the capacitors or associated circuits. 


There are two 500 VDC analog meters at the points marked by the circles with +V and 
—V inside them. These are not used for data collection, but serve to show that things are 
working. 


Inductors Ly and Ls represent the stray inductance of the wires connecting the components 
just discussed (located outside the Faraday cage) to the components in the lower portion of 
the figure, which are located inside the Faraday cage. If tests show that high frequency ripple 
currents are excessive on C, and C4, ferrite toroids can be installed over the wires, to provide 
a little isolation. 


Capacitors Cy and Cs are nominally 72 uF, formed of four banks of twelve 1.5 wF capacitors 
rated at 400 V, located immediately adjacent to each IGBT. They have a higher operating 
frequency than the electrolytics used for C, and C2, so provide some voltage support when 
the IGBTs are switching. 


Variable inductors Lg to Lg have a range of 17 to 180 wH. The typical setting is less 
than 40 uH in operation, so a somewhat different construction from what I used might be 
advantageous. I wound 58 turns of 12 gauge copper wire (not magnet wire) on a length of 
0.5 inch PVC plastic pipe. The resultant winding was about 9.5 inches long by about 1 inch 
in diameter. I then mounted three Amidon R61-050-300 ferrite rods inside a smaller PVC 
pipe that would slide inside the 0.5 inch pipe. The ferrite rods are 0.5 inch in diameter and 3 
inches long. The purpose of these inductors will be explained in the next section. 


Inductor Ly19 represents the stray inductance of the isolation transformers, while Cg rep- 
resents the static capacitance of the feedline to the Tesla coil as measured by a capacitance 
meter. Rapid changes in voltage cause ringing at the 1, node. The frequency of this oscillation 
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was used to estimate Ly. 


The current to the Tesla coil is measured as the voltage across R4, a nominal 0.02 2 resistor, 
but now measuring 0.024 2 after years of operation. Ry, is fabricated from 11 surface mount 
resistors, 0.22 Q and 2 W each. They are mounted in parallel across a gap in two sections of 
2 inch wide copper foil. Where possible, the ground path is made of wide conductors to lower 
the stray inductance of the circuit. 


If the rms current just before discharge would reach 40 A, the power dissipated in Ry 
would be 


P = I?R = (40)?(0.024) = 38.4 W (7.7) 


which is above the rating of a nominal 22 W resistor. Again, the hope is that the resistor will 
survive if the duty cycle is low. If the IGBTs fail in a shorted state, then the stored energy 
in the capacitor bank is available to turn the 0.22 2 resistors, traces on the printed circuit 
board, and anything else in the way, into fuses. An IGBT that fails shorted will eventually 
be blown apart and become an open circuit. A metal cover over the inverter is important at 
high power testing to keep shrapnel from flying around the room. Hearing protection is not a 
bad idea either. 


In addition to this means of measuring Tesla coil current (with R4), we need a voltage 
waveform at the controller that looks just like the current waveform, except cleaner and with 
sharper transitions. That is, we want to convert the sine wave of current into the resonant 
Tesla coil into a square wave. At an rms current of 40 A, the voltage across the nominal 0.02 
Q resistor has a peak value on the order of 100 mV, far too low to clean up with diodes. We 
really need a square wave with amplitude of a volt or two. We can get such a wave with the 
use of a current transformer, shown just to the left of R4. 


The current transformer is fabricated with two stacked ferrite toroids, what I believe to 
be the Phillips TX22/14/6.4-3E2A. Ferrites typically arrive with no markings whatsoever, so 
there is no simple method of determining the exact material used by looking at the device. 
In this case, I believe I received a ferrite material different from the type I ordered. But the 
materials were similar enough that the toroids received worked satisfactorily. The primary 
of the current transformer is a straight section of 14 gauge copper wire connecting Ry, to the 
bottom isolation transformer, passing through the center of the two stacked toroids. The 
secondary is 25 turns of 20 gauge magnet wire wound over the two toroids. If the two ends 
of this winding were shorted (and if everything were properly designed), a current of 25 A to 
the Tesla coil should result in a 1 A current in this winding. 


A current transformer is always operated with a short across the secondary for two reasons. 
One is linearity, which we are not particularly interested in here. The other is that if the 
secondary is open circuited, the secondary voltage can rise to extremely high values, enough 
to damage the insulation in current transformers built for 60 Hz service. The second reason 
is used to advantage here. The secondary is loaded with two chains of two diodes each (Dy, 
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Dz and Ds, D4), the chains pointed in opposite directions. As the Tesla coil current passes 
through zero, the secondary voltage also goes through zero but very rapidly. Two diodes in 
series have a very high impedance when the voltage across each diode is less than 0.5 V or so. 
Then while the Tesla coil current is still near zero, the secondary voltage will be clamped at 
two diode drops, say 1.4 V. The secondary voltage is not a perfect flat top square wave, but 
will reach 1 V for a Tesla coil current of hundreds of mA, and then rise slowly to 1.4 V ora 
little more for Tesla coil currents of tens of A. This approach works surprisingly well for the 
purpose. 


7.2 Gate Driver and IGBTs 


The gate driver for one pair of the IGBTs is shown in Fig. 7.2. This schematic represents the 
components on one PC board. There are four of these boards in Fig. 7.1. The top trace is 
connected to +V and to the drain of IGBT @Q2. The bottom trace is connected to —V and 
to the source of IGBT @Q1. Inputs T’P2 and T’P3 are connected to the corresponding outputs 
of the controller board, to be discussed in the next section. The junction of the two IGBTs, 
TP1, is connected to the unity isolation transformer. 


The MIC4420 (or TSC4420) is a non-inverting driver for gates of MOSFETs and IGBTs. 
It can supply up to 6 A while charging a capacitive load. It turns on with an input of 2.4 V. 
In some of my early development, I suspected that electrical noise on the wiring was causing 
false gate drive signals. As previously mentioned, if gate signals are present at both IGBTs 
in the H-bridge at the same time, we have a short across a rather large capacitor bank which 
very quickly causes the IGBTs to explode. As a patch, I put the voltage divider RC circuit 
shown on the left of the driver. For a nominal input voltage of 15 V at TP2 or TP3, the 
nominal voltage at the input of the 4420 is 5 V. The capacitors add some pulse shaping and 
filtering to the waveform. A voltage transient at the inputs must exceed 7.2 V for a long 
enough time period to charge C3 (or Cg) to 2.4 V before the 4420 turns on. The divider may 
not be necessary in the present version, but it works nicely, and I have not felt led to remove 
it. The two gate driver circuits are identical, so I will discuss only the top circuit. 


Each IGBT gate is charged and discharged through a pulse transformer. I used a ferrite 
toroid (probably the Philips TX22/6.4/6.4-3E2A, which I also used for the current transformer 
of the previous section) with 5 turns of 20 ga. magnet wire for the primary and 10 turns for 
the secondary. 


The numbers 22, 14, and 6.4 in the part number refer to the outside diameter in mm, 
the inside diameter, and the height or thickness of the toroid, respectively. The windings 
need to be physically separated enough to withstand up to 500 V, so I put the primary on 
one side of the toroid and the secondary on the other side. The leakage inductance would be 
smaller if the windings are on top on one another, but then I would need to place additional 
insulation between windings. This can be a challenge in relatively small toroids, so I opted 
for air insulation between totally separate windings. 
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Figure 7.2: Tesla Coil Gate Driver 
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The IGBT gate needs at least 10 V to turn on, and 12 to 15 V is better, to assure some 
noise immunity. The gate voltage should not exceed +20 V on a repetitive basis, and the 
absolute maximum is +380 V. When the first pulse is applied, the capacitor C7 is discharged, 
and its voltage does not change instantaneously. If 15 V is applied to the left capacitor 
terminal, the right side of the capacitor will also try to have 15 V on it. However, the two 
back-to-back 1N4737 zeners will clamp the voltage at about 8 V. There will then be 7 V across 
the impedance of the 4420, the resistance of circuit traces, and C7. The pulse transformer 
will put out about double the input voltage or 16 V. If there are transients or oscillations, the 
bidirection zener, P6KE18, clamps the gate voltage at 18 V. 


While current is flowing out the pulse transformer secondary to the IGBT gate and the 
bidirectional zener, it is also flowing into C7. With a pulse train of approximately 50% duty 
cycle, C7 will charge to 7.5 V after the first few cycles. Once C7 is charged, voltage excursions 
of +15 to 0 at the 4420 output will be seen as a pulse train of +7.5 V at the pulse transformer 
input. The gate would then see +15 V. The bidirectional zener and the 1N4737 zeners will 
not conduct after the first few cycles. 


There are many tradeoffs in the gate driver circuit. C7 and the leakage inductance of the 
pulse transformer form a resonant circuit which is excited by the first pulse. As energy trades 
back and forth between C7 and the inductance, the gate voltage can drop below the threshold 
value, causing the IGBT to try to turn off when it should be on. This effect seems to be helped 
by keeping the leakage inductance small, that is, by having as few turns as possible on the 
pulse transformer. On the other hand, the time that a pulse transformer can hold a voltage 
is proportional to the number of turns, so for low frequency operation, one needs many turns. 


One should not use ceramic capacitors for C7. Some of them change capacitance value 
with temperature, and the initial pulse will heat the capacitor enough to cause a thermal 
resonance in the gate voltage in the 3 to 10 kHz range. A polyethylene capacitor rated at 50 
V is much better than a ceramic capacitor. 


If it were not for the initial transient, design of the driver would be very straightforward. 
I could not find a book or paper which discussed the design of this particular configuration. 
Modeling would be difficult because the current and voltage output of the 4420 under the 
overload conditions of the first pulse would be strongly dependent on the exact construction 
of the printed circuit board. The leakage inductance of the pulse transformer would vary with 
the specific toroid and also with the hand winding technique. I did tests with various values 
of C7 and R7, and basically quit when I got something to work. There could easily be much 
better designs than this one. 


Circuit board layout of the drivers and IGBTs is also critical to success. The system needs 
to be as compact as possible, consistent with adequate heat sinking and voltage separations. 
The three dimensional aspects are difficult to present in two dimensional drawings, so I will 
attempt a verbal description. The IGBTs are mounted on separate heat sinks and placed on 
the board where the ground foil has been removed. The heat sinks are finned black aluminum, 
approximately 1 by 1.6 by 2 inches outside dimensions. There are two pins about an inch apart 
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that go through the board along with the IGBT pins to provide some mechanical stability. 


The IGBTs are mounted on their heat sinks without an insulating pad. This maximizes 
the heat transfer, but with the disadvantage that the heat sinks are electrically charged. I 
think one set of IGBTs was destroyed due to a discharge from a heat sink to the grounded 
case. I now set the driver/IGBT board, oriented with heat sinks down, on a 1/2 inch thick 
layer of polyethylene. 


Next to the MOSFET or IGBT symbol in Fig. 7.2 is an upward directed diode. There is 
such a diode built into these particular IGBTs (the Harris HGTG18N120BND), but I added 
these MUR4100 Fast Recovery diodes in parallel to the internal diodes. The external diodes 
should not be needed, but I thought they helped a little with wave shaping. These diodes 
(internal or internal plus external) are an important part of the circuit. They allow for current 
to continue flowing after one IGBT turns off and before the other one turns on. They also 
allow for energy stored in the Tesla coil to flow back into the capacitors when the controller 
is turned off and all gate pulses cease. The waveforms of voltage and current look almost the 
same after the controller is turned off, except for a 180° phase shift in voltage. While the gate 
pulses are applied, power flows from IGBTs to Tesla coil. When pulses are removed, power 
flows from coil to IGBTs. This is most obvious when the operating voltage is too low for a 
discharge to occur. The Tesla coil is then fully energized when the end of the gate pulses 
occurs. 


As mentioned, this inverter circuit provides a nominal square wave voltage to the Tesla 
coil input. A square wave can always be composed into a fundamental and a series of odd 
harmonics of sine waves. The Tesla coil is not resonant at exactly three (or five, or seven) 
times the fundamental frequency, so the harmonics always face a very high surge impedance. 
The current will build up at the resonant frequency but not at the harmonics. This means 
that a square wave of applied voltage will produce only a sine wave of current. This sine 
wave will be in phase with the voltage at resonance, will lag above resonance, and lead below 
resonance. 


IGBT switching at resonance means the IGBTs turn on and off when the current through 
them is small. This is a desirable feature, in that it reduces the losses in the IGBTs. It cannot 
be too small, however. 


It is often not emphasized in power electronics books, but power MOSFETs and IGBTs in 
this circuit configuration essentially require an inductive load. With a lagging current, when 
one IGBT is turned off, current will continue to flow, charging and discharging the drain-to- 
source capacitances (and capacitors Cg and Co if installed), until the voltage across the other 
IGBT has been reduced to zero. If the voltage reaches zero while the devices are off, current 
will start to flow through the diode(s) of the other IGBT. When a gate pulse arrives at this 
IGBT, it already has zero voltage across it, so the turn-on process works very smoothly. On 
the other hand, a leading current will result in each IGBT turning on with full voltage across 
it. Very large transient currents will flow, with steep rise times. Stray inductances will cause 
substantial ringing in voltage due to these transient currents. A leading current is a disaster 
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on its way to happen in an inverter circuit like this. 


We need to be able to tune the Tesla coil across resonance. One way to assure the IGBTs 
are always operating into a lagging load, even if the Tesla coil happens to be drawing a leading 
current, is to add additional inductance in parallel with the Tesla coil. This is accomplished 
with inductors Lg through Lg in Fig. 7.1. This inductor current also provides for relatively 
soft switching of the IGBTs when the Tesla coil current is very small. 


The difference between hard and soft switching is shown in Fig. 7.3. We have two plots of 
an (approximately) square wave of voltage applied to a Tesla coil, with the resulting (approx- 
imately) sinusoidal current. In the Dec 19 plot, the voltage takes about 300 ns to make the 
transition. There is little ripple and the current waveform is reasonably smooth. In the Sep 
30 plot, however, the voltage transition lasts only about 100 ns. This means that the voltage 
waveform contains more high frequency components, which can excite various resonances in 
the Tesla coil itself, the transmission line, and the IGBT driver. There is some ripple on 
the voltage and more on the current. The Dec 19 plot is obviously the preferred method of 
operation. 


There are two necessary conditions to get soft switching. The current flowing through the 
conducting IGBT must be sufficient at the time of turn-off, and the dead band (the time lapse 
before the other IGBT turns on) must be adequate. The controller chip (discussed later) can 
have a dead band as low as 80 ns, but is set for about 400 ns in the Dec 19 plot. 


Note that the peak inductor current is the critical factor, not the Tesla coil impedance. 
At higher frequencies, the peak current in a given inductor will decrease, requiring the use of 
a smaller inductor. If 40 wH works at 150 kHz, then we would expect a 300 kHz coil to need 
a 20 wH inductor. 


7.3. Tesla Coil Controller 


The controller supplies the gate drive signals to the gate driver discussed in the previous 
section. It is adjusted so the initial frequency is near resonance. It senses the current and 
strives to maintain the proper phase relationship between voltage and current by relatively 
small changes in frequency of the Motorola 34066 controller chip. The circuit is shown in 
Fig. 7.4. 
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Figure 7.3: Soft switching (Dec 19), harder switching (Sep 30). 
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Figure 7.4: Tesla Coil Controller 
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The frequency of the 34066 (IC1) is determined by the current into pin 3. This current is 
supplied by a LM324N single supply op amp (IC2) in a voltage follower mode. The voltage 
of the 34066 pin 3 must be at least two diode drops and must be less than 6 V. The source 
voltage at pin 3 of IC2 is obtained by voltage division of R1, D2, R2, and either RV1 or RV2 
depending on which solid state switch (IC3A or IC3B) is closed. D2 is a 1N4734 zener diode, 
with a nominal voltage of 5.6 V. This 5.6 V is further reduced by the voltage division between 
R2 and RV1 or RV2. The resulting current into the 34066 will not be able to produce output 
frequencies above 250 or 300 kHz, but the Tesla coils of interest are resonant at frequencies 
below 200 kHz, so this frequency range is quite adequate. 


We never want both legs of the H-bridge on at the same time, and, because of various 
lag times involved in turn-on and turn-off of the IGBTs, we actually want some dead time 
between when pin 14 of the 34066 turns off and pin 12 turns on. This dead time is controlled 
by R4 between pins 1 and 2 of the 34066. A value of 430 2 for R4 gives a dead band on the 
order of 400 ns. 


This version of the controller has Enable (pin 9 of the 34066) always held high, so the 
controller wakes up sending gate drive pulses out pins 12 and 14. These pulses appear at solid 
state switches IC3C and IC3D, which are held open until the controller has stabilized and a 
request is made for the gate drive pulses to actually leave the controller. The board is capable 
of three different operating modes: 


1. Continuous 
2. One Shot 


3. Repetitive Pulse 


There are three manual switches that control the operating mode. Switch SW1 is a 
momentary contact toggle switch that initiates each One Shot, and also initiates Continuous 
as well as Repetitive mode. Switch SW2 is up(on) for either One Shot or Repetitive mode and 
down(off) for Continuous mode. Switch SW3 is up(on) for Repetitive mode and down(off) 
for One Shot or Continuous mode. For One Shot operation, SW2 is up, SW3 is down, and 
SW1 is toggled momentarily. Both switches SW2 and SW3 are up for Repetitive and down 
for Continuous mode. 


The Continuous mode is used for initial low power testing. The variac driving the main 
power supply is set at 5 to 10% of rating and the control voltage to switches IC3C and IC3D 
remains high once set high. Frequency is controlled by pot RV1. There is no feedback to 
control the phase between voltage and current at the input to the Tesla coil, so the phase 
difference may drift with time. We would use this mode to determine the resonant frequency 
for a new coil, or a coil with a different top load, and also to adjust the inductors on the 
IGBTs for minimum ringing. It is useful to do basic circuit tests, such as continuity and the 
functionality of all components. 
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The One Shot mode is initiated by switch SW1. This Sets the flip flop IC4A so frequency 
control of the 34066 is passed from RV1 to RV2 for a fixed number of gate drive pulses. The 
number is determined by the particular output connection of the 4040 counter IC5. The 
schematic shows pin 12 of IC5 connected back to Reset of IC4A through switch SW2, which 
provides a Reset pulse to IC4A after 256 pulses into IC5 from the 34066 chip. This can be 
changed with a soldering iron as desired. Using pin 14 of IC5 would result in 512 pulses 
applied to the Tesla coil, etc. Current into the Tesla coil builds with each pulse until a spark 
is emitted (if the variac is set high enough), after which the current drops and remains more 
or less constant. The resonant frequency drops slightly because of the space charge in the 
spark. Applying the same Tesla coil input current for a longer period of time makes the spark 
look fuller and richer. The spark length remains the same. A sequence of 128 pulses is likely 
to result in a thin, wimpy, violet colored spark, for example. Once the system is demonstrated 
to be stable and robust, a sequence even longer than 128 or 256 will probably be desirable. 


The Repetitive Pulse mode is used to get multiple sparks from the Tesla coil at some fixed 
interval up to a few seconds. Switch SW3 must be closed, so the output of the second counter, 
IC6, is fed back to the Clock of IC4A. The D input (pin 5 of IC4A) is always high, so a clock 
signal acts to set the Q output (pin 1) high, the same as the Set pulse applied by switch 
SW1. Assume a frequency of 143 kHz, with a period of 7 ws. A pulse train of 128 pulses 
would require (7)(128) = 896s = 0.896 ms. The maximum delay between sparks is when 
both counters use Q12 as the output. Q12 changes state for every 2'! = 2024 input pulses, 
or 14.168 ms for the first counter. The second counter then changes state at 2024(14.168) = 
28676 ms = 28.7 s. This is somewhat too long for entertainment (everyone gets bored), so 
something less would be appropriate. The adjustment is by a factor of 2”, so the possible gaps 
in seconds between sparks are 28.676, 14.338, 7.169, 3.584, 1.792, 0.896, etc. The schematic 
shows Q12 used by IC5, and Q7 by IC6, with an interval between sparks of about a second. 


We will now discuss the One Shot mode in more detail. Switch SW1 Sets IC4A. IC4A pin 
1 Resets the two counters IC5 and IC6 and applies a logic high to the D input of IC4B. The 
next pulse from the 34066 pin 14 Clocks IC4B pin 13 high, which turns on the switches IC3C 
and IC3D. The 4066 switch IC3A is still on, so pot RV1 is still controlling frequency for the 
first few pulses to the Tesla coil. IC4B also applies a logic high to the D input of IC9B. One of 
the lower outputs of counter IC5, Q4 shown in the schematic, then Clocks IC9B, which turns 
off switch IC3A. IC9B applies a logic high to the D input of IC9A, which is then Clocked by 
the next pin 12 pulse from the 34066. When pin 1 of IC9A goes high, switch IC3B is turned 
on and frequency control is passed to RV2. 


During the time when neither RV1 nor RV2 are in the circuit, the capacitors C3 and C12 
charge through R2. The frequency of the 34066 skews upward at the rate of about 2 kHz 
per cycle. This version of the controller causes RV1 to turn off with a pin 14 pulse from the 
34066, and RV2 to turn on with the next pin 12 pulse from the 34066, or one-half cycle. If 
this slew rate is too large, either C3 or C12 will need to be increased in size. 


The RC circuit C16 and R5 is a patch to make sure that IC9B wakes up in the proper 
state. 
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As mentioned earlier, we monitor the Tesla coil current with a small current transformer 
terminated with two pairs of diodes. This current transformer terminated by diodes acts more 
like a comparator in that a few mV swing on the input would result in maximum swing on 
the output. It presents a fairly good square wave to port TP1, starting with the first pulse to 
the Tesla coil, many pulses before any current is visible on a scope. 


The next step is to convert this small analog signal (plus and minus two diode drops) to 
a clean digital signal, +15 to zero, a nontrivial task. The solution is clever, even elegant. I 
wish I could say it was original, but I saw it in an application note. Two pins of the 4046 
phase-lock loop are used, one for input and the other for output. There is no feedback. The 
4046 converts the analog into a digital signal very nicely. It has an internal biasing network 
so it will work if a capacitor C8 is used to allow a dc offset. It was my experience that the 
duty cycle of the resulting square wave was not close enough to 50%, so I added an external 
bias network R10, RV3, and R1l. RV3 needs to be adjusted so the output at pin 2 of the 
4046 is near a 50% duty cycle. This is only done once, in initial testing of the circuit. 


The current square wave is sent through four stages of a 4050N non inverting buffer, to 
provide a little extra delay. This delay allows one to tune through resonance a greater amount. 
The waveform is then sent to the Reset of IC9A. The other two stages of the 4050N, IC7E 
and IC7F, are not used. 


Controller frequency is controlled in the following manner. Pin 12 of the 34066 turns on 
half the IGBTs, initiating one edge of the voltage waveform applied to the base of the Tesla 
coil. Phase difference between voltage and current can be measured with respect to this edge. 
This leading edge also Clocks pin 1 of IC9A high which turns on switch IC3B. RV2 starts to 
discharge C3 and C12, reducing the voltage at pin 3 of IC2A, and lowering the frequency of 
the 34066. Then the leading edge of the current waveform comes to Reset of IC9A, which 
turns off switch IC3B. RV2 is out of the picture and frequency starts to drift back up. When 
RV2 is properly adjusted, the frequency will be drifting back and forth across the resonant 
frequency. 


At resonance, the duty cycle on pin 1 of IC9A should be on the order of 0.3 to 0.5. As 
the frequency drifts below resonance, the duty cycle decreases, RV2 is not connected as large 
a fraction of the time, the voltage on C3 and C12 increases, increasing the current to pin 3 of 
the 34066 and thus increasing the frequency. The opposite happens when the frequency drifts 
above resonance. If the duty cycle at resonance is 0.4, things need to work over a duty cycle 
range of say 0.2 to 0.6. If we start off with a resonance duty cycle of 0.1 or 0.9, we just do 
not have the tuning range that we need. 


The voltage on C3 and C12 is buffered through IC2A, a LM324 op amp in a voltage 
follower configuration. This is an old, slow op amp, but great speed is not needed in this 
application. A precision op amp would probably be better, but a quick look failed to turn up 
a single supply op amp that met all the other requirements. I did not want to add a —15 V 
supply to the printed circuit board just for the op amp if this single supply version would do 
the job, and it does reasonably well. 
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The 34066 has a double frequency ripple at pin 1, one hump controlling gate pulse length 
at pin 12 and the following hump the length at pin 14. Our technique of frequency control 
results in a single frequency ripple on pin 3 of the 34066, which will cause one of the gate pulses 
to be longer than the other by a few percent. The Tesla coil sees a slightly non symmetrical 
waveform with the time spent at V+ somewhat different from the time at V—. This is not 
a particular problem to the Tesla coil, but the inductors Lg to Lg in Fig. 7.1 will integrate 
this non symmetry over the duration of the discharge. Instead of an inductor current varying 
symmetrically from +1 to —1 A, it might vary from +0.6 to —1.4 A. We might have not 
enough current to properly switch the IGBTs at one extreme, and too much at the other 
extreme. The current waveform might show more noise on the positive going edge than the 
negative going edge, or vice versa. 


A value of 4.3 MQ for R12 was found empirically to balance the current into the current 
mirror of the 34066 and give equal length output pulses. 


Another interesting feature of the inductor currents is the following. In steady state 
operation, the IGBTs switch when the inductor current is at a peak. The voltage across 
each variable inductor is approximately a square wave, so the inductor current is basically a 
saw tooth wave. This saw tooth should be symmetric about zero. The current for the first 
pulse starts at zero rather than a negative peak, so if the first pulse is the same length as 
the following ones, the inductor current will rise to about double the steady-state peak. The 
current waveform will then drift down with time as the losses in the circuit cause a change 
toward symmetry. The IGBT switching current is too high for half the IGBTs and too low for 
the other half during the first few cycles. This can easily cause ringing and IGBT destruction. 
The current waveform is helped immensely by shortening the first pulse until the peak inductor 
current is about equal to the long term peak (rather than double). It is not essential that the 
peaks be exactly those of steady state during the first few cycles. Peak values between 70% 
and 130% of the long term peaks seem to work in an acceptable fashion. 


The IGBT gate driver circuit has capacitors that will be charged up during the first few 
pulses. The output of the pulse transformer in Fig. 7.2 will be somewhat longer when C7 is 
discharged than when it is charged to approximately half the supply voltage. If this additional 
duration is greater than the dead band of the 34066, we can have both IGBTs on at the same 
time. So this is a second reason to have the 34066 generate shorter pulses (greater effective 
deadband) for the first few pulses of each discharge. 


Frequency and pulse width output of the 34066 are controlled by external components 
attached to pins 1 and 2 (oscillator) or pin 16 (one shot). The circuit in Fig. 7.4 has component 
values selected so the oscillator input does the controlling and the one shot input just floats 
along, do basically nothing, except for the first few pulses. The components C7, R8, D6, RV4, 
and C20 perform this function of shortening the first few pulses. The capacitor C7 is charged 
by internal transistors of the 34066 to 5.1 V, starting the one shot period. The one shot 
period ends when the capacitor is discharged to 3.6 V. R8 is selected so the discharge period 
is greater than for the components attached to pins 1 and 2 of the 34066, so the oscillator 
is normally in control. But when SW1 is pressed to initiate a pulse train, pin 12 of IC4B 
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goes low. This pulls the top plate of C20 low (actually below ground, but the excursion is 
limited to one diode drop by the zener D6). C7 now discharges through both R8 and RV4, 
so the pulse width is controlled by this circuit until C20 gets charged to a new voltage. The 
duration is basically five or six pulses if the first pulse in the pulse train is adjusted to be half 
the length of the pulses controlled by the oscillator. 


The 34066 is able to drive IGBT gates directly, and certainly the four inputs of the gate 
driver circuit of Fig. 7.2, but the 4066 switches have too great a series impedance to drive 
the capacitance of the lines and these four inputs in parallel effectively. Therefore I added a 
IR2110 gate driver, shown on the right side of Fig. 7.4. The pulse train outputs go to TP2 
and TP3, on a terminal strip mounted on the controller box. This chip allows the use of two 
different power supplies for input and output and I used that feature. There will be significant 
current flow in the power supply connected to the gate drivers, which causes some switching 
noise and some non trivial voltage dips at the leading edge of each pulse. This could feed 
back into some of the voltage sensitive circuits of the remainder of the controller, such as the 
voltage input to the input of the LM324. This sort of problem is difficult to trouble shoot, 
so it is simpler to just have two power supplies. The controller board runs on a 1 A, 15 V 
supply, while the gate drivers work on a 3 A, 15 V supply. 


7.4 Waveform Plots 


Fig. 7.5 shows the shortening of the first few pulses to the gate driver. These waveforms were 
measured at TP2 and TP3 of Fig. 7.4 with the X10 probes of the HP54645D oscilloscope. 
The first pulse is 1.55 us in length, the second 1.62 ys, the third 2.00 ps, etc. By pulse 7, we 
are basically at the final value of 2.80 us. 


Fig. 7.6 shows the Tesla coil input voltage and current during the first 30 cycles. The 
Al trace is the primary voltage of the isolation transformer in Module 1, measured with a 
1:200 probe. It shows a peak value of about 0.5 V, so the actual primary voltage is 200 times 
this value, or about 100 V. The actual measured voltage on the power supply was +107 V. 
The Tesla coil input voltage would be four times this value because there are four isolation 
transformers with the secondaries in series, or about 400 V. This is less than the breakout 
voltage, so no streamer occurs. 


The A2 trace is the current. It looks like noise for the first few pulses. DO is the control 
voltage for RV2 while D1 is the control voltage for RV1. RV1 needs to be adjusted so that the 
pulse frequency is close to resonance, which produces a recognizable current waveform before 
control is passed to RV2. Once a small current is establised at resonant frequency, the circuit 
in Fig. 7.4 is able to switch RV2 in and out of the circuit at the necessary rate to maintain 
resonance. 


Voltage and current for 4.5 ms are shown in Fig. 7.7. The applied voltage is the same 
as Fig. 7.6. The current continues to build up for the entire 3.2 ms that gate signals are 
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Figure 7.5: Shorter Pulses at Start 


applied, but at a progressively slower rate toward the end. This current was not adequate to 
produce a streamer. If one looks closely, the voltage is decreasing during this 3.2 ms period, 
then increases when the gate signals are removed. At 0.1 ms, the rms voltage was 96.6 V, 
decreasing to 89.0 V at 3.18 ms, and increasing to 100.6 V at 3.22 ms. The decrease is due to 
the discharge of the power supply capacitors. As discussed earlier, the transformers and other 
components in the power supply were sized more for average power rather than peak power. 
A voltage decrease of less than 8% (for this example) is just not a problem. 


The increase of applied voltage when the gate drive is removed requires a bit more ex- 
planation. The resonant circuit of the Tesla coil has become almost fully energized (for this 
particular input voltage) at 3.2 ms. This energy is swapping back and forth between electric 
and magnetic forms. The Tesla coil is physically large with significant inductance. The de- 
scription may not be precise, but it appears that the current has a large inertia, just like a 
very large and heavy pendulum. The current really wants to continue flowing! If an inductor 
carrying a large current is completely disconnected from its source, the voltage at its terminals 
will rise to extreme values, often enough to spark across the terminals and dissipate the stored 
energy. In this case, the IGBTs are still connected. Current can flow through the built-in 
diodes and recharge the capacitor bank. The capacitors basically limit the voltage. All the 
voltage drops are reversed in sign because the capacitor bank is being charged rather than 
discharged. This causes the measured voltage across this transformer winding to increase. 
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Figure 7.6: Voltage and current at startup. 


Voltage and current waveforms around the time of the end of the gate signals are shown 
in Fig. 7.8. The left portion of the figure shows voltage and current in phase, as one would 
expect for a basically resistive load at resonance. The right side portion shows voltage and 
current out of phase, as one would expect when the Tesla coil becomes a source rather than 
a sink. The current never misses a beat! The voltage waveform has to do the adjusting. The 
transition is best handled if the voltage waveform is low for about half the usual period and 
high for about half the usual period. After this event, lasting for half a cycle of the current 
waveform, the phase relationships are proper, so no additional phase shifting is necessary. 


The lower portion of Fig. 7.8 shows the gate drive waveforms and also the signal that 
terminates the pulse train. DO is signal HO (High Out) of the IR2110. D3 is LO (Low Out). 
D1 is pin 13 of IC4B, which controls the switches IC3C and IC3D that pass the gate signals 
to IR2110. Signal HO is connected to the negative half of the H-bridge of IGBTs, so when it 
is high, the output voltage to the Tesla coil is negative. D1 has a delay relative to the gate 
signals HO and LO because of the because of the ripple effect in the counter IC5. D1 does 
not go low until after HO has gone high for the last pulse. One might expect the last pulse to 
just be a blip because switch IC3D has been turned off just after its output (pin 10) has gone 
high. However, the circuit board trace between pin 10 of IC3D and pin 10 of IC10 does not 
instantly go to zero voltage just because the switch is off. The chips have very high input and 
output impedances, so the voltage on the circuit board trace capacitance takes some time to 
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Figure 7.7: Voltage and current for 4.5 ms, without a streamer. 


bleed off, to get low enough that HO will change state. In my case, this time was excessive, 
so I had to add R18, empirically determined to cause a state change in about a quarter-cycle 
of the current waveform. My power supply is driving the Tesla coil until the end of the short 
HO pulse. At that point the Tesla coil becomes the source and the voltage transitions are 
determined entirely by the current transitions. 


After a period of time, the energy in the Tesla coil has decreased until it is no longer 
able to recharge the power supply capacitors. We now have what we might call an unloaded 
series RLC circuit where we are measuring the voltage across and the current through the 
capacitances of the IGBTs. From Circuit Theory I we recall that the voltage across and the 
current through a capacitor has to be 90 degrees out of phase. This is shown in Fig. 7.9. The 
left portion of the figure shows some power supply recharging occurring and the right side 
portion shows two sine waves out of phase by 90 degrees. 
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Figure 7.8: Voltage and current at end of gate signals. 
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Figure 7.9: Voltage and current after recharging the power supply is finished. 
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STREAMERS AND SPARKS 


Our results for sparks and streamers apply to the case of a vertical coil above a ground 
plane, driven at the base with a more-or-less square wave of voltage from a solid state driver. 
Let me present some definitions found in the IEEE Standard dictionary [1]. 


arc. (1) A discharge of electricity through a gas, normally characterized by a voltage drop 
in the immediate vicinity of the cathode approximately equal to the ionization potential of the 
gas. (2) A continuous luminous discharge of electricity across an insulating medium, usually 
accompanied by the partial volatilization of the electrodes. 


corona. (1)(air) A luminous discharge due to ionization of the air surrounding a conductor 
caused by a voltage gradient exceeding a certain critical value. (2)(gas) A discharge with 
slight luminosity produced in the neighborhood of a conductor, without greatly heating it, 
and limited to the region surrounding the conductor in which the electric field exceeds a 
certain value. 


spark. A brilliantly luminous phenomenon of short duration that characterizes a disruptive 
discharge. Note: A disruptive discharge is the sudden and large increase in current through 
an insulating medium due to the complete failure of the medium under electric stress. 


streamer. An incomplete disruptive discharge in a gaseous or liquid dielectric that does 
not completely bridge the test piece or gap. 


In general usage, arc is for big currents. An arc welder may run hundreds of amps through 
a welding rod to make a weld. It is possible for a power arc to develop across a spark gap in 
a conventional Tesla coil, a highly undesirable event. About the closest to an arc that occurs 
in my solid state drive is when the IGBTs in both legs of the H-bridge are turned on at the 
same time, forming a short across a large capacitor bank. This event produces a very loud 
thump, smoke, and flying pieces of silicon. The traces on the printed circuit board volatilize. 
If not a real arc, it is a good enough imitation for me! 


Corona is a common occurrence around Tesla coils. Leads and connections from the high 
voltage terminals may be in corona, although it may be difficult to see in the light produced by 
the spark gap and discharges from the toroid. It will be especially noticeable on the lead from 
the top of the secondary to the bottom of the extra coil in a magnifier connection. Corona 
is a violet or purple discharge extending out from the conductor a centimeter or two. Power 
losses are not excessive, but using larger diameter conductors to avoid corona will also reduce 
resistance losses and may help performance of the coil. 


It would seem from these definitions that one could use either spark or streamer for the 
discharge from the top of a Tesla coil. There are some differences between discharges to air 
and discharges to a nearby object, so separate nomenclature is appropriate. I think I will 
refer to discharges from a toroid to air as streamers, and discharges to some object as sparks. 
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Those interested in more information about these terms should refer to one of the books on 
lightning by Uman [2], [3]. 


My measurements of sparks and streamers have been in two different locations (Manhattan, 
Kansas and Canon City, Colorado). In Kansas I used a daisy chain of IGBTs in series to get 
the necessary drive voltage, while in Colorado I used four modules, each with a pair of IGBTs 
driving the input of a unity pulse transformer, and the transformer outputs connected in series 
to get the necessary drive voltage. The module approach simplifies the power supply and allows 
for greater input voltage at the base of the Tesla coil by simply adding more modules. It may 
also be more robust than the Kansas approach. Unfortunately, the pulse transformers add 
stray inductance to the Tesla coil input circuit, and droop under load conditions. Breakout 
still occurs from the Tesla coil top load, but the waveforms are just not as pretty. 


The facilities are also much different. In Kansas, the clearance to walls and ceiling was 
much greater. The input impedance for coil 145 was on the order of 25 2, mostly due to ohmic, 
skin, and proximity losses of the coil itself. In Colorado, the input impedance was above 50 
Q until I installed aluminum flashing on the ceiling above the Tesla coil. This lowered the 
impedance to fairly close to the Kansas value, which was low enough to allow adequate current 
to flow for a given voltage to get sparks. I could perhaps get a lower input impedance by 
covering the concrete floor with copper sheeting, but am not sure that the improvement would 
justify the cost and effort. 


Of course, the elevations are different. The Kansas facility was at about 1000 ft elevation 
while the Colorado facility is at about 5600 ft. The higher elevation should make streamers 
easier to initiate and possibly extend their length. The humidity is also much different. The 
Kansas facility had a dirt floor covered with asphalt millings. The soil was quite hydrophilic 
so soil moisture would be pulled up from deep underground and evaporate into the facility. 
Humidity readings of 60% were not uncommon when the outside humidity was at a typical 
Kansas value of 40%. In Colorado, moisture also comes in from the soil through the old 
concrete floor, but the humidity levels are more typically in the range of 20% to 30% inside 
the facility. 


8.1 Waveforms 


I apply a square wave of voltage at the bottom of a Tesla coil, which causes a sinusoidal 
current to flow, as shown in Fig. 8.1. These waveforms are for coil 12T operating at a 
voltage insufficient for a discharge to occur. (The variac was set at 30%). The oscilloscope 
is the HP54645D, which contains the capability to calculate quantities like rms voltage and 
frequency of a waveform, and also to print out the screen image to a printer. 


Voltage is actually the voltage output of a 10:1 voltage divider. The vertical scale for 
channel Al is thus 500 V/div rather than 50. At the lower left of the screen image we see 
Vrms(A1)=44.34 V. The rms value would actually be 443.4 V. The rms of a perfect square 
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Figure 8.1: Square wave of voltage, sine wave of current at base of coil 12T below breakout. 
Kansas drive configuration. 


wave is the same as the peak value, so I am applying a voltage of approximately +443.4 V to 
the coil. 


The current waveform is the voltage across a 20 mQ resistor. The current corresponding 
to a voltage of 293.3 mV is 


Vrc —-393.3mV 
Tea = 


= ———— = 19. A 
R 20mQ ee 


We saw in the previous chapter that when the voltage is square wave and the current is 
in phase with the voltage but is sinusoidal, the integral for average power becomes 


ape 2 2 
Pie | Vplp sin 040 = —Vplp = —Vp(W2lac) = 0.9Vaclac (8.1) 
T JO T T 


At resonance, when the input impedance is real, it can be defined as 


P. 2/Wplac _ 2V2V, V, 
Zin = Rro = 2V Vp loc 2 2V2Vp 215.9 ¥ee (8.2) 
Té. TL. Tae lac 
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In the following, I will refer to the rms voltage applied to the coil as Vro and the rms 
current as Ipc. The average power for Fig. 8.1 is approximately 


P =0.9VrcIrc = (0.9)(443.4) (19.66) = 7846 W 


The input resistance is 


443.4 
Rea O0- =903 
ro = 0970 6 


The frequency is calculated from the time lapse between two adjacent leading edges of the 
voltage at the points indicated by the two vertical dashed lines. It is 176.7 kHz for this case. 


The horizontal scale is 2 ys/div. The coil was driven at resonance for about 3 ms. The 
scope records the entire waveform. This particular screen is a section of that waveform at 
2.43 ms after initiation. 


When the variac is turned up to 40%, a streamer occurs when the toroid reaches a sufficient 
voltage. Fig. 8.2 shows the current building up slowly until streamer initiation at about 1.4 
ms and then a more rapid decline. The jagged appearance of the wave envelope is due to the 
scope’s algorithm to select representative points from the long waveform data set and does not 
imply that the current is experiencing rapid excursions. When one spreads out the waveform, 
it can be seen that the current looks as well behaved as the current in Fig. 8.1 except that it 
grows by a small amount each cycle. 


One cycle at 176.7 kHz lasts for about 5.66 ys. If the spark starts at 1400 ps, then it took 
1400/5.66 = 247 cycles to get to breakout. 


As charge leaves the toroid in the streamer, the toroid effectively gets a little larger due 
to the space charge in the streamer. This increases the capacitance and lowers the resonant 
frequency. The solid state drive must be able to slew with this change in frequency to keep 
driving the coil at its (changing) resonant frequency. My driver will maintain a fixed resonant 
frequency to within 50 Hz and then slew at 50-100 Hz/js to follow changes. These are 
nontrivial specifications, and efforts to develop a good controller have occupied considerable 
time over the past 15 years. 


Fig. 8.3 shows the voltage and current waveforms at 1.40 ms, just before breakout. They 
look the same as in Fig. 8.1 except larger. The voltage is now 604.2 V, the current is 17.1 A, 
and the instantaneous power (the average power at this instant) is about 9.3 kW. 


Fig. 8.4 shows the voltage and current waveforms during the streamer. Voltage has in- 
creased slightly, as might be expected during lower current flow. Current is holding steady at 
4.1 A. The controller is not holding the optimum phase angle between voltage and current, 
hence there is some ripple at the point of switching. If the phase angle were optimum, the 
current would be a little larger, perhaps in the 4.5 to 5 A range. Note that the frequency has 
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Figure 8.2: Current waveform at variac = 40% 


dropped only slightly, to 176.1 kHz. People with classical spark gap coils will notice a greater 
frequency shift because they have bigger sparks for smaller coils and toroids than my case. 


When the input voltage is increased, the current will build more rapidly, and streamer 
initiation is sooner. The current waveform for a variac setting of 60% is shown in Fig. 8.5. 
The streamer now starts at about 0.93 ms instead of 1.4 ms for the variac setting of 40%. Just 
before breakout the voltage is 913.5 V, and current is 21.5 A, and the instantaneous power is 


about 17.7 kW. 


The voltage and current during the spark for the variac set to 60% are shown in Fig. 8.6. 
The voltage has rebounded slightly from 913.5 to 936.7 V. The spark current is still about 4 
A. 


The coil appears to operate as a constant-current device in breakout. Both the input 
de current and the ac current into the coil stay about constant, perhaps even decreasing a 
little, as the variac is increased. The input impedance must increase proportional to the input 
voltage, in order for the current to remain constant. The reason for this is still obscure to me. 


The previous data were all collected in Kansas. The following data were collected in 
Colorado. I was surprised at how close the data matched between the two places. I installed 
coil 145 with the half spun toroid open side up on top, and the 30 inch aluminum sphere sitting 
on the flat ‘bottom’ of the toroid. The coil was raised by polyethylene spacers 1 inch above 
one of my copper ground planes (3 ft by 4 ft) resting on the concrete floor. The clearance 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 8—Sparks 8-6 


Obi too AZ {200%/ + 1.402 2.0027 Sngl4Ai STOP 


Figure 8.3: Voltage and current just before breakout at variac = 40% 


between the top of the aluminum sphere and the aluminum flashing on the ceiling is 36 inches. 


Some typical data for streamers to air are shown in Table 8.1. The time to firing will vary 
over a range. The longer the streamer takes to get started, the greater the input current and 
power at the time of firing. The numbers in the table are in the middle of the observed range. 
The first four rows of data were collected with a 1 inch diameter brass knob extending out 
the side of the sphere 0.5 inch. This provided a nice breakout point. The knob was removed 
for the last three rows. This increased the time to firing and also increased the variability of 
the data. 


8.2 Streamer Length 


Estimating length of a streamer ending in air is a challenge, more than sparks to a ground 
point where the distance can be checked later with a tape measure. My record on distance 
to a ground point is 143 cm, in front of about 60 members and guests of the local ham club. 
For streamers (discharges to air), I put a piece of Styrofoam behind a bump on the top load 
with a series of concentric circles drawn on it. I would get in a position to observe where 
the bump appeared to be in the center of the circles. I would set the controller in repetition 
mode (about one pulse per second) and watch the streamers. An average observed distance 
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Figure 8.4: Voltage and current during streamer at variac = 40% 


would be recorded and corrected for parallax. The corresponding power just before breakout 
would also be recorded. I found the streamer length to be related to the square root of power 
according to the formula 


é, =0.17VP inches (8.3) 


where P is measured in watts. John Freau has measured streamer lengths on his coils and 
determined the relationship should be @, = 1.7,/P, inches, where P, is the average wattmeter 
reading on the 60 Hz, 120/240 V side of his circuit. There is obviously a factor of 10 difference 
between the two formulas. I have watched his videos and am convinced his equation is valid. I 
think the difference is in the definition of power. Suppose John supplies 1 kW to his capacitor 
for 1 cycle of the 60 Hz waveform, then dumps all this stored energy to the secondary in 1/100 
of the 60 Hz waveform or 1/6000th of a second. The power during this impulse would then 
be 100 kW. Take the square root and we have the factor of 10. 


A classical disruptive coil operates at very high power levels when the spark gap fires. The 
capacitor may be charged to 15 kV and it is connected to the input of a step-up transformer 
which raises the effective applied voltage to the secondary even more. The toroid voltage rings 
up to its firing voltage in a few cycles (at most) of the coil resonant frequency, so the energy 
supplied to the capacitor over a relatively long time period will be dumped to the streamer 
quickly. Ionization of air requires some time, so the quicker the power is applied, the higher 
the voltage can rise, and the longer the streamer. 


Solid State Tesla Coil by Dr. Gary L. Johnson March 15, 2016 


Chapter 8—Sparks 8-8 


#- 0.005 2008/ 3ng14Ai STOP 


a -—t— SS SS 


¥YrmsCAi> chan off Yrms (AZ) =402 .6mv FreqcAl) chan off 


Figure 8.5: Current buildup for variac at 60% 


I have checked my formula to about 30 kW, with streamer lengths about 30 inches. Even 
modest classical coils can reach 54 inches of streamer, which would require 100 kW peak 
according to my formula. A coil with a 14 ft streamer (big by my standards) would require a 
peak power of 1 MW. Reaching these power levels at Tesla coil frequencies with silicon devices 
will be a challenge! 


8.3. Space Charge 


I was curious about the space charge near the Tesla coil top load during discharges, both for 
streamers to air and sparks to ground or other objects. I built a voltage divider with a 8 inch 
diameter copper sphere on top, as shown in Fig. 8.7. The copper sphere would be located 
within one to four feet of the Tesla coil. It was supported by a horizontal 42 inch piece of 
2 inch copper tubing, attached to a vertical section of PVC pipe with resistors inside. The 
top portion was composed of three physically large resistors (10.5 inches long by 1.125 inches 
outside diameter) of resistance 375 kQ each. The measured series resistance of the string was 
1.22 MQ. The bottom resistor of the divider was a 710 2 non inductive resistor. The large 
resistors were not marked non inductive, but did not appear to be wire wound (and therefore 
of high inductance). From various tests I concluded that the inductance of the 375 kQ resistors 
was not a major problem up to the operating frequency of the Tesla coil (in this case about 
158 kHz). 
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Figure 8.6: Voltage and current during spark for variac at 60% 


I wanted charge to flow through the divider chain rather than through a shielding system, 
so I did not try to shield the resistors. I put them inside a section of 8 inch PVC pipe, just 
for structural support. It turned out that the capacitance between the bodies of the 375 kQ 
resistors and the Tesla coil caused the divider chain to behave as if it were series resonant 
at the Tesla coil frequency. The effective impedance of the 1.22 MQ chain was substantially 
lowered from the dc value, making the output voltage higher than what would have been 
measured under ideal conditions by a factor of perhaps five to ten. The magnitude of the 
output divider voltage is therefore of little use in trying to determine actual voltage values. 


However, the variation of the divider output v,,, may be of use in a qualitative fashion 
to get a ‘feel’ for the space charge behavior. We will redraw Fig. 8.7 into something closer 
to a standard electrical schematic, shown in Fig. 8.8. The top load of the Tesla coil (a 30 
inch diameter aluminum sphere in this case) is represented by the larger circle at the top left 
of the schematic. The top of the voltage divider (a 8 inch copper sphere) is represented by 
the smaller circle at the top right. Between the two spheres is an effective capacitance. The 
magnitude of the capacitance can be estimated with the formulas given in the chapter on 
capacitance. Although its value is only a few pF, its impedance at the Tesla coil resonant 
frequency is low enough to allow a measurable current i,,, to flow in the divider circuit. This 
current flowing through the 710 Q resistor produces the voltage v,,, actually measured by the 
oscilloscope. 


As already stated, there are additional capacitances from all portions of the Tesla coil to 
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Table 8.1: Streamer Data 


Vro Itc P time knob 
Vrms Irms kW ms 
502. «14.1 = 7.0 1.26 yes 
742 14.2 94 0.82 yes 
9388 14.4 12.2 0.63 yes 
1118 14.5 144 0.51 yes 
734 20.3 13.4 1.50 no 
923 21.5 17.8 1.12 no 
1105 21.0 21.9 0.78 no 


Tesla 
coil 


Figure 8.7: Voltage Divider 


all portions of the divider circuit, which are not shown in this schematic. These do not change 
the arguments about the observed waveforms. 


The Tesla coil input current 7rc flows through a 0.024 2 resistor. The voltage v;r¢ across 
this resistor is measured by the oscilloscope. This current is associated with a resonant voltage 
rise through the Tesla coil. This voltage rise is represented by a voltage source Vjop. 


There are additional capacitances from the top load to ground that are not shown in this 
schematic. These additional paths to ground mean that 7,,,, will be substantially smaller than 
irc. But since these additional paths are similar in character to the one shown, it would seem 
reasonable that 75, is basically in phase with irc. Indeed, this is exactly what is measured, 
when allowance is made for the fact that virco and voy are 180 degrees out of phase (irc is 
entering the v;rc terminal while i.,4 is leaving the vou, terminal). 


Waveforms for Ugy¢ and v;rc before and during a spark are shown in Fig. 8.9. The upper 
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Figure 8.8: Voltage Divider Circuit 


waveform is Upy,z. The Invert button on the oscilloscope was not pushed for this plot so the two 
voltages appear out of phase. I will refer to the lower waveform as iro, although technically 
one would need to divide the mV values by 24 mQ to get the actual current. For this particular 
event there are two distinct portions of the discharge. First we see several cycles where Upyz 
has a strong offset above zero. During this portion of the discharge, irc decreases but at 
a relatively slow rate. Then something happens and v,,, abruptly returns to a waveform 
symmetric about the zero axis and starts a rapid decay toward zero. In the meantime, irc 
also starts a more rapid decline toward zero. 


What is a reasonable explanation for this plot? One possibility is the following: At some 
point in the voltage buildup of the Tesla coil, the electric field gets large enough to pull 
electrons from one or both of the Tesla coil top load and the divider top load. These electrons 
start to form an electron cloud between the two spheres. The spark between spheres does not 
occur until later. In the meantime, the electron cloud builds up in size and density. In this 
example the buildup takes perhaps 30 us. 


While the electron cloud is forming, there is a flow of electrons up through the 710 Q 
resistor, and also up through the 0.024 2 resistor. Electron flow up is the same as conventional 
current flow down, so Uy; will be positive. The charge in the electron cloud can come mostly 
from the voltage divider sphere or mostly from the Tesla coil top load, depending on local 
factors like the size and shape of the spark emitting points, so there is no requirement that the 
corresponding charge flows through the two routes be the same. If most of the electrons bound 
for the electron cloud flowed through the divider, then we would expect to see a substantial 
offset in Voyz and little offset in i7c. In looking at hundreds of traces similar to Fig. 8.9 I 
never saw a trace where irc had any significant offset. The current irc is several orders of 
magnitude larger than i,,,4, so even if the Tesla coil top load was contributing half the charge 
of the electron cloud, and the cloud was narrow in cross section and limited to the space 
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Figure 8.9: Top Waveform: vgyz, Bottom Waveform: irc. 


between spheres, the lack of electron flow in other directions might mean that any offset in 
irc would be difficult to see. Or it could mean that almost all the charge in the electron cloud 
was coming through the divider and not the Tesla coil. 


Once the electron cloud is fully formed, there is a sphere-to-sphere fault. Current flows 
directly from one conducting surface to the other, as opposed to just flowing into air. This 
clears out the electrons in the cloud, and lowers the electric field to where no new net charge 
flows into the space, at least from the divider. There may still be some charge flow into space 
from the Tesla coil top load along the fault path. The energy in the Tesla coil decays rapidly 
to zero, as does irc. If the Tesla coil is still sending out charge toward the divider, the divider 
‘sees’ a region of space at the same potential as the top load, but closer, so a given vo, will 
produce a larger Upyt. 


The variety in spark phenomena is quite large. The next three traces show some of the 
variety. Fig. 8.10 shows a greater offset than Fig. 8.9, a more rapid discharge, and a rebuilding 
of Vout but peaking at a lower level than before. Perhaps there are enough stray electrons 
hanging around for hundreds of ps to always cause discharges to occur at lower voltages than 
when the discharge channel has been off for hundreds of ms. 


Fig. 8.11 shows the offset in v,,,¢ but not an abrupt discharge between spheres. If the 
sphere-to-sphere fault does not occur, the Tesla coil voltage will still drop enough to inhibit 
additional electron flow into space. The offset will gradually and monotonically decrease to 
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Figure 8.10: voy large offset. 


Zero. 


Fig. 8.12 shows a case where the offset was much less pronounced, but there was still an 
abrupt return to symmetry about the axis, as though there was a sphere-to-sphere fault. The 
voltage rebuilds as before, but to a smaller terminal value. 


8.4 Conclusions 


It has been demonstrated that a Tesla coil secondary can be driven with IGBTs to produce 
streamers up to a meter in length. Streamers look the same as those from classical coils. 
Streamers get fatter the longer the coil is driven after streamer initiation. A streamer initiated 
at 1 ms and terminated at 2 ms will be much thinner and weaker looking than one terminated 
at 10 ms. 


Once the system is adjusted, it will produce several streamers per second for 24 hours 
per day and 7 days per week. No fan is necessary and there is no spark gap to get hot and 
corroded. 
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Figure 8.11: voy¢ small offset. 
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Figure 8.12: vo, no sphere-to-sphere discharge. 
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